Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,450,065 G→A 6.4% Q4* (CAA→TAA)  PP_1268 ← HtrA suppressor protein SohA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,450,0650GA 6.4% 332.3 / 8.6 109Q4* (CAA→TAA) PP_1268HtrA suppressor protein SohA
Reads supporting (aligned to +/- strand):  ref base G (53/49);  new base A (2/5);  total (55/54)
Fisher's exact test for biased strand distribution p-value = 2.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 125 total aligned reads displayed.

ACGTTTTTTCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAGGTGGCGGTGGCCCATGATGAT  >  NC_002947/1449916‑1450214
                                                                                                                                                     |                                                                                                                                                     
aCGTTTTTTCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGattgcattg                                                                                                                                                       >  2:4228821/1‑150 (MQ=255)
  gTTTTTTCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGcc                                                                                                                                                     <  1:2014927/150‑1 (MQ=255)
  gTTTTTTCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGcc                                                                                                                                                     >  2:4303921/1‑150 (MQ=255)
   ttttttCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCt                                                                                                                                                    >  1:33913/1‑150 (MQ=255)
    tttttCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGCAATGCAATCGCATTACCtt                                                                                                                                                   <  2:3965070/150‑1 (MQ=255)
      tttCAAGCCCAAGGCGCGCCGCTCCGTTTCGGGAACTGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGt                                                                                                                                                 <  2:5110168/150‑1 (MQ=255)
      tttCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGCAATGCAATCGCATTACCTTga                                                                                                                                                 <  1:1586900/150‑2 (MQ=255)
      tttCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGCAATGCAATCGCATTACCTTga                                                                                                                                                 <  1:1169367/150‑2 (MQ=255)
      tttCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGCAATGCAATCGCATTACCTTga                                                                                                                                                 <  1:1169246/150‑2 (MQ=255)
          aaGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAAt                                                                                                                                             <  1:4354722/150‑1 (MQ=255)
            gCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGc                                                                                                                                           >  1:4789444/1‑150 (MQ=255)
             cccAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCa                                                                                                                                          <  1:5040744/150‑1 (MQ=255)
              ccAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAg                                                                                                                                         >  1:123377/1‑150 (MQ=255)
              ccAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAg                                                                                                                                         <  2:74770/150‑1 (MQ=255)
               cAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAgg                                                                                                                                        >  2:2181380/1‑150 (MQ=255)
                aaGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGt                                                                                                                                       >  2:1011847/1‑150 (MQ=255)
                aaGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGt                                                                                                                                       >  1:3774542/1‑150 (MQ=255)
                aaGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGt                                                                                                                                       >  1:2294154/1‑150 (MQ=255)
                aaGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGACAAGGCAATGCaatc                                                                                                                                       <  1:3054738/150‑4 (MQ=255)
                     gcgcCGCACCGTTTCGGGATCAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCAtt                                                                                                                                  <  2:2439206/150‑1 (MQ=255)
                       gccgcACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGAGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTAGTCAATGCAGGTGCATTAc                                                                                                                                >  1:4939247/1‑150 (MQ=255)
                          gcACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCAc                                                                                                                             >  1:4299956/1‑150 (MQ=255)
                           cACCGTTTCGGGAACAGTGTTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACg                                                                                                                            <  2:4715680/150‑1 (MQ=255)
                            acattttCGGGAACAGTTGTTTGATAGCGATCGTTCATCGTGGACTAGACGTCGAAGGTTGAAGCCATGAGTAATGTCCTGTGGGTTACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACgg                                                                                                                           <  2:3487216/146‑1 (MQ=255)
                            aCCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTATTGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGACAAGGCAATGCAATCGCATTAccttga                                                                                                                           <  2:3589538/150‑5 (MQ=38)
                                    gggAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGc                                                                                                                   <  1:820757/150‑1 (MQ=255)
                                    gggAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGc                                                                                                                   >  2:2681178/1‑150 (MQ=255)
                                     ggAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCg                                                                                                                  >  1:2811474/1‑150 (MQ=255)
                                      gAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCgg                                                                                                                 >  1:2385742/1‑150 (MQ=255)
                                        aCAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTc                                                                                                               >  1:1498777/1‑150 (MQ=255)
                                               tttGATAGCGATCGTTCTGCGTGGACTCGACGTCGAAGTTGGCAGCCATGAGTTATTTCATGTTGTTGACAGCGATACAAATCAATGTAATTGGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCaga                                                                                                        <  1:302389/150‑1 (MQ=255)
                                                 tGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGc                                                                                                      >  1:4543901/1‑150 (MQ=255)
                                                   aTAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGccc                                                                                                    >  2:434965/1‑150 (MQ=255)
                                                    tAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCcca                                                                                                   <  2:4234033/150‑1 (MQ=255)
                                                     aGCGATCGGTCAGCGTGGACTCGTCGTCGAAGGTGGCAGGCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCccac                                                                                                  <  2:4209888/150‑1 (MQ=255)
                                                     aGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCccac                                                                                                  <  1:2681178/150‑1 (MQ=255)
                                                            ggTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACt                                                                                           <  1:633982/150‑1 (MQ=255)
                                                               cAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGc                                                                                        <  1:4708737/150‑1 (MQ=255)
                                                                 gCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCAt                                                                                      <  1:2306999/150‑1 (MQ=255)
                                                                 gCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCAt                                                                                      <  2:647460/150‑1 (MQ=255)
                                                                 gCGTGGACTAGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACATCTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCAt                                                                                      <  2:3099618/150‑1 (MQ=255)
                                                                      gACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTa                                                                                 <  2:2296348/150‑1 (MQ=255)
                                                                       aCTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTAc                                                                                >  2:827972/1‑150 (MQ=255)
                                                                        cTCGACGTCGAAGGTGGCATACATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTAcc                                                                               >  2:22922/1‑150 (MQ=255)
                                                                        cTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTAcc                                                                               >  1:4894633/1‑150 (MQ=255)
                                                                        cTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTAcc                                                                               >  2:1549008/1‑150 (MQ=255)
                                                                           gACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCgg                                                                            >  1:3473609/1‑150 (MQ=255)
                                                                           gACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCgg                                                                            >  1:236084/1‑150 (MQ=255)
                                                                           gACGTCGAAGGTGACAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCgg                                                                            <  1:3190683/150‑1 (MQ=255)
                                                                            aCGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGt                                                                           <  1:827972/150‑1 (MQ=255)
                                                                            aCGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGt                                                                           <  1:1022818/150‑1 (MQ=255)
                                                                              gTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCa                                                                         >  1:4317112/1‑150 (MQ=255)
                                                                               tCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCaa                                                                        <  1:1598354/150‑1 (MQ=255)
                                                                               tCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCaa                                                                        <  1:1598509/150‑1 (MQ=255)
                                                                                cGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAAt                                                                       >  2:4078158/1‑150 (MQ=255)
                                                                                 gAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATc                                                                      >  2:3609067/1‑150 (MQ=255)
                                                                                 gAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATc                                                                      >  2:2793988/1‑150 (MQ=255)
                                                                                 gAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATc                                                                      <  2:2181925/150‑1 (MQ=255)
                                                                                      tGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCtt                                                                 <  2:4000800/150‑1 (MQ=255)
                                                                                              aTGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGa                                                         >  2:4725415/1‑150 (MQ=255)
                                                                                              aTGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGa                                                         <  1:4228821/150‑1 (MQ=255)
                                                                                              aTGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGa                                                         >  1:3051357/1‑150 (MQ=255)
                                                                                              aTGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGa                                                         >  1:812947/1‑150 (MQ=255)
                                                                                                     aTGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACgg                                                  <  2:4533643/150‑1 (MQ=255)
                                                                                                        tCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACggggt                                               <  1:3609067/150‑3 (MQ=255)
                                                                                                        tCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTc                                               >  1:1184250/1‑150 (MQ=255)
                                                                                                        tCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTc                                               <  2:4646171/150‑1 (MQ=255)
                                                                                                          cTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACGCTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCa                                             <  2:4920069/150‑1 (MQ=255)
                                                                                                          cTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCa                                             >  2:2980594/1‑150 (MQ=255)
                                                                                                          cTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCa                                             <  1:1458719/150‑1 (MQ=255)
                                                                                                              gggTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCaa                                         >  2:4282296/1‑150 (MQ=255)
                                                                                                                     aGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGcccc                                  >  1:4268488/1‑150 (MQ=255)
                                                                                                                     aGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGcccc                                  >  2:1829381/1‑150 (MQ=255)
                                                                                                                      gCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCa                                 >  2:4892725/1‑150 (MQ=255)
                                                                                                                       cGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCaa                                <  1:4460841/150‑1 (MQ=255)
                                                                                                                           aCAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACcgc                            <  2:3639453/150‑1 (MQ=255)
                                                                                                                           aCAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGAAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACcgc                            >  2:4916975/1‑150 (MQ=255)
                                                                                                                             aaCTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCACCCGCCCCAACCGCGa                          >  2:3148649/1‑150 (MQ=255)
                                                                                                                             aaCTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGa                          >  2:1173541/1‑150 (MQ=255)
                                                                                                                             aaCTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGa                          >  1:982025/1‑150 (MQ=255)
                                                                                                                               cTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACa                        >  2:4235812/1‑150 (MQ=255)
                                                                                                                                  aaGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAggt                     <  2:3656492/150‑1 (MQ=255)
                                                                                                                                  aaGGCAATGCAATCGCATTACCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAggt                     >  1:2924143/1‑150 (MQ=255)
                                                                                                                                  aaGGCAATGCAATCGCATTACCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAggt                     >  2:2611217/1‑150 (MQ=255)
                                                                                                                                     gTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAggtggc                  <  1:946422/150‑1 (MQ=255)
                                                                                                                                         tgcacctgcattgaCAAGGCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAggtggcggtg              >  1:2061342/7‑150 (MQ=255)
                                                                                                                                           cGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAggtggcggtggc            <  2:1483475/150‑1 (MQ=255)
                                                                                                                                            gATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAGGTGGCGGTGGcc           >  2:4333251/1‑150 (MQ=255)
                                                                                                                                            gATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAGGTGGCGGTGGcc           >  2:2955394/1‑150 (MQ=255)
                                                                                                                                            acctgcattgaCAAGGCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAGGTGGCGGTGGcc           >  1:429797/4‑150 (MQ=255)
                                                                                                                                             attgcattgcCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAGGTGGCGGTGGccc          <  2:2584624/150‑1 (MQ=255)
                                                                                                                                              ttgcattgcCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAGGTGGCGGTGGCCCa         >  2:1454353/1‑150 (MQ=255)
                                                                                                                                                gcattgcCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAGGTGGCGGTGGCCCatg       <  2:2811474/150‑1 (MQ=255)
                                                                                                                                                 cattgcCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAGGTGGCGGTGGCCCatga      >  1:4069976/1‑150 (MQ=255)
                                                                                                                                                 cattgcCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAGGTGGCGGTGGCCCatga      >  1:1370379/1‑150 (MQ=255)
                                                                                                                                                  attgcCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAGGTGGCGGTGGCCCatgat     >  1:720054/1‑150 (MQ=255)
                                                                                                                                                   ttgcCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAGGTGGCGGTGGCCCatgatg    >  2:2705830/1‑150 (MQ=255)
                                                                                                                                                    tgcCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTAATTGACGGGTCCAGCAACCGCCCCAACCGCGACAGGTGGCGGTGGCCCAtgatga   <  2:397499/150‑1 (MQ=255)
                                                                                                                                                     gcCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAGGTGGCGGTGGCCCATgatgat  <  2:4939247/150‑1 (MQ=255)
                                                                                                                                                     gcCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAGGTGGCGGTGGCCCATgatgat  <  1:2642507/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACGTTTTTTCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAGGTGGCGGTGGCCCATGATGAT  >  NC_002947/1449916‑1450214

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: