Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 689,079 C→A 10.1% E445* (GAA→TAA)  miaB ← tRNA‑i(6)A37 methylthiotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273689,0790CA10.1% 1045.1 / 63.7 318E445* (GAA→TAA) miaBtRNA‑i(6)A37 methylthiotransferase
Reads supporting (aligned to +/- strand):  ref base C (156/127);  new base A (20/15);  total (176/142)
Fisher's exact test for biased strand distribution p-value = 8.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 321 total aligned reads displayed.

AAAGATGGGGATAAAGAGAGAAAAAACAAGGCCCACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTCGAAGTTGACCACGCGGTTATTTTCCGTACGCCCGGAAAGCTCCATGATGCTCT  >  CP009273/688931‑689228
                                                                                                                                                    |                                                                                                                                                     
aaaGATGGGGATAAAGAGAGAAAAAACAAGGCCCACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCa                                                                                                                                                      >  2:4536150/1‑150 (MQ=255)
  aGATGGGGATAAAGAGAGAAAAAACAAGGCCCACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATc                                                                                                                                                    <  1:2164689/150‑1 (MQ=255)
   gATGGGGATAAAGAGAGAAAAAACAAGGCCCACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCt                                                                                                                                                   >  2:2696462/1‑150 (MQ=255)
    aTGGGGATAAAGAGAGAAAAAACAAGGCCCACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTTCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCtt                                                                                                                                                  >  1:190942/1‑150 (MQ=255)
     tGGGGATAAAGAGAGAAAAAACAAGGCCCACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTc                                                                                                                                                 >  1:2589070/1‑150 (MQ=255)
      ggggATAAAGAGAGAAAAAACAAGGCCCACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCa                                                                                                                                                >  2:3572045/1‑150 (MQ=255)
       gggATAAAGAGAGAAAAAACAAGGCCCACCGGAACGGCAGGCCTGAGAATTACGGCTTATAATAACACACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAg                                                                                                                                               <  1:1792003/150‑1 (MQ=255)
         gATAAAGAGAGAAAAAACAAGGCCCACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTa                                                                                                                                             <  1:4299224/150‑1 (MQ=255)
          aTAAAGAGAGAAAAAACAAGGCCCACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTAc                                                                                                                                            >  2:954842/1‑150 (MQ=255)
          aTAAAGAGAGAAAAAACAAGGCCCACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTAc                                                                                                                                            >  2:855246/1‑150 (MQ=255)
           tAAAGAGAGAAAAAACAAGGCCCACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACg                                                                                                                                           >  1:892782/1‑150 (MQ=255)
           tAAAGAGAGAAAAAACAAGGCCCACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTATGCCACGCGCAGACCCAATTCATCTTCAGTACg                                                                                                                                           >  1:4054077/1‑150 (MQ=255)
              agagagAAAAAACAAGGCCCACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGcac                                                                                                                                        >  2:1466525/1‑150 (MQ=255)
              agagagAAAAAACAAGGCCCACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGcac                                                                                                                                        <  1:4182714/150‑1 (MQ=255)
              agagagAAAAAACAAGGCCCACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTGTTCTTTGCGGGTACGGGCAATCACTGATTCCGGGGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGcac                                                                                                                                        >  1:4139614/1‑150 (MQ=255)
                  agaAAAAACAAGGCCCACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTTAGGGTACGGGCAATAACTGATTACGGTTTTTATGCCACGCGCAGACCCATTTCATCATAAGTAAGcaccaca                                                                                                                                    >  1:4012804/1‑149 (MQ=255)
                    aaaaaaCAAGGCCCACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACttt                                                                                                                                  >  2:21663/1‑150 (MQ=255)
                    aaaaaaCAAGGCCCACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACttt                                                                                                                                  >  2:1110971/1‑150 (MQ=255)
                      aaaaCAAGGCCCACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCAATTCATCTTCAGTAAGCACCACTTTaa                                                                                                                                >  2:2803161/1‑149 (MQ=255)
                       aaaCAAGGCCCACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTAcc                                                                                                                               >  1:4720936/1‑150 (MQ=255)
                        aaCAAGGCCCACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACcg                                                                                                                              >  1:3558960/1‑150 (MQ=255)
                               cccACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCagag                                                                                                                       <  2:2963921/150‑1 (MQ=255)
                                  aCCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTAATCTTCAGTACGCACCACTTTACCGCGCAGAGAGt                                                                                                                    <  2:4091021/150‑1 (MQ=255)
                                         cGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAATGTCGATTTATTTTCTGGTACGGGCAAACACTGATTCACGTGTTTATGCCACGCGCAGACCCATTTAATCTTCAGAAAGCACCACTAAACACCGCAGACAATACGGGTa                                                                                                             >  1:671541/1‑150 (MQ=255)
                                         cGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTTCGGGTACGGGCCAACACTGATTCCGGTGTTTGTGCCACGCGCCGACCCATTTCATCTTCAGGAAGCACCCCATTAACGCGCAGAGAGTTAGGGTa                                                                                                             >  1:833044/1‑150 (MQ=255)
                                         cGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTa                                                                                                             <  1:2492902/150‑1 (MQ=255)
                                         cGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACAGATTCCGGTGTTTCTGCCACGCGCAGACCCAATTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTa                                                                                                             >  1:3060744/1‑150 (MQ=255)
                                          ggCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCCGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAg                                                                                                            >  1:4584832/1‑150 (MQ=255)
                                          ggCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGATTTCGGGTAg                                                                                                            >  1:3975902/1‑150 (MQ=255)
                                          ggCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTGTTCTTTGCGGGTACGGGCCATCACTGATTCCGGTGTTTCTGCCACGCGCCGACCCCATTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAg                                                                                                            >  1:3696513/1‑150 (MQ=255)
                                           gCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGa                                                                                                           <  1:3277686/150‑1 (MQ=255)
                                           gCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTATTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTAACGCGCAGAGAGTTCGGGTAGa                                                                                                           >  1:2240792/1‑150 (MQ=255)
                                             aGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTAATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACg                                                                                                         >  2:1195357/1‑150 (MQ=255)
                                             aGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTATTTGCGGGGACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCCGACCCAATTCATCTTCAGTACGCACCACATTACCCCGCAGAGAGTTCGGGTAGACg                                                                                                         >  1:2920813/1‑150 (MQ=255)
                                                 cTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCgg                                                                                                     >  2:1643490/1‑150 (MQ=255)
                                                        ttACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTc                                                                                              <  2:2279592/150‑1 (MQ=255)
                                                           cGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTTCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCAATTCATCTTCAGTACGCACCACTTTAACGCGCAGAGAGTTAGGGTAGACGTCGGTAATTTCTAc                                                                                           >  1:227628/1‑150 (MQ=255)
                                                             gCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACaa                                                                                         >  2:517819/1‑149 (MQ=255)
                                                             gCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACAt                                                                                         <  1:2199903/150‑1 (MQ=255)
                                                              cTGATAATAACCCACGCCAAGGTCGTGTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATc                                                                                        >  2:3038963/1‑150 (MQ=255)
                                                                     taaCCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTACGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAAt                                                                                 >  2:1025326/1‑150 (MQ=255)
                                                                      aaCCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAAtt                                                                                <  2:3550796/150‑1 (MQ=255)
                                                                       aCCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAAttt                                                                               <  2:4351467/150‑1 (MQ=255)
                                                                       aCCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTAATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAAttt                                                                               <  1:1665856/150‑1 (MQ=255)
                                                                       aCCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGACAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAAttt                                                                               <  2:696687/150‑1 (MQ=255)
                                                                        cccACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTa                                                                              <  1:513525/150‑1 (MQ=255)
                                                                         ccACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTAc                                                                             <  2:1818994/150‑1 (MQ=255)
                                                                            cGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGa                                                                          <  1:3428237/150‑1 (MQ=255)
                                                                             gCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGAt                                                                         <  1:468850/150‑1 (MQ=255)
                                                                               cAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCa                                                                       <  1:3883622/150‑1 (MQ=255)
                                                                                 aGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCata                                                                     <  1:1405256/150‑1 (MQ=255)
                                                                                 aGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCata                                                                     >  2:4279720/1‑150 (MQ=255)
                                                                                  ggTCGTTTTCTTTGCGGGAACGGGCAATCCCTGATTCCGGGGGGTCTGCCACGCGCAGACCCAATTCAACGTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCatat                                                                    <  2:3465269/150‑1 (MQ=255)
                                                                                    tCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATcc                                                                  <  1:855246/150‑1 (MQ=255)
                                                                                    tCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATcc                                                                  <  2:4188077/150‑1 (MQ=255)
                                                                                    tCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTAATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATcc                                                                  <  2:1352906/150‑1 (MQ=255)
                                                                                      gTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCgg                                                                >  2:4066254/1‑150 (MQ=255)
                                                                                       ttttCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTAATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGc                                                               <  2:639622/150‑1 (MQ=255)
                                                                                        tttCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCg                                                              >  2:1039680/1‑150 (MQ=255)
                                                                                        tttCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCg                                                              >  1:946983/1‑150 (MQ=255)
                                                                                         ttCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGt                                                             >  1:2003222/1‑150 (MQ=255)
                                                                                             ttGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCt                                                         <  1:3479476/150‑1 (MQ=255)
                                                                                               gCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTc                                                       <  1:1161964/150‑1 (MQ=255)
                                                                                                 gggTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTCGa                                                     >  2:3947801/1‑150 (MQ=255)
                                                                                                    tACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTCGAAGt                                                  <  2:2148644/150‑1 (MQ=255)
                                                                                                           aaTCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTCGAAGTTGACCAc                                           >  2:1980281/1‑150 (MQ=255)
                                                                                                           aaTCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTCGAAGTTGACCAc                                           >  1:2782291/1‑150 (MQ=255)
                                                                                                               aCTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTCGAAGTTGACCACGCgg                                       <  2:1649968/150‑1 (MQ=255)
                                                                                                                cTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTCGAAGTTGACCACGCGGt                                      >  2:812849/1‑150 (MQ=255)
                                                                                                                 tGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTCGAAGTTGACCACGCGGtt                                     <  1:3671009/150‑1 (MQ=255)
                                                                                                                  gATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTCGAAGTTGACCACGCGGTTa                                    >  1:185515/1‑150 (MQ=255)
                                                                                                                   aTTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTCGAAGTTGACCACGCGGTTAt                                   <  1:1972648/150‑1 (MQ=255)
                                                                                                                       cGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTCGAAGTTGACCACGCGGTTATTTTc                               <  1:2378/150‑1 (MQ=255)
                                                                                                                       cGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTCGAAGTTGACCACGCGGTTATTTTc                               <  1:3323/150‑1 (MQ=255)
                                                                                                                        ggTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTCGAAGTTGACCACGCGGTTATTTTcc                              >  1:2206375/1‑150 (MQ=255)
                                                                                                                        ggTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTCGAAGTTGACCACGCGGTTATTTTcc                              >  2:2898636/1‑150 (MQ=255)
                                                                                                                          tgtTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTCGAAGTTGACCACGCGGTTATTTTCCGt                            <  2:2399898/150‑1 (MQ=255)
                                                                                                                           gtTTCTGCCACGCGCAGACCCATTTAATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTCGAAGTTGACCACGCGGTTATTTTCCGTa                           <  2:4386383/150‑1 (MQ=255)
                                                                                                                            tttCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTCGAAGTTGACCACGCGGTTATTTTCCGTAc                          <  2:2222898/150‑1 (MQ=255)
                                                                                                                             ttCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTCGAAGTTGACCACGCGGTTATTTTCCGTACg                         <  2:3339921/150‑1 (MQ=255)
                                                                                                                               cTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTCGAAGTTGACCACGCGGTTATTTTCCGTACGcc                       >  2:3281972/1‑150 (MQ=255)
                                                                                                                                tGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTCGAAGTTGACCACGCGGTTATTTTCCGTACGccc                      >  1:627682/1‑150 (MQ=255)
                                                                                                                                tGCCACGCGCAGACCCATTTAATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTCGAAGTTGACCACGCGGTTATTTTCCGTACGccc                      >  1:3025279/1‑150 (MQ=255)
                                                                                                                                  ccACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTCGAAGTTGACCACGCGGTTATTTTCCGTACGCCCgg                    >  2:943075/1‑150 (MQ=255)
                                                                                                                                   cACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTCGAAGTTGACCACGCGGTTATTTTCCGTACGCCCGGa                   <  1:3817455/150‑1 (MQ=255)
                                                                                                                                   cACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTCGAAGTTGACCACGCGGTTATTTTCCGTACGCCCGGa                   <  2:190942/150‑1 (MQ=255)
                                                                                                                                         cAGACCCATTTAATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTCGAAGTTGACCACGCGGTTATTTTCCGTACGCCCGGAAAGCTc             >  1:523692/1‑150 (MQ=255)
                                                                                                                                          aGACCCATTTAATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTCGAAGTTGACCACGCGGTTATTTTCCGTACGCCCGGAAAGCTcc            <  1:1195357/150‑1 (MQ=255)
                                                                                                                                           gACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTCGAAGTTGACCACGCGGTTATTTTCCGTACGCCCGGAAAGCTCCa           <  1:73642/150‑1 (MQ=255)
                                                                                                                                           gACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTCGAAGTTGACCACGCGGTTATTTTCCGTACGCCCGGAAAGCTCCa           >  2:4629701/1‑150 (MQ=255)
                                                                                                                                             cccATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTCGAAGTTGACCACGCGGTTATTTTCCGTACGCCCGGAAAGCTCCatg         >  1:1498013/1‑150 (MQ=255)
                                                                                                                                             cccATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTCGAAGTTGACCACGCGGTTATTTTCCGTACGCCCGGAAAGCTCCatg         >  2:1502358/1‑150 (MQ=255)
                                                                                                                                                aTTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTCGAAGTTGACCACGCGGTTATTTTCCGTACGCCCGGAAAGCTCCatgatg      >  1:1626039/1‑150 (MQ=255)
                                                                                                                                                aTTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTCGAAGTTGACCACGCGGTTATTTTCCGTACGCCCGGAAAGCTCCatgatg      >  1:1540158/1‑150 (MQ=255)
                                                                                                                                                  ttCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTCGAAGTTGACCACGCGGTTATTTTCCGTACGCCCGGAAAGCTCCATGATGct    <  1:2662933/150‑1 (MQ=255)
                                                                                                                                                   tCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTCGAAGTTGACCACGCGGTTATTTTCCGTACGCCCGGAAAGCTCCATGATGctc   <  2:4134748/150‑1 (MQ=255)
                                                                                                                                                   tCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTCGAAGTTGACCACGCGGTTATTTTCCGTACGCCCGGAAAGCTCCATGATGctc   >  2:3214135/1‑150 (MQ=255)
                                                                                                                                                   tCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTCGAAGTTGACCACGCGGTTATTTTCCGTACGCCCGGAAAGCTCCATGATGctc   <  1:100595/150‑1 (MQ=255)
                                                                                                                                                    cATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTCGAAGTTGACCACGCGGTTATTTTCCGTACGCCCGGAAAGCTCCATGATGctct  <  1:1340223/150‑1 (MQ=255)
                                                                                                                                                    cATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTCGAAGTTGACCACGCGGTTATTTTCCGTACGCCCGGAAAGCTCCATGATGctct  <  1:747779/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AAAGATGGGGATAAAGAGAGAAAAAACAAGGCCCACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTCGAAGTTGACCACGCGGTTATTTTCCGTACGCCCGGAAAGCTCCATGATGCTCT  >  CP009273/688931‑689228

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