Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,339,813 T→C 7.8% K85E (AAA→GAA)  mlaF ← ABC transporter maintaining OM lipid asymmetry, ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,339,8130TC 7.8% 364.6 / 15.8 115K85E (AAA→GAA) mlaFABC transporter maintaining OM lipid asymmetry, ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (56/50);  new base C (5/4);  total (61/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

Only 100 of 119 total aligned reads displayed.

AGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTG  >  NC_000913/3339665‑3339958
                                                                                                                                                    |                                                                                                                                                  
aGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGttt                                                                                                                                                    <  1:908380/149‑1 (MQ=255)
       cGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACTTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCAc                                                                                                                                              <  3:555299/148‑1 (MQ=255)
            tCCAGCTTCATTATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGtata                                                                                                                                        >  1:798990/1‑149 (MQ=255)
             ccAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGtata                                                                                                                                        >  4:337538/1‑148 (MQ=255)
             ccAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGtata                                                                                                                                        >  3:750371/1‑148 (MQ=255)
              cAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACa                                                                                                                                      >  1:461365/1‑149 (MQ=255)
               aGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACa                                                                                                                                      >  3:769213/1‑148 (MQ=255)
                 cTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑gc                                                                                                                                  <  2:821614/149‑1 (MQ=255)
                  ttCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑gcg                                                                                                                                 >  2:632228/1‑149 (MQ=255)
                  ttCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑gc                                                                                                                                  >  3:290491/1‑148 (MQ=255)
                     atcatcaCCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGa                                                                                                                              >  2:794987/1‑149 (MQ=255)
                      tcatcaCCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGa                                                                                                                              <  3:160070/148‑1 (MQ=255)
                       catcaCCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGaa                                                                                                                             <  3:527976/148‑1 (MQ=255)
                       catcaCCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGaa                                                                                                                             <  4:750371/148‑1 (MQ=255)
                         tcaCCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACga                                                                                                                          <  1:552741/149‑1 (MQ=255)
                         tcaCCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACg                                                                                                                           <  4:307227/148‑1 (MQ=255)
                           aCCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACgaga                                                                                                                        >  2:594100/1‑149 (MQ=255)
                             cGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGAc                                                                                                                       >  3:364453/1‑148 (MQ=255)
                                 aaaTGCCACAATGGCGCCGGCAGTTTTGTATGTTCGCGCAATCGATAGGCGGCGATGTCAAAAACTTTCATAGCAGTGAACAACGCCCCGGACTGAAAAATCATTGTCATCCGTTTGCGCACTGTATCCCG‑GCGCGAACGCGACACCgt                                                                                                                  <  2:502132/147‑2 (MQ=255)
                                   aTGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGcc                                                                                                                 <  4:729432/148‑1 (MQ=255)
                                       aacaaTGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGaa                                                                                                             >  2:386495/1‑148 (MQ=255)
                                        acaaTGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAat                                                                                                            <  4:107821/148‑1 (MQ=255)
                                            tGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATtct                                                                                                       <  2:662059/149‑1 (MQ=255)
                                             ggCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATtct                                                                                                       <  3:961426/148‑1 (MQ=255)
                                              gcgcGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATtctc                                                                                                      <  3:836768/148‑1 (MQ=255)
                                                        ttGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTCTACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGaaaa                                                                                           <  2:913327/149‑1 (MQ=255)
                                                        ttGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGaaaa                                                                                           <  1:616882/149‑1 (MQ=255)
                                                         tGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGaaaa                                                                                           <  3:122331/148‑1 (MQ=255)
                                                           ggTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGa                                                                                        >  1:894143/1‑149 (MQ=255)
                                                            gTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGAt                                                                                       <  1:843812/149‑1 (MQ=255)
                                                                 ttCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCa                                                                                   >  3:353361/1‑148 (MQ=255)
                                                                  tCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCAc                                                                                  >  4:1009891/1‑148 (MQ=255)
                                                                      gcAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGa                                                                             <  2:1009781/149‑1 (MQ=255)
                                                                        aGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGAt                                                                            <  3:75921/148‑1 (MQ=255)
                                                                           ggATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTgg                                                                        <  1:501761/149‑1 (MQ=255)
                                                                           ggATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTgg                                                                        <  2:1060175/149‑1 (MQ=255)
                                                                            gATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGt                                                                       >  1:203546/1‑149 (MQ=255)
                                                                             aTAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGt                                                                       >  3:1008968/1‑148 (MQ=255)
                                                                             aTAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTg                                                                      <  2:882768/149‑1 (MQ=255)
                                                                              tAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGc                                                                     >  1:248820/1‑149 (MQ=255)
                                                                                 gCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGa                                                                   <  4:409746/148‑1 (MQ=255)
                                                                                 gCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGAt                                                                  <  2:391106/149‑1 (MQ=255)
                                                                                  cGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGAt                                                                  >  4:671415/1‑148 (MQ=255)
                                                                                       ttGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGccc                                                            <  2:1018985/149‑1 (MQ=255)
                                                                                          tCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGc                                                          <  3:84738/148‑1 (MQ=255)
                                                                                           cAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGcc                                                         >  3:328529/1‑148 (MQ=255)
                                                                                              aTACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAAt                                                      >  4:14725/1‑148 (MQ=255)
                                                                                                aCGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCa                                                    >  3:430182/1‑148 (MQ=255)
                                                                                                 cGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGa                                                  <  2:38559/149‑1 (MQ=255)
                                                                                                  gTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCg                                                                                                                >  2:458753/1‑86 (MQ=255)
                                                                                                  gTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCg                                                                                                                <  1:458753/86‑1 (MQ=255)
                                                                                                       tatCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTCGCGC‑CTGTATACAGTGCGCAAACgga                                                                                                                         <  1:1063353/72‑3 (MQ=21)
                                                                                                       tatCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTCGCGC‑CTGTATACAGTGCGCAAACgga                                                                                                                         >  2:1063353/1‑70 (MQ=21)
                                                                                                         tCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAg                                          >  2:629448/1‑149 (MQ=255)
                                                                                                           aGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGcg                                        >  1:809051/1‑149 (MQ=255)
                                                                                                            gTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGcgt                                       >  1:479632/1‑149 (MQ=255)
                                                                                                             tGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGcgtc                                      >  1:141302/1‑149 (MQ=255)
                                                                                                                acaacGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGcgtcgt                                    >  4:463290/1‑148 (MQ=255)
                                                                                                                acaacGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGtt                                   <  2:479632/149‑1 (MQ=255)
                                                                                                                 caacGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGtt                                   <  3:1009891/148‑1 (MQ=255)
                                                                                                                 caacGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGtt                                   <  2:717593/148‑1 (MQ=255)
                                                                                                                 caacGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGtt                                   >  1:717593/1‑148 (MQ=255)
                                                                                                                  aacGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGtttt                                 >  2:829078/1‑149 (MQ=255)
                                                                                                                   acGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGtttt                                 >  3:744473/1‑148 (MQ=255)
                                                                                                                    cGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTa                                <  4:1008968/148‑1 (MQ=255)
                                                                                                                    cGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTa                                <  4:766289/148‑1 (MQ=255)
                                                                                                                       cccGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCa                                                                                >  1:176963/1‑97 (MQ=255)
                                                                                                                       cccGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGa                            <  2:980714/149‑1 (MQ=255)
                                                                                                                        ccGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGa                            >  3:84943/1‑148 (MQ=255)
                                                                                                                        ccGGACCGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGAt                           >  2:679444/1‑149 (MQ=255)
                                                                                                                            aCTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGccc                       >  2:231868/1‑149 (MQ=255)
                                                                                                                            aCTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGcc                        <  4:438374/148‑1 (MQ=255)
                                                                                                                             cTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGccc                       >  3:698445/1‑148 (MQ=255)
                                                                                                                               gAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCGGCCAATCAGACGGATTAGCGTCGTTTTACCGATGCCCGa                     >  3:697443/1‑148 (MQ=255)
                                                                                                                               gAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGa                            <  4:795612/141‑1 (MQ=255)
                                                                                                                                aaaTAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGAt                    >  4:704898/1‑148 (MQ=255)
                                                                                                                                 aaTAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATgg                  >  1:537025/1‑149 (MQ=255)
                                                                                                                                 aaTAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATgg                  >  1:46830/1‑149 (MQ=255)
                                                                                                                                  aTAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATgg                  >  4:718211/1‑148 (MQ=255)
                                                                                                                                        tGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCAt            >  3:731905/1‑148 (MQ=255)
                                                                                                                                          cTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGa          >  3:290497/1‑148 (MQ=255)
                                                                                                                                           tCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATc        >  2:313474/1‑149 (MQ=255)
                                                                                                                                             aTCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGc      >  1:946157/1‑149 (MQ=255)
                                                                                                                                              tCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCt                                                                          >  4:278393/1‑80 (MQ=255)
                                                                                                                                              tCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCt                                                                          <  3:278393/80‑1 (MQ=255)
                                                                                                                                              tCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGcc     <  1:679444/149‑1 (MQ=255)
                                                                                                                                              tCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGc      >  3:63202/1‑148 (MQ=255)
                                                                                                                                              ctcGTTCGCGC‑CTGTATACAGTGCGCAAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGcc     >  1:383886/3‑149 (MQ=255)
                                                                                                                                               tcGTTCGCGC‑CTGTATACAGTGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGa          <  2:568613/142‑1 (MQ=255)
                                                                                                                                               ccGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCg    >  2:103086/1‑149 (MQ=255)
                                                                                                                                                cGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCg    <  3:463290/148‑1 (MQ=255)
                                                                                                                                                cGTTCGCGC‑CTGTATACAGTGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGcgtcg                                     >  4:25289/1‑115 (MQ=255)
                                                                                                                                                cGTTCGCGC‑CTGTATACAGTGCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGcgtcg                                     <  3:25289/115‑1 (MQ=255)
                                                                                                                                                cGTTCGCGC‑CTGTATACAGTGCGCAAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCg    >  3:601156/1‑148 (MQ=255)
                                                                                                                                                 gTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGcgt                                       <  3:79013/112‑1 (MQ=255)
                                                                                                                                                 gTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGcgt                                       >  4:79013/1‑112 (MQ=255)
                                                                                                                                                 gTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTg  >  2:908851/1‑149 (MQ=255)
                                                                                                                                                  tttGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGa                            <  3:902739/122‑1 (MQ=255)
                                                                                                                                                  tttGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGa                            >  4:902739/1‑122 (MQ=255)
                                                                                                                                                  tttGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTg  >  4:622655/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
AGCCCCACGGCCTCCAGCTTCATCATCACCGTACTATGCAACAATGGCGCGGGAAGTTGGGTATGTTCGCGCAGTGGATAGGCGACGTTGTCAAATACGTTCATATCAGTGAACAACGCCCCGGACTGAAATAACATGCTCATCCGTTTGCGCACTGTATACAG‑GCGCGAACGAGACATCGCCGGAATATTCTCACCATCGAAAAGGATCTCACCATGATCTGGTGCGATTTGCCCGCCAATCAGACGGAGTAGCGTCGTTTTACCGATGCCCGATGGCCCCATGATCGCCGTG  >  NC_000913/3339665‑3339958

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