Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,177,715 C→T 5.1% A284T (GCT→ACT)  fbaB ← fructose‑bisphosphate aldolase class I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,177,7150CT 5.1% 248.7 / 2.8 79A284T (GCT→ACT) fbaBfructose‑bisphosphate aldolase class I
Reads supporting (aligned to +/- strand):  ref base C (31/43);  new base T (2/2);  total (33/46)
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.19e-01

CACGCCGTCAGCCATCGATTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTTTGTTGATAACCGCAGTACGCACTGCATCGCTGAGGTCAGTTTCACCGCCCGCAGCACCGCCGGAGTTTATCAACCCAGCCCGACCCATATAGCAGTTAGCTAACTGATAACGCACCAGATCAATCGGGTTTTCGCTGGTCAATTTGCTGTAAACACGATCGTCGGTGTAACCGTAATTAATTGCTTTATAGCCGCCGTTATTTTCCGCCATTTTTTGTTTGACGATA  >  NC_000913/2177569‑2177863
                                                                                                                                                  |                                                                                                                                                    
cACGCCGTCAGCCATCGATTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTTTGTTGATAACCGCAGTACGCACTGCATCGCTGAGGTCAGTTTCACCGCCCGCAGCACCGCCGGAGTTTATCAACCCAGcc                                                                                                                                                     <  5:321918/148‑1 (MQ=255)
 aCGCCGTCAGCCATCGATTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTTTGTTGATAACCGCAGTACGCACTGCATCGCTGAGGTCAGTTTCACCGCCCGCAGCACCGCCGGAGTTTATCAACCCAGccc                                                                                                                                                    >  1:184818/1‑148 (MQ=255)
 aCGCCGTCAGCCATCGATTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTTTGTTGATAACCGCAGTACGCACTGCATCGCTGAGGTCAGTTTCACCGCCCGCAGCACCGCCGGAGTTTATCAACCCAGccc                                                                                                                                                    <  2:184818/148‑1 (MQ=255)
   gCCGTCAGCCATCGATTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTTTGTTGATAACCGCAGTACGCACTGCATCGCTGAGGTCAGTTTCACCGCCCGCAGCACCGCCGGAGTTTATCAACCCAGCCCGAc                                                                                                                                                 <  4:288147/149‑1 (MQ=255)
      gTCAGCCATCGATTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTTTGTTGATAACCGCAGTACGCACTGCATCGCTGAGGTCAGTTTCACCGCCCGCAGCACCGCCGGAGTTTATCAACCCAGCCCGACCCa                                                                                                                                              <  1:152974/149‑1 (MQ=255)
       tCAGCCATCGATTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTTTGTTGATAACCGCAGTACGCACTGCATCGCTGAGGTCAGTTTCACCGCCCGCAGCACCGCCGGAGTTTATCAACCCAGCCCGACCCat                                                                                                                                             <  6:359207/149‑1 (MQ=255)
       tCAGCCATCGATTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTTTGTTGATAACCGCAGTACGCACTGCATCGCTGAGGTCAGTTTCACCGCCCGCAGCACCGCCGGAGTTTATCAACCCAGCCCGACCCat                                                                                                                                             <  1:330281/149‑1 (MQ=255)
         aGCCATCGATTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTTTGTTGATAACCGCAGTACGCACTGCATCGCTGAGGTCAGTTTCACCGCCCGCAGCACCGCCGGAGTTTATCAACCCAGCCCGACCCatat                                                                                                                                           <  4:313403/149‑1 (MQ=255)
          gCCATCGATTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTTTGTTGATAACCGCAGTACGCACTGCATCGCTGAGGTCAGTTTCACCGCCCGCAGCACCGCCGGAGTTTATCAACCCAGCCCGACCCAtata                                                                                                                                          <  7:395366/149‑1 (MQ=255)
            cATCGATTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTTTGTTGATAACCGCAGTACGCACTGCATCGCTGAGGTCAGTTTCACCGCCCGCAGCACCGCCGGAGTTTATCAACCCAGCCCGACCCATATAg                                                                                                                                         >  7:292065/1‑148 (MQ=255)
             aTCGATTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTTTGTTGATAACCGCAGTACGCACTGCATCGCTGAGGTCAGTTTCACCGCCCGCAGCACCGCCGGAGTTTATCAACCCAGCCCGACCCATATAg                                                                                                                                         >  1:97628/1‑147 (MQ=255)
              tCGATTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTTTGTTGATAACCGCAGTACGCACTGCATCGCTGAGGTCAGTTTCACCGCCCGCAGCACCGCCGGAGTTTATCAACCCAGCCCGACCCata                                                                                                                                            >  5:233824/1‑143 (MQ=255)
              tCGATTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTTTGTTGATAACCGCAGTACGCACTGCATCGCTGAGGTCAGTTTCACCGCCCGCAGCACCGCCGGAGTTTATCAACCCAGCCCGACCCata                                                                                                                                            <  6:233824/143‑1 (MQ=255)
                gATTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTTTGTTGATAACCGCAGTACGCACTGCATCGCTGAGGTCAGTTTCACCGCCCGCAGCACCGCCGGAGTTTATCAACCCAGCCCGACCCATATAGCAGt                                                                                                                                     <  6:281978/148‑1 (MQ=255)
                 aTTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTTTGTTGATAACCGCAGTACGCACTGCATCGCTGAGGTCAGTTTCACCGCCCGCAGCACCGCCGGAGTTTATCAACCCAGCCCGACCCATATAGCAGt                                                                                                                                     >  1:222481/1‑147 (MQ=255)
                 aTTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTTTGTTGATAACCGCAGTACGCACTGCATCGCTGAGGTCAGTTTCACCGCCCGCAGCACCGCCGGAGTTTATCAACCCAGCCCGACCCATATAGCAGTTa                                                                                                                                   <  6:227956/149‑1 (MQ=255)
                    tCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTTTGTTGATACCCGCAGTACGCACTGCATCGCTGAGGTCAGTTTCACCGCCCGCAGCACCGCCGGAGTTTATCAACCCAGCCCGACCCATATAGCAGTTAGc                                                                                                                                 >  3:237774/1‑148 (MQ=255)
                      ttGATCGCTTTACGTCCAATTATCAGCCCCATTCCGCCTGCTAGTTTGTTGATAACCTCAGTAAGCACTGCATCGCTGACGTCAGTTTCAGCGCCCGCCTCAACGCCTGAGTTTATCAACCAATCCCGAACCATATAGCAGTTAGCTaa                                                                                                                              <  2:25479/149‑1 (MQ=255)
                       tGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTTTGTTGATAACCGCAGTACGCACTGCATCGCTGAGGTCAGTTTCACCGCCCGCAGCACCGCCGGAGTTTATCAACCCAGTCCGACCCATATAGCAGTTAGCTAAc                                                                                                                             >  7:162076/1‑149 (MQ=255)
                       tGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTTTGTTGATAACCGCAGTACGCACTGCATCGCTGAGGTCAGTTTCACCGCCCGCAGCACCGCCGGAGTTTATCAACCCAGCCCGACCCATATAGCAGTTAGCTAAc                                                                                                                             >  8:1538/1‑149 (MQ=255)
                             cTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTTTGTTGATAACCGCAGTACGCACTGCATCGCTGAGGTCAGTTTCACCGCCCGCAGCACCGCCGGAGTTTATCAACCCAGCCCGACCCATATAGCAGTTAGCTAACTGATaa                                                                                                                       <  8:355982/149‑1 (MQ=255)
                             cTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTTTGTTGATAACCGCAGTACGCACTGCATCGCGGAGGGCAGTTTCACCGCCCGCAGCACCGCCGGAGTTTATCAACCCAGCCCGACCCATATAGCAGTTAGCTAACTGATaa                                                                                                                       <  8:328569/149‑1 (MQ=255)
                              tttACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTTTGTTGATAACCGCAGTACGCACTGCATCGCTGAGGTCAGTTTCACCGCCCGCAGCACCGCCGGAGTTTATCAACCCAGCCCGACCCATATAGCAGTTAGCTAACTGATAAc                                                                                                                      >  2:358147/1‑149 (MQ=255)
                              tttACGTCCAAGAATCAGCCCCATTCAGCCTGCGCGTTTGTTGATAACCGCAGTACGCACTGCATCGCTGAGGTCAGTTTCACCGCCCGCAGCACCGCCGGAGTTTATCAACCCAGCCCGACCCATATAGCAGTTAGCTAACTGATAAc                                                                                                                      >  8:73476/1‑149 (MQ=255)
                                 aCGTCCAAGAATCAGCCCCATTCCGCCTGCTCGTTTGTTGATAACCGCAGTACGCACTGCATCGCTGAGGTCAGTTTCACCTCACGCAGCACCGCCGGAGTTTATCAACCCAGCCCGACCCATATAGCAGTTAGCTAACTGATATCGCa                                                                                                                   >  6:18135/1‑149 (MQ=255)
                                 aCGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTTTGTTGATAACCGCAGTACGCACTGCATCGCTGAGGTCAGTTTCACCGCCCGCAGCACCGCCGGAGTTTATCAACCCAGCCCGACCCATATAGCAGTTAGCTAACTGATAACGc                                                                                                                    <  5:341972/148‑1 (MQ=255)
                                  cGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTTTGTTGATAACCGCAGTACGCACTGCATCGCTGAGGTCAGTTTCACCGCCCGCAGCACCGCCGGAGTTTATCAACCCAGCCCGACCCATATAGCAGTTAGCTAACTGATAACGCAc                                                                                                                  <  5:76905/149‑1 (MQ=255)
                                   gTCCAAGAATCAGCCCCATTCCGCCTGCGCGTTTGTTGATAACCGCAGTACGCACTGCATCGCTGAGGTCAGTTTCACCGCCCGCAGCACCGCCGGAGTTTATCAACCCAGCCCGACCCATATAGCAGTTAGCTAACTGATAACGCAcc                                                                                                                 <  2:157132/149‑1 (MQ=255)
                                    tCCAAGAATCAGCCCCATTCCGCCTGCGCGTTTGTTGATAACCGCAGTACGCACTGCATCGCTGAGGTCAGTTTCACCGCCCGCAGCACCGCCGGAGTTTATCAACCCAGCCCGACCCATATAGCAGTTAGCTAACTGATAACGCACCa                                                                                                                <  7:280805/149‑1 (MQ=255)
                                       aaGAATCAGCCCCATTCCGCCTGCGCGTTTGTTGATAACCGCAGTACGCACTGCATCGCTGAGGTCAGTTTCACCGCCCGCAGCACCGCCGGAGTTTATCAACCCAGCCCGACCCATATAGCAGTTAGCTAACTGATAACGCAc                                                                                                                  >  7:356672/1‑144 (MQ=255)
                                        aGAATCAGCCACATTCCGCCTGCGCGTTTGTTGATAACCGCAGTACGCACTGCATCGCTGAGGTCAGTTTCACCGCCCGCAGCACCGCCGGAGTTTATCAACCCAGCCCGACCCATATAGCAGTTAGCTAACTGATAACGCACCAGAt                                                                                                             >  2:121447/1‑148 (MQ=255)
                                         gAATCAGCCCCATTCCGCCTGCGCGTTTGTTGATAACCGCAGTACGCACTGCATCGCTGAGGTCAGTTTCACCGCCCGCAGCACCGCCGGAGTTTATCAACCCAGCCCGACCCATATAGCAGTTAGCTAACTGATAACGCACCAGATc                                                                                                            >  7:125427/1‑148 (MQ=255)
                                              aGCCCCATTCCGCCTGCGCGTTTGTTGATAACCGCAGTACGCACTGCATCGCTGAGGTCAGTTTCACCGCCCGCAGCACCGCCGGAGTTTATCAACCCAGCCCGACCCATATAGCAGTTAGCTAACTGATAACGCACCAGATCAATCgg                                                                                                      <  1:200765/149‑1 (MQ=255)
                                              aGCCCCATTCCGCCTGCGCGTTTGTTGATAACCGCAGTACGCACTGCATCGCTGAGGTCAGTTTCACCGCCCGCAGCACCGCCGGAGTTTATCAACCCAGCCCGACCCATATAGCAGTTAGCTAACTGATAACGCACCAGATCAATCgg                                                                                                      <  4:237774/149‑1 (MQ=255)
                                                 cccATTCCGCCTGCGCGTTTGTTGATAACCGCAGTACGCACTGCATCGCTGAGGTCAGTTTCACCGCCCGCAGCACCGCCGGAGTTTATCAACCCAGCCCGACCCATATAGCAGTTAGCTAACTGATAACGCACCAGATCAATCGGGtt                                                                                                   <  8:292065/149‑1 (MQ=255)
                                                   cATTCCGCCTGCGCGTTTGTTGATAACCGCAGTACGCACTGCATCGCTGAGGTCAGTTTCACCGCCCGCAGCACCGCCGGAGTTTATCAACCCAGCCCGACCCATATAGCAGTTAGCTAACTGATAACGAACCAGATCAATCGGGtttt                                                                                                 >  6:181832/1‑149 (MQ=255)
                                                        cGCCTGCGCGTTTGTTTATAACCGCAGTACGCACTGCATCGCTGAGGTCAGTTTCACCGCCCGCAGCACCGCCGGAGTTTATCAACCCAGCCCGACCCATATAGCAGTTAGCTAACTGATAACGCACCAGATCAATCGGGTTTTCGCTg                                                                                            <  6:140045/149‑1 (MQ=255)
                                                        cGCCTGCGCGTTTGTTGATAACCGCAGTACGCACTGCATCGCTGAGGTCAGTTTCACCGCCCGCAGCACCGCCGGAGTTTATCAACCCAGCCCGACCCATATAGCAGTTAGCTAACTGATAACGCACCAGATCAATCGGGTTTTCGCt                                                                                             <  6:305953/148‑1 (MQ=255)
                                                           cTGCGCGTTTGTTGATAACCGCAGTACGCACTGCATCGCTGAGGTCAGTTTCACCGCCCGCAGCACCGCCGGAGTTTATCAACCCAGCCCGACCCATATAGCAGTTAGCTAACTGATAACGCACCAGATCAATCGGGTTTTCGCTGGt                                                                                          >  7:38688/1‑148 (MQ=255)
                                                            tGCGCGTTTGTTGATAACCGCAGTACGCACTGCATCGCTGAGGTCAGTTTCACCGCCCGCAGCACCGCCGGAGTTTATCAACCCAGCCCGACCCATATAGCAGTTAGCTAACTGATAACGCACCAGATCAATCGGGTTTTCGCTGGTc                                                                                         >  6:236647/1‑148 (MQ=255)
                                                                 gTTTGTTGATAACCGCAGTACGCACTGCATCGCTGAGGTCAGTTTCACCGCCCGCAGCACCGCCGGAGTTTATCAACCCAGCCCGACCCATATAGCAGTTAGCTAACTGATAACGCACCAGATCAATCGGGTTTTCGCTGGTCAATTTg                                                                                   <  1:169401/149‑1 (MQ=255)
                                                                 gTTTGTTGATAACCGCAGTACGCACTGCATCGCTGAGGTCAGTTTCACCGCCCGCAGCACCGCCGGAGTTTATCAACCCAGCCCGACCCATATAGCAGTTAGCTAACTGATAACGCACCAGATCAATCGGGTTTTCGCTGGTCAATTTg                                                                                   >  6:53452/1‑149 (MQ=255)
                                                                  tttGTTGATAACCGCAGTACGCACTGCATCGCTGAGGTCAGTTTCACCGCCCGCAGCACCGCCGGAGTTTATCAACCCAGCCCGACCCATATAGCAGTTAGCTAACTGATAACGCACCAGATCAATCGGGTTTTCGCTGGTCAATTTGc                                                                                  <  5:53452/149‑1 (MQ=255)
                                                                        gATAACCGCAGTACGCACTGCATCGCTGAGGTCAGTTTCACCGCCCGCAGCACCGCCGGAGATTATCAACCCAGCCCGACCCATATAGCAGTTAGCTAACTGATAACGCACCAGATCAATCGGGTTTTCGCTGGTCAATTTGCTGTaa                                                                             <  1:315752/148‑1 (MQ=255)
                                                                         aTAACCGCAGTACGCACTGCATCGCTGAGGTCAGTTTCACCGCCCGCAGCACCGCCGGAGTTTATCAACCCAGCCCGACCCATATAGCAGTTAGCTAACTGATAACGCACAAGATCAATCGGGTTTTCGCTGGTCAATTTGCTGTaaa                                                                            >  8:165508/1‑148 (MQ=255)
                                                                          tAACCGCAGTACGCACTGCATCGCTGAGGTCAGTTTCACCGCCCGCAGCACCGCCGGAGTTTATCAACCCAGCCCGACCCATATAGCAGTTAGCTAACTGATAACGCACCAGATCAATCGGGTTTTCGCTGGTCAATTTGCTGTAAaca                                                                          <  8:125427/149‑1 (MQ=255)
                                                                                  gTACGCACTGCATCGCTGAGGTCAGTTTCACCGCCCGCAGCACCGCCGGAGTTTATCAACCCAGCCCGACCCATATAGCAGTTAGCTAACTGATAACGCACCAGATCAATCGGGTTTTCGCTGGTCAATTTGCTGTAAACACGAtcgtc                                                                  <  5:236647/149‑1 (MQ=255)
                                                                                  gTACGCACTGCATCGCTGAGGTCAGTTTCACCGCCCGCAGCACCGCCGGAGTTTATCAACCCAGCCCGACCCATATAGCAGTTAGCTAACTGATAACGCACCAGATCAATCGGGTTTTCGCTGGTCAATTTGCTGTAAACACGAtcgtc                                                                  <  7:1538/149‑1 (MQ=255)
                                                                                   tACGCACTGCATCGCTGAGGTCAGTTTCACCGCCCGCAGCACCGCCGGAGTTTATCAACCCAGCCCGACCCATATAGCAGTTAGCTAACTGATAACGCACCAGATCAATCGGGTTTTCGCTGGTCAATTTGCTGTAAACACGAt                                                                      >  5:365313/1‑144 (MQ=255)
                                                                                   tACGCACTGCATCGCTGAGGTCAGTTTCACCGCCCGCAGCACCGCCGGAGTTTATCAACCCAGCCCGACCCATATAGCAGTTAGCTAACTGATAACGCACCAGATCAATCGGGTTTTCGCTGGTCAATTTGCTGTAAACACGAt                                                                      <  6:365313/144‑1 (MQ=255)
                                                                                               cGCTGAGGTCAGTTTCACCGCCCGCAGCACCGCCGTAGTTTATCAACCCCGCCCGACCCATATAGCAGTTAGCTACCTGATACCGCACCAGATCAATCGGGTTTTCGCTGGTCAATTTGCTGTAAACACGATCGTCGGTGTAACCGtaa                                                     >  3:31378/1‑149 (MQ=255)
                                                                                                    aGGTCAGTTTCACCGCCCGCAGCACCGCCGGAGTTTATCAACCCAGCCCGACCCATATAGCAGTTAGCTAACTGATAACGCACCAGATCAATCGGGTTTTCGCTGGTCAATTTGCTGTAAACACGATCGTCGGTGTAACCGtaattaa                                                 <  2:222481/148‑1 (MQ=255)
                                                                                                      gTCAGTTTCACCGCCCGCAGCACCGCCGGAGTTTATCAACCCAGCCCGACCCATATAGCAGTTAGCTAACTGATAACGCACCAGATCAATCGGGTTTTCGCTGGTCAATTTGCTGTAAACACGATCGTCGGTGTAACCGTAATTAATTg                                              >  1:195953/1‑149 (MQ=255)
                                                                                                       tCAGTTTCACCGCCCGCAGCACCGCCGGAGTTTATCAACCCAGTCCGACCCATATAGCAGTTAGCTAACTGATAACGCACCAGATCAATCGGGTTTTCGCTGGTCAATTTg                                                                                   <  1:169198/111‑1 (MQ=255)
                                                                                                       tCAGTTTCACCGCCCGCAGCACCGCCGGAGTTTATCAACCCAGTCCGACCCATATAGCAGTTAGCTAACTGATAACGCACCAGATCAATCGGGTTTTCGCTGGTCAATTTg                                                                                   >  2:169198/1‑111 (MQ=255)
                                                                                                         aGTTTCACCGCCCGCAGCACCGCCGGAGTTTATCAACCCAGCCCGACCCATATAGCAGTTAGCTAACTGATAACGCACCAGATCAATCGGGTTTTCGCTGGTCAATTTGCTGTAAACACGATCGTCGGTGTAACCGTAATTAATTGCt                                            <  2:195953/148‑1 (MQ=255)
                                                                                                         aGTTTCACCGCCCGCAGCACCGCCGGAGTTTATCAACACAGCCCGTACCATATAGCAGTTAGCTAACTGATAACGCACCAGATCAATCGGGTTTTCGCTGGTCAATTTGCTGTAAACACGATCGTCGGTGTAACCGTAATTAATTGCtt                                           <  1:358147/149‑1 (MQ=255)
                                                                                                          gTTTCACCGCCCGCAGCACCGCCGGAGTTTATCAACCCAGCCCGACCCATATAGCAGTTAGCTAACTGATAACGCACCAGATCAATCGGGTTTTCGCTGGTCAATTTGCTGTAAACACGATCGTCGGTGTAACCGTAATTAATTGCttt                                          <  7:165508/149‑1 (MQ=255)
                                                                                                          gTTTCACCGCCCGCAGCACCGCCGGAGTTTATCAACCCAGCCCGACCCATATAGCAGTTAGCTAACTGATAACGCACCAGATCAATCGGGTTTTCGCTGGTCAATTTGCTGTAAACACGATCGTCGGTGTAACCGTAATTAATTGCttt                                          <  7:73476/149‑1 (MQ=255)
                                                                                                           tttCACCGCCCGCAGCACCGCCGGAGTTTATCAACCCAGCCCGACCCATATAGCAGTTAGCTAACTGATAACGCACCAGATCAATCGGGTTTTCGCTGGTCAATTTGCTGTAAACACGATCGTCGGTGTAACCGTAATTAATTGCTTta                                         <  8:356672/149‑1 (MQ=255)
                                                                                                           tttCACCGCCCGCAGCACCGCCGGAGTTTATCAACCCAGCCCGACCCATATAGCAGTTAGCTAACTGATAACGCACCAGATCAATCGGGTTTTCGCTGGTCAATTTGCTGTAAACACGATCGTCGGTGTAACCGTAATTAATTGCTTta                                         <  8:38688/149‑1 (MQ=255)
                                                                                                           tttCACCGCCCGCAGCACCGCCGGAGTTTATCAACCCAGCCCGACCCATATAGCAGTTAGCTAACTGATAACGCACCAGATCAATCGGGTTTTCGCTGGTCAATTTGCTGTAAACACGATCGTCGGTGTAACCGTAATTAATTGCTTta                                         <  2:97628/149‑1 (MQ=255)
                                                                                                                  gcccgcAGCGCCGCCGGAGGTTATCACCCCAGCCCGAGCCATATAGCAGTTAGCTAACTGATACCGCCCCAGATCTATCGGGTGTTCGCTGGTCAATTTGGTGTAAACACGATCGTCGGTGTTACCGTAATTAATTGGTTTATAgccgc                                  >  2:260883/1‑149 (MQ=255)
                                                                                                                      gcagcaCCGCCGGAGTTTATCAACCCAGCCCGCCCCATATAGCAGTTAGCTAACTGATAACGCACCAGATCAATCGGGTTTTCGCTGGTCAATATGCTGTAAACACGCTCGTCGGTGTAACCGTAATAAAGTGCTTTATAGccgccgg                               >  4:319874/1‑147 (MQ=255)
                                                                                                                      gcagcaACGCCGGAGTTTATCAACCCAGCCCGACCCATATAGCAGTTAGCTAACTGATAACGCACCAGATCAATCGGGTTTTCGCTGGTCAATTTGCTGTAAACACGATCGTCGGTGTAACCGTAATTAATTGCTTTATAGCCGCCGt                               >  8:145228/1‑148 (MQ=255)
                                                                                                                          caCCGCCGGAGTTTATCAACCCAGCCCGACCCATATAGCAGTTAGCTAACTGATAACGCACCAGATCAATCGGGTTTTCGCTGGTCAATTTGCTGTAAACACGATCGTCGGTGTAACCGTAATTAATTGCTTTATAGCCGCCGTTAttt                          <  7:87134/149‑1 (MQ=255)
                                                                                                                            ccgactgAGTTTATCAACCCAGACCGACCCATATAGAAGTTAGCTAACTGATAACTCACCAGATGAATCGTGTTTTCGCTGGTCAATTTGCTGTATACACGATCGTCGGTGTAACCGTAATTAATTGCTTTATAGCCGCCGTTATTTTc                        <  5:18135/143‑1 (MQ=255)
                                                                                                                             cgccgGTGTTTATCAACCCAGCCCGACCCATATAGCAGTTAGCTAACTGATAACGCACCAGATCAATCGGGTTTTCGCTGGTCAATTTGCTGTAAACACGATCGTCGGTGTAACCGTAATTAATTGCTTTATAGCCGCCGTTATTTTcc                       <  7:145228/149‑1 (MQ=255)
                                                                                                                              gccgGAGTTTATCAACCCAGCCCGACCCATATAGCAGTTAGCTAACTGATAACGCACCAGATCAATCGGGTTTTCGCTGGTCAATTTGCTGTAAACACGATCGTCGGTGTAACCGTAATTAATTGCTTTATAGCCGCCGTTATTTTCCg                      >  2:299227/1‑149 (MQ=255)
                                                                                                                                  gAGTTTATCAACCCAGCCCGACCCATATAGCAGTTAGCTAACTGATAACGCACCAGATCAATCGGGTTTTCGCTGGTCAATTTGCTGTAAACACGATCGTCGGTGTAACCGTAATTAATTGCTTTATAGCCGCCGTTATTTTCCGCCAt                  >  2:126392/1‑149 (MQ=255)
                                                                                                                                     tttATCAACCCAGCCCGACCCATATAGCAGTTAGCTAACTGATAACGCACCAGATCAATCGGGTTTTCGCTGGTCAATTTGCTGTAAACACGATCGTCGGTGTAACCGTAATTAATTGCTTTATAGCCGCCGTTATTTTCCGCCAtttt               >  8:44781/1‑149 (MQ=255)
                                                                                                                                        aTCAACCCAGCCCGACCCATATAGCAGTTAGCTAACTGATAACGCACCAGATCAATCGGGTTTTCGCTGGTCAATTTGCTGTAAACACGATCGTCGGTGTAACCGTAATTAATTGCTTTATAGCCGCCGTTATTTTCCGCCAtttttt             <  1:121447/148‑1 (MQ=255)
                                                                                                                                        aTCAACCCAGCCCGACCCATATAGCAGTTAGCTAACTGATAACGCACCAGATCAATCGGGTTTTCGCTGGTCAATTTGCTGTAAACACGATCGTCGGTGTAACCGTAATTAATTGCTTTATAGCCGCCGTTATTTTCCGCCAtttttt             <  7:321878/148‑1 (MQ=255)
                                                                                                                                        aTCAACCCAGCCCGACCCATATAGCAGTTAGCTAACTGATAACGCACCAGATCAATCGGGTTTTCGCTGGTCAATTTGCTGTAAACACGATCGTCGGTGTAACCGTAATTAATTGCTTTATAGCCGCCGTTATTTTCCGCCATTTtttg            >  3:22891/1‑149 (MQ=255)
                                                                                                                                           aaCCCAGTCCGACCCATATAGCAGTTAGCTAACTGATAACGCACCAGATCAATCGGGTTTTCGCTGGTCAATTTGCTGTAAACACGATCGTCGGTGTAACCGTAATTAATTGCTTTATAGCCGCCGTTATTTTCCGCCATTTtttg            <  8:162076/146‑1 (MQ=255)
                                                                                                                                           aaCCCAGCCCGACCCATATAGCAGTTAGTTAACTGATAACGCACCAGATCAATCGGGTTTTCGCTGGTCAATTTGCTGTAAACACGATCGTCGGTGTAACCGTAATTAATTGCTTTATAGCCGCCGTTATTTTCCGCCATTTtttgttt         >  2:77788/1‑149 (MQ=255)
                                                                                                                                             cccAGCCCGACCCATATAGCAGTTAGCTAACTGATAACGCACCAGATCAATCGGGTTTTCGCTGGTCAATTTGCTGTAAACACGATCGTCGGTGTAACCGtaa                                                     >  1:193180/1‑103 (MQ=255)
                                                                                                                                             cccAGCCCGACCCATATAGCAGTTAGCTAACTGATAACGCACCAGATCAATCGGGTTTTCGCTGGTCAATTTGCTGTAAACACGATCGTCGGTGTAACCGtaa                                                     <  2:193180/103‑1 (MQ=255)
                                                                                                                                              ccAGCCCGACCCATATAGCAGTTAGCTAACTGATAACGCACCAGATCAATCGGGTTTTCGCTGGTCAATTTGCTGTAAACACGATCGTCGGTGTAACCGTAATTAATTGCTTTATAGCCGCCGTTATTTTCCGCCATTTTTTGTTTGAc      <  1:260883/149‑1 (MQ=255)
                                                                                                                                              ccAGCCCGACCCATATAGCAGTTAGCTAACTGATAACGCACCAGATCAATCGGGTTTTCGCTGGTCAATTTGCTGTAAACACGATCGTCGGTGTAACCGTAATTAATTGCTTTATAGCCGCCGTTATTTTCCGCCATTTTTTGTTTGAc      <  3:158406/149‑1 (MQ=255)
                                                                                                                                                aGCCCGACCCATATAGCAGTTAGCTAACTGATAACGCACCAGATCAATCGGGTTTTCGCTGGTCAATTTGCTGTAAACACGATCGTCGGTGTAACCGTAATTAATTGCTTTATAGCCGCCGTTATTTTCCGCCATTTTTTGTTTGACGa    >  3:62076/1‑149 (MQ=255)
                                                                                                                                                aGCCCGACCCATATAGCAGTTAGCTAACTGATAACGCACCAGATCAATCGGGTTTTCGCTGGTCAATTTGCTGTAAACACGATCGTCGGTGTAACCGTAATCAATTGCTTTATAGCCGCCGTTATTTTCCGCCATTTTTTGTTTGACGa    >  2:303919/1‑149 (MQ=255)
                                                                                                                                                  cccGACCCATATAGCAGTTAGCTAACTGATAACGCACCAGATCAATCGGGTTTTCGCTGGTCAATTTGCTGTAAACACGATCGTCGGTGTAACCGTAATTAATTGCTTTATAGCCGCCGTTATTTTCCGCCATTTTTTGTTTGACGata  >  4:24022/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CACGCCGTCAGCCATCGATTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTGCGCGTTTGTTGATAACCGCAGTACGCACTGCATCGCTGAGGTCAGTTTCACCGCCCGCAGCACCGCCGGAGTTTATCAACCCAGCCCGACCCATATAGCAGTTAGCTAACTGATAACGCACCAGATCAATCGGGTTTTCGCTGGTCAATTTGCTGTAAACACGATCGTCGGTGTAACCGTAATTAATTGCTTTATAGCCGCCGTTATTTTCCGCCATTTTTTGTTTGACGATA  >  NC_000913/2177569‑2177863

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