Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,660,199 T→C 5.4% M31V (ATG→GTG)  iscU ← iron‑sulfur cluster assembly scaffold protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,660,1990TC 5.4% 236.6 / 3.1 74M31V (ATG→GTG) iscUiron‑sulfur cluster assembly scaffold protein
Reads supporting (aligned to +/- strand):  ref base T (36/34);  new base C (2/2);  total (38/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.57e-01

TTCCCTTTCACCCATTCGGTGACCAGGGAGCTGGAAGCGATAGCGGAACCGCAGCCGTAAGTTTTAAAACGCGCGTCTTCAATGATACCTTCATCGTTGACTTTAATCTGCAACTTCATCACGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTAATGATGAGCCCATTCGATGCTGTTCA  >  NC_000913/2660051‑2660345
                                                                                                                                                    |                                                                                                                                                  
ttCCCTTTCACCCATTCGGTGACCAGGGAGCTGGAAGCGATAGCGGAACCGCAGCCGTAAGTTTTAAAACGCGCGTCTTCAATGATACCTTCATCGTTGACTTTAATCTGCAACTTCATCACGTCGCCACAGGCCGGTGCCCCCACCAt                                                                                                                                                    >  2:392266/1‑149 (MQ=255)
     tttCACCCATTCGGTGACCAGGGAGCTGGAAGCGATAGCGGAACCGCAGCCGTAAGTTTTAAAACGCGCGTCTTCAATGATACCTTCATCGTTGACTTTAATCTGCAACTTCATCACGTCGCCACAGGCCGGTGCCCCCACCATGCCGc                                                                                                                                               >  8:93197/1‑149 (MQ=255)
     tttCACCCATTCGGTGACCAGGGAGCTGGAAGCGATAGCGGAACCGCAGCCGTAAGTTTTAAAACGCGCGTCTTCAATGATACCTTCATCGTTGACTTTAATCTGCAACTTCATCACGTCGCCACAGGCCGGTGCCCCCACCATGCCGc                                                                                                                                               >  8:307857/1‑149 (MQ=255)
         aCCCATTCGGTGACCAGGGAGCTGGAAGCGATAGCGGAACCGCAGCCGTAAGTTTTAAAACGCGCGTCTTCAATGATACCTTCATCGTTGACTTTAATCTGCAACTTCATCACGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGc                                                                                                                                            <  6:320526/148‑1 (MQ=255)
           ccATTCGGTGACCAGGGAGCTGGAAGCGATAGCGGAACCGCAGCCGTAAGTTTTAAAACGCGCGTCTTCAATGATACCTTCATCGTTGACTTTAATCTGCAACTTCATCACGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGa                                                                                                                                         >  8:316458/1‑149 (MQ=255)
              ttCGGTGACCAGGGAGCTGGAAGCGATAGCGGAACCGCAGCCGTAAGTTTTAAAACGCGCGTCTTCAATGATACCTTCATCGTTGACTTTAATCTGCAACTTCATCACGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGACGt                                                                                                                                      >  3:152585/1‑149 (MQ=255)
                  gTGACCAGGGAGCTGGAAGCGATAGCGGAACCGCAGCCGTAAGTTTTAAAACGCGCGTCTTCAATGATACCTTCATCGTTGACTTTAATCTGCAACTTCATCACGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGAc                                                                                                                                        <  4:91516/143‑1 (MQ=255)
                  gTGACCAGGGAGCTGGAAGCGATAGCGGAACCGCAGCCGTAAGTTTTAAAACGCGCGTCTTCAATGATACCTTCATCGTTGACTTTAATCTGCAACTTCATCACGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGAc                                                                                                                                        >  3:91516/1‑143 (MQ=255)
                  gTGACCAGGGAGCTGGAAGCGATAGCGGAACCGCAGCCGTAAGTTTTAAAACGCGCGTCTTCAATGATACCTTCATCGTTGACTTTAATCTGCAACTTCATCACGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGACGTtctc                                                                                                                                  <  7:307857/149‑1 (MQ=255)
                          ggAGCTGGAAGCGATAGCGGAACCGCAGCCGTAAGTTTTAAAACGCGCGTCTTCAATGATACCTTCATCGTTGACTTTAATCTGCAACTTCATCACGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCgttgt                                                                                                                          >  4:112150/1‑149 (MQ=255)
                           gAGCTGGAAGCGATAGCGGAACCGCAGCCGTAAGTTTTAAAACGCGCGTCTTCAATGATACCTTCATCGTTGACTTTAATCTGCAACTTCATCACGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCgttgtt                                                                                                                         <  5:262346/149‑1 (MQ=255)
                               tGGAAGCGATAGCGGAACCGCAGCCGTAAGTTTTAAAACGCGCGTCTTCAATGATACCTTCATCGTTGACTTTAATCTGCAACTTCATCACGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTc                                                                                                                      >  7:135316/1‑148 (MQ=255)
                                 gAAGCGATAGCGGAACCGCAGCCGTAAGTTTTAAAACGCGCGTCTTCAATGATACCTTCATCGTTGACTTTAATCTGCAACTTCATCACGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCa                                                                                                                     >  4:296308/1‑147 (MQ=255)
                                 gAAGCGATAGCGGAACCGCAGCCGTAAGTTTTAAAACGCGCGTCTTCAATGATACCTTCATCGTTGACTTTAATCTGCAACTTCATCACGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTTTGTCaaa                                                                                                                   <  7:259825/148‑1 (MQ=255)
                                     cGATAGCGGAACCGCAGCCGTAAGTTTTAAAACGCGCGTCTTCAATGATACCTTCATCGTTGACTTTAATCTGCAACTTCATCACGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGaa                                                                                                               <  1:392266/149‑1 (MQ=255)
                                       aTAGCGGAACCGCAGCCGTAAGTTTTAAAACGCGCGTCTTCAATGATACCTTCATCGTTGACTTTAATCTGCAACTTCATCACGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAAcc                                                                                                             >  8:88322/1‑149 (MQ=255)
                                              aaCCGCAGCCGTAAGTTTTAAAACGCGCGTCTTCAATGATACCTTCATCGTTGACTTTAATCTGCAACTTCATCACGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTa                                                                                                      >  7:372640/1‑149 (MQ=255)
                                                 cGCAGCCATAAGTTTTAAAACGCGCGTCTTCAATGATACCTTCATCGTTGACTTTAATCTGCAACTTCATCACGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTAcgc                                                                                                   <  8:246249/149‑1 (MQ=255)
                                                    aGCCGTAAGTTTTAAAACGCGCGTCTTCAATGATACCTTCATCGTTGACTTTAATCTGCAACTTCATCACGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGa                                                                                                <  2:15232/149‑1 (MQ=255)
                                                       cGTAAGTTTTAAAACGCGCGTCTTCAATGATACCTTCATCGTTGACTTTAATCTGCAACTTCATCACGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATtc                                                                                             <  8:35288/149‑1 (MQ=255)
                                                         tAAGTTTTAAAACGCGCGTCTTCAATGATACCTTCATCGTTGACTTTAATCTGCAACTTCATCACGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATtctc                                                                                           <  3:28407/149‑1 (MQ=255)
                                                             ttttAAAACGCGCGTCTTCAATGATACCTTCATCGTTGACTTTAATCTGCAACTTCATCACGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAAccc                                                                                                            <  1:244484/128‑1 (MQ=255)
                                                             ttttAAAACGCGCGTCTTCAATGATACCTTCATCGTTGACTTTAATCTGCAACTTCATCACGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAAccc                                                                                                            >  2:244484/1‑128 (MQ=255)
                                                             ttttAAAACGCGCGTCTTCAATGATACCTTCATCGTTGACTTTAATCTGCAACTTCATCACGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATtctc                                                                                           <  1:248855/145‑1 (MQ=255)
                                                             ttttAAAACGCGCGTCTTCAATGATACCTTCATCGTTGACTTTAATCTGCAACTTCATCACGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATtctc                                                                                           >  2:248855/1‑145 (MQ=255)
                                                               ttAAAACGCGCGTCTTCAATGATACCTTCATCGTTGACTTTAATCTGCAACTTCATCACGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATg                                                                                     >  3:178420/1‑149 (MQ=255)
                                                               ttAAAACGCGCGTCTTCAATGATACCTTCATCGTTGACTTTAATCTGCAACTTCATCACGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATg                                                                                     >  1:381846/1‑149 (MQ=255)
                                                                  aaaCGCGCGTCTTCAATGATACCTTCATCGTTGACTTTAATCTGCAACTTCATCACGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGt                                                                                   <  3:112150/148‑1 (MQ=255)
                                                                  aaaCGCGCGTCTTCAATGATACCTTCATCGTTGACTTTAATCTGCAACTTCATCACGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGt                                                                                   <  1:14412/148‑1 (MQ=255)
                                                                  aaaCGCGCGTCTTCAATGATACCTTCATCGTTGACTTTAATCTGCAACTTCATCACGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGt                                                                                   <  7:316458/148‑1 (MQ=255)
                                                                    aCGCGCGTCTTCAATGATACCTTCATCGTTGACTTTAATCTGCAACTTCATCACGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTAGTTGTCAAAGGAACCCACCTTACGCGGATTCTCGTAATGTTCg                                                                                 >  5:374055/1‑148 (MQ=255)
                                                                      gcgcgTCTTCAATGATACCTTCATCGTTGACTTTAATCTGCAACTTCATCACGTCGCCACAGGCCGGAGCCCCCACCATGCCGCTGCCGCCGGTCTCGTCGGTGTGGTCAACGGCACCCCCGTTCCGAGGCGTCTCGTAAGTGTCGAt                                                                               >  8:286099/1‑148 (MQ=255)
                                                                           tCTTCAATGATACCTTCATCGTTGACTTTAATCTGCAACTTCATCACGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACttt                                                                         <  7:88322/149‑1 (MQ=255)
                                                                             ttCAATGATACCTTCATCGTTGACTTTAATCTGCAACTTCATCACGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATCGTCGATAACttttt                                                                       >  3:305648/1‑149 (MQ=255)
                                                                               cAATGATACCTTCATCGTTGACTTTAATCTGCAACTTCATCACGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACttttt                                                                       <  2:381846/147‑1 (MQ=255)
                                                                               cAATGATACCTTCATCGTTGACTTTAATCTGCAACTTCATCACGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCg                                                                     >  7:362012/1‑149 (MQ=255)
                                                                                 aTGATACCTTCATCGTTGACTTTAATCTGCAACTTCATCACGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATgg                                                                                    >  5:18754/1‑132 (MQ=255)
                                                                                 aTGATACCTTCATCGTTGACTTTAATCTGCAACTTCATCACGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATgg                                                                                    <  6:18754/132‑1 (MQ=255)
                                                                                  tGATACCTTCATCGTTGACTTTAATCTGCAACTTCATCACGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTg                                                                  >  5:152135/1‑149 (MQ=255)
                                                                                     tACCTTCATCGTTGACTTTAATCTGCAACTTCATCACGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTaa                                                               <  4:178420/149‑1 (MQ=255)
                                                                                        cTTCATCGTTGACTTTAATCTGCAACTTCATCACGTCGCCACAGGCCGGTGCCCCCACCACGCCGCTGCCGACGTTCTCGTCGTTGTTGTCaaa                                                                                                                   >  5:273666/1‑94 (MQ=255)
                                                                                        cTTCATCGTTGACTTTAATCTGCAACTTCATCACGTCGCCACAGGCCGGTGCCCCCACCACGCCGCTGCCGACGTTCTCGTCGTTGTTGTCaaa                                                                                                                   <  6:273668/94‑1 (MQ=255)
                                                                                        cTTCATCGTTGACTTTAATCTGCAACTTCATCACGTCGCCACAGGCCGGTGCCCCCACCACGCCGCTGCCGACGTTCTCGTCGTTGTTGTCaaa                                                                                                                   <  6:273666/94‑1 (MQ=255)
                                                                                        cTTCATCGTTGACTTTAATCTGCAACTTCATCACGTCGCCACAGGCCGGTGCCCCCACCACGCCGCTGCCGACGTTCTCGTCGTTGTTGTCaaa                                                                                                                   >  5:273668/1‑94 (MQ=255)
                                                                                         ttCATCGTTGACTTTAATCTGCAACTTCATCACGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGt                                                                 <  2:67897/143‑1 (MQ=255)
                                                                                         ttCATCGTTGACTTTAATCTGCAACTTCATCACGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGt                                                                 >  1:67897/1‑143 (MQ=255)
                                                                                            aTCGTTGACTTTAATCTGCAACTTCATCACGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATta                                                        <  3:322843/149‑1 (MQ=255)
                                                                                               gTTGACTTTAATCTGCAACTTCCTCACGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATaa                                                     <  6:260534/149‑1 (MQ=255)
                                                                                               gTTGACTTTAATCTGCAACTTCATCACGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATtata                                                      >  2:191291/1‑148 (MQ=255)
                                                                                               gTTGACTTTAATCTGCAACTTCATCACGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGATGTCAAAGGAACCCACGATACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATtata                                                      >  8:121423/1‑148 (MQ=255)
                                                                                                ttGACTTTAATCTGCAACTTCATCACGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATaaa                                                    >  5:155381/1‑149 (MQ=255)
                                                                                                   aCTTTAATCTGCAACTTCATCACGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAAtt                                                  <  1:25512/148‑1 (MQ=255)
                                                                                                   aCTTTAATCTGCAACTTCATCACGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATtc                                                 >  4:349189/1‑149 (MQ=255)
                                                                                                     tttAATCTGCAACTTCATCACGTCGCCACAGGCCGGTGCCCCCCCCATGCCGCTGCCGACGGTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATtctc                                               >  5:16621/1‑149 (MQ=255)
                                                                                                     tttAATCTGCAACTTCATCACGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATtctc                                               >  7:411608/1‑149 (MQ=255)
                                                                                                         aTCTGCAACTTCATCACGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGa                                           <  7:93197/149‑1 (MQ=255)
                                                                                                         aTCTGCAACTTCATCACGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGa                                           >  7:386590/1‑149 (MQ=255)
                                                                                                             gCAACTTCATCACGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTcc                                       >  6:29459/1‑149 (MQ=255)
                                                                                                               aaCTTCATCACGTCGCCTCAGGCCGGCGGCCCCAACTTGCCGCTGGCTACGTTCTCGTCGTTGTTGTCCTAGGAACCCAGGTTACGCGGATTCTCTTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAACTTCTCCTGATTCCGa                                     <  1:276597/149‑1 (MQ=255)
                                                                                                               aaCTTCATCACGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGa                                     <  4:152585/149‑1 (MQ=255)
                                                                                                                  ttCATCACGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATAc                                  >  2:219746/1‑149 (MQ=255)
                                                                                                                   tcatcaCGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGAtt                                         >  6:266059/1‑141 (MQ=255)
                                                                                                                   tcatcaCGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGAtt                                         <  5:266059/141‑1 (MQ=255)
                                                                                                                    catcaCGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATAc                                  <  3:296308/147‑1 (MQ=255)
                                                                                                                             gCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTAatgat                       <  3:349189/149‑1 (MQ=255)
                                                                                                                               cacaGGCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTAAtgatga                     <  8:135316/149‑1 (MQ=255)
                                                                                                                                    gCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTAATGATGAGCCCa                >  1:67703/1‑149 (MQ=255)
                                                                                                                                       ggTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTAATGATGAGCCCATTc             <  8:372640/149‑1 (MQ=255)
                                                                                                                                            ccccACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTAATGATGAGCCCATTCGATGc        >  3:4039/1‑149 (MQ=255)
                                                                                                                                             cccACAATGCCGCTGCCGACGTTCTCGTCGATGTTGTAAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACATTTTACCTGAAAGCAATTATAAATACTCCTGATTCAGATACCGATTAATAATCACCCCATTCGAtac        >  2:84927/1‑146 (MQ=255)
                                                                                                                                              ccaccaTGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATaaa                                                    <  8:197487/103‑1 (MQ=255)
                                                                                                                                              ccaccaTGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATaaa                                                    >  7:197487/1‑103 (MQ=255)
                                                                                                                                                 ccaTGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACtttt                                                                        >  1:292301/1‑80 (MQ=255)
                                                                                                                                                 ccaTGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACtttt                                                                        <  2:292301/80‑1 (MQ=255)
                                                                                                                                                  caTGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTGTCCTGATTCCGATACCGATTAATGATGAGCCCATTCGATGCTGTTCa  <  4:305648/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TTCCCTTTCACCCATTCGGTGACCAGGGAGCTGGAAGCGATAGCGGAACCGCAGCCGTAAGTTTTAAAACGCGCGTCTTCAATGATACCTTCATCGTTGACTTTAATCTGCAACTTCATCACGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTAATGATGAGCCCATTCGATGCTGTTCA  >  NC_000913/2660051‑2660345

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: