Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,438,586 C→T 6.6% G733G (GGG→GGA)  pfo ← pyruvate‑flavodoxin oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,438,5860CT 6.6% 164.4 / 3.0 61G733G (GGG→GGA) pfopyruvate‑flavodoxin oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (25/32);  new base T (2/2);  total (27/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.75e-01

TGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTGCGCGAATAGCTGAGTGTGGGCAAGCGGCAACGCAGTGGTTACATTGGGTACAGAGTTCCTCTTTCC  >  NC_000913/1438438‑1438734
                                                                                                                                                    |                                                                                                                                                    
tGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGc                                                                                                                                                      <  2:991437/149‑1 (MQ=255)
           tcttcGACATGTTCCAGGCGAGACGTCATATTGACGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGTCCGCGCATATc                                                                                                                                           >  1:657861/1‑149 (MQ=255)
           tcttcGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGTCCGCGCATATc                                                                                                                                           >  1:1021019/1‑149 (MQ=255)
           tcttcGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATc                                                                                                                                           <  1:284454/149‑1 (MQ=255)
                 aCATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCAGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATAGTGCTGCCCGCGCATGTCAcgcgc                                                                                                                                     >  2:201211/1‑148 (MQ=255)
                      ttCCACGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCCTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCCCGCGATTTca                                                                                                                                <  2:695890/149‑1 (MQ=255)
                        ccAGGCGAGAAATCATATTGAGGGCTTTAATCTCTGGATTCTGACGGTATTTCGCCGGGCAAACTTCAACGCACAGGTTAGAACCGGTGCAATCTTCCGGTGCAACCTGCACGACATATTTCTGCCAGCGCAGAACACGCGAGTTgaca                                                                                                                              >  2:326555/1‑149 (MQ=255)
                          aGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCGGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATc                                                                                                                            >  2:914403/1‑149 (MQ=255)
                                 aCATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACCGGTTACACCCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGaa                                                                                                                     <  2:47452/149‑1 (MQ=255)
                                   atcatATTGATGGCTTTAATCTCTGGCTTCAGACGGTCTTTCGCCGGGCAAACTTCGACGCACATGTTACAACCGGTGCAATCTTCCGGTGCCGCCTGCAAGCCATATTTCTGCCCGCGCATCTCACGCGATTTCCCATCCCGCGACTg                                                                                                                   >  2:1143132/1‑149 (MQ=255)
                                    tcatATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACCTCCTGCGAGTGc                                                                                                                  >  1:519500/1‑149 (MQ=255)
                                              gccTTTAATCCCTCGACTCTGACGCACTTTCCCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAg                                                                                                        <  1:592341/147‑1 (MQ=255)
                                                  ttAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACCTCCAGCGAATGCAGGCTGGCCGGGGc                                                                                                    >  1:1015812/1‑149 (MQ=255)
                                                  ttAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCCGCGAATGCAGGCTCGCAGGGGc                                                                                                    >  1:873984/1‑149 (MQ=255)
                                                  ttAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGc                                                                                                    >  1:706293/1‑149 (MQ=255)
                                                   tcagcgctGCATTCTGCCGGTCCTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCg                                                                                                   <  2:508597/143‑1 (MQ=255)
                                                   tAATCTCGTGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCg                                                                                                   <  2:296944/149‑1 (MQ=255)
                                                   gaagctctGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCg                                                                                                   <  1:451682/145‑1 (MQ=255)
                                                        tctGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTc                                                                                              >  2:804069/1‑149 (MQ=255)
                                                         cgggATTCTGCCGGTCTTTGGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTcc                                                                                             <  2:240570/147‑1 (MQ=255)
                                                           ggATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCAt                                                                                           <  1:1045237/149‑1 (MQ=255)
                                                            gATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATc                                                                                          >  2:1199594/1‑149 (MQ=255)
                                                              ttCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGc                                                                                        <  2:609634/149‑1 (MQ=255)
                                                                 tGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCGTCGCTTc                                                                                     >  2:636601/1‑149 (MQ=255)
                                                                  gACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCACGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCa                                                                                    >  2:1165966/1‑149 (MQ=255)
                                                                  gACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCa                                                                                    >  1:183768/1‑149 (MQ=255)
                                                                  gACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTCCCATCGCTTCa                                                                                    >  1:183769/1‑149 (MQ=255)
                                                                  gACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAAAGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCa                                                                                    >  2:726870/1‑149 (MQ=255)
                                                                    cGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAgg                                                                                  <  2:1033840/149‑1 (MQ=255)
                                                                        cgttCGCCGTGCAAACTTCGACGCACAGGTTACAGCCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCTCGCGATTTCCCATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAggcggc                                                                              <  2:809877/147‑1 (MQ=255)
                                                                         tttCGCCGGGCACACTTCGACGAACAGGTTACACCCGGTGCAATCTTCCGGTGCCACCGGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGca                                                                             <  2:997517/149‑1 (MQ=255)
                                                                         gttCGCCGCGCAAACTTCGACGCCCAGGTTACAACCGGTGCAATCTTCCGGGGCCACCTGGAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACAGCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGca                                                                             <  2:29148/148‑1 (MQ=255)
                                                                              ccGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGGTTTCCATCGCTTCAGGCGGCACCACt                                                                        >  1:688244/1‑149 (MQ=255)
                                                                              ccGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATCTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCCCt                                                                        >  2:1281253/1‑149 (MQ=255)
                                                                              ccGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGAGTTCACATCCAGCGACTGCAGGCTGGCAGGGGCGTTTTCCATCCCTTCCGGCGGCCCCACt                                                                        >  1:768320/1‑149 (MQ=255)
                                                                                 ggCAAACTTCGTCGACCAGGTAACAACCGGTGCAATCCTCCGGGGCCACCTGCAAGTCATTAGTCTGCCCGCGCTTATCACGCTTTTTCACATGCAGCGTAGGCGGTTTTGCAGGGGCTGTTGCCATCGCTTCAGGCGGCACCACttgt                                                                     <  2:49596/149‑3 (MQ=255)
                                                                                 ggCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGTCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACtttt                                                                     <  2:657861/149‑1 (MQ=255)
                                                                                 ggCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGTCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACtttt                                                                     <  2:1021019/149‑1 (MQ=255)
                                                                                        ttCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTGCGCGaa                                                              <  1:420999/149‑1 (MQ=255)
                                                                                        ttCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTGCGCGaa                                                              <  2:682468/149‑1 (MQ=255)
                                                                                              gCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTGCGCGAATAGCTg                                                        >  2:1124109/1‑149 (MQ=255)
                                                                                                    gTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTGCGCGAATAGCTGAGtgtg                                                  <  1:1237119/149‑1 (MQ=255)
                                                                                                     ttACCACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTGCGCGAATAGCTGAGTGTgg                                                 <  1:151574/149‑1 (MQ=255)
                                                                                                     ttACAACCGGGGGAATCTCCCGGTGCCACCTGCAGGCCACTTGTCGGCTCGGGCATATCACGCGATTCTACATCTATTGAATGCCGGTTGGTAGCGGCATTTTTCCTCTTTGCAGGTGTCGCCACTTTTGCGCGAATAGCTGAGTGTgg                                                 <  2:376556/149‑1 (MQ=255)
                                                                                                           ctgggtcAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACTTCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTGCGCGAATAGCTGAGTGTGGGCAAGc                                           <  2:966182/143‑1 (MQ=255)
                                                                                                              gTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGACTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTGCGCGAATAGCTGAGTGTGGGCAAGCGGc                                        <  1:653517/149‑1 (MQ=255)
                                                                                                                  aaTCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTGCGCGAATAGCTGAGGGTGGGCAAGCGGCAACg                                    >  1:387252/1‑149 (MQ=255)
                                                                                                                     cTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACAGCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTGCGCGAATAGCTGAGTGTGGGCAAGCGGCAACGCAg                                 >  1:690399/1‑149 (MQ=255)
                                                                                                                       tCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTGCGCGAATAGCTGAGTGTGGGCAAGCGGCAACGCAGTg                               <  1:294902/149‑1 (MQ=255)
                                                                                                                          ggTGCCACCTGCAAGCCATATTTCTGCCCGCGCATATCCCGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTGCGCGAATAGCTGAGTGTGGGCAAGCGGCAACGCAGTGGtt                            <  1:548432/149‑1 (MQ=255)
                                                                                                                            tGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTGCGCGAATAGCTGAGTGTGGGCAAGCGGCAACGCAGTGGTTAc                          <  2:319516/149‑1 (MQ=255)
                                                                                                                                aCCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTGCGCGAATAGCTGAGTGTGGGCAAGGGGCAACGCAGTGGTTACATTg                      >  1:145154/1‑149 (MQ=255)
                                                                                                                                aCCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACAACCAGAGAAGTCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTGCGCGAATAGCTGAGTGTGGGCAAGCGGCAACGCAGTGGTTACATTg                      >  1:1013757/1‑149 (MQ=255)
                                                                                                                                aCCTGCAAGACAGATTTCTGCCCGCGCATATCACGCGATATCACATCCAGCGAATGCAGGCTGGCAGGGGGGTTTTCCATCGCTTCAGGCGGCACCCCTTTTGCACGCATAGCTGACTGCGGGCAACCGCCAACGCAAGGGTTCCaacg                      >  2:130288/1‑146 (MQ=255)
                                                                                                                                 ccTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTGCGCGAATAGCTGAGTGTGGGCAAGCGGCAACGCAGTGGTTACATTgg                     <  1:804069/149‑1 (MQ=255)
                                                                                                                                      aaGACATATTTCTGCCCGCGCATATCACTCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTGCCCGAAGAGCTGAGTGTGGGCAAGCGGCAACGCAGTGGTTACATTGGGTACa                <  1:256447/149‑1 (MQ=255)
                                                                                                                                      aaGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTGCGCGAATAGCTGGGTGTGGGCAAGCGGCAACGCAGTGGTTACATTGGGTACa                <  2:440817/149‑1 (MQ=255)
                                                                                                                                      aaGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACCCATGCGCGAATAGCTGAGTGTGGGCAAGCGGCAACGCAGTGGTTACATTGGGTACa                <  1:894240/149‑1 (MQ=255)
                                                                                                                                      aaGACATATGTCTGCCCGCGCATATCGCGCGATTTCAGATCCAGGTAATGCGGGCGGTCAGGGGCGTTTTCCACGGCTTGAGGCGGGAGCACGTTTGCGCGAATAGCGGCGTGTGGGGAAGCGGCAACGCAGTGGTTACCTTGGGGACa                <  2:554149/149‑1 (MQ=255)
                                                                                                                                           atatTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTGCGCGAATAGCTGAGTGTGGGCAAGCGGCAACGCAGTGGTTACATTGGGTACAGGGtt           >  2:694495/1‑149 (MQ=255)
                                                                                                                                                   gCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTGCGCGAATAGCTGAGTGTGGGCAAGCGGCAACGCAGTGGTTACATTGGGTACAGAGTTCCTCTTTc   <  2:232181/149‑1 (MQ=255)
                                                                                                                                                   gCCCCCGTATATCACGCGCTTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACCTTTGCGCGAATAGCTGAGTGTGCGCAAGCGGCAACGCAGTGGGTACATTGGGTACAGAGTTCCTCTTTc   <  1:1143132/149‑1 (MQ=255)
                                                                                                                                                    cccGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTGCGCGAATAGCTGAGTGTGGGCAAGCGGCAACGCCGTGATTAATTCAGGTACAGAGTTCCTCGTTcc  >  2:648137/1‑149 (MQ=255)
                                                                                                                                                    cccGCGCATATCACGCGATTTAACATCCAGAGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTGCGCGAATAGCTGAGTGTGGGCAAGCGGCAACGCAGTGGTTACATTGGGTACAGAGTTCCTCATTcc  >  1:680806/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TGATTTTCTCTTCTTCGACATGTTCCAGGCGAGACATCATATTGATGGCTTTAATCTCTGGATTCTGACGGTCTTTCGCCGGGCAAACTTCGACGCACAGGTTACAACCGGTGCAATCTTCCGGTGCCACCTGCAAGACATATTTCTGCCCGCGCATATCACGCGATTTCACATCCAGCGAATGCAGGCTGGCAGGGGCGTTTTCCATCGCTTCAGGCGGCACCACTTTTGCGCGAATAGCTGAGTGTGGGCAAGCGGCAACGCAGTGGTTACATTGGGTACAGAGTTCCTCTTTCC  >  NC_000913/1438438‑1438734

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑