Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,880,202 G→T 6.5% S126* (TCA→TAA)  ygcB ← R‑loop helicase‑annealase Cas3 needed for Cascade anti‑viral activity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,880,2020GT 6.5% 1034.7 / 32.2 269S126* (TCA→TAA) ygcBR‑loop helicase‑annealase Cas3 needed for Cascade anti‑viral activity
Reads supporting (aligned to +/- strand):  ref base G (113/138);  new base T (7/11);  total (120/149)
Fisher's exact test for biased strand distribution p-value = 8.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 270 total aligned reads displayed.

TTTCCCAACGCGACTTATCTTGATCCTGGGAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGATTTATATTGGAATCGTATATCAAACTTTCC  >  CP009273/2880053‑2880350
                                                                                                                                                      |                                                                                                                                                    
 ttCCCAACGCGACTTATCTTGATCCTGTGAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGa                                                                                                                                                     <  1:4617456/150‑1 (MQ=255)
 ttCCCAACGCGACTTATCTTGATCCTGGGAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGa                                                                                                                                                     >  1:1229532/1‑150 (MQ=255)
 ttCCCAACGCGACTTATCTTGATCCTGGGAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGa                                                                                                                                                     <  1:1364379/150‑1 (MQ=255)
 ttCCCAACGCGACTTATCTTGATCCTGGGAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGa                                                                                                                                                     <  2:3300967/150‑1 (MQ=255)
 ttCCCAACGCGACTTATCTTGATCCTGGGAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGa                                                                                                                                                     >  1:2412252/1‑150 (MQ=255)
   cccAACGCGACTTATCTTGATCCTGGGAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAAt                                                                                                                                                   >  2:748494/1‑150 (MQ=255)
      aaCGCGACTTATCTTGATCCTGGGAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCt                                                                                                                                                <  2:1629158/150‑1 (MQ=255)
      aaCGCGACTTATCTTGATCCTGGGAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCt                                                                                                                                                >  2:218538/1‑150 (MQ=255)
       aCGCGACTTATCTTGATCCTGGGAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTg                                                                                                                                               <  1:4124211/150‑1 (MQ=255)
       aCGCGACTTATCTTGATCCTGGGAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTg                                                                                                                                               <  2:4221855/150‑1 (MQ=255)
         gcgACTTATCTTGATCCTGGGAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGt                                                                                                                                             <  2:1036233/150‑1 (MQ=255)
          cgACTTATCTTGATCCTGGGAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGtt                                                                                                                                            <  1:1170229/150‑1 (MQ=255)
           gACTTATCTTGATCCTGGGAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTa                                                                                                                                           >  2:1466814/1‑150 (MQ=255)
           gACTTATCTTGATCCTGGGAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTa                                                                                                                                           >  1:394924/1‑150 (MQ=255)
             cTTATCTTGATCCTGGGAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTaaa                                                                                                                                         >  1:3032032/1‑150 (MQ=255)
             cTTATCTTGATCCTGGGAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTaaa                                                                                                                                         >  1:643700/1‑150 (MQ=255)
                  cTTGATCCTGGGAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTTAATCCTGGTTAAACCAAt                                                                                                                                    <  1:4153643/150‑1 (MQ=255)
                      aTCCTGGGAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAg                                                                                                                                >  1:4034757/1‑150 (MQ=255)
                      aTCCTGGGAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAg                                                                                                                                <  2:3178815/150‑1 (MQ=255)
                        ccTGGGAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGAc                                                                                                                              >  2:4696758/1‑150 (MQ=255)
                         cTGGGAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGAcc                                                                                                                             <  2:1514176/150‑1 (MQ=255)
                         cTGGGAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGAcc                                                                                                                             <  2:3046254/150‑1 (MQ=255)
                           gggAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCgg                                                                                                                           >  1:690183/1‑150 (MQ=255)
                              aaTGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTg                                                                                                                        >  1:3248893/1‑150 (MQ=255)
                                tGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCa                                                                                                                      <  1:887146/150‑1 (MQ=255)
                                   aaTATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCa                                                                                                                   >  2:3481995/1‑150 (MQ=255)
                                   aaTATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCa                                                                                                                   >  2:1913969/1‑150 (MQ=255)
                                         aaaCCATGATGTCCTGTAACGGCCTCTACCCATTGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTa                                                                                                             <  2:4640290/150‑1 (MQ=255)
                                         aaaCCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTa                                                                                                             <  2:2371365/150‑1 (MQ=255)
                                          aaCCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTaa                                                                                                            <  1:3826078/150‑1 (MQ=255)
                                            ccATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAAt                                                                                                          <  1:15472/150‑1 (MQ=255)
                                             cATTATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCTGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAAtt                                                                                                         <  2:1343906/150‑1 (MQ=255)
                                             cATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAAtt                                                                                                         >  2:41339/1‑150 (MQ=255)
                                              atgatgTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAAttt                                                                                                        >  1:2021032/1‑150 (MQ=255)
                                              atgatgTCCTGTAACGGCCTCTACACATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAACAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAAACAATACAGACCGGCTGCACCATGATTAAAttt                                                                                                        <  2:3996270/150‑1 (MQ=255)
                                               tgatgTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTa                                                                                                       >  2:3580607/1‑150 (MQ=255)
                                                gatgTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTAc                                                                                                      <  2:1085525/150‑1 (MQ=255)
                                                 atgTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACg                                                                                                     >  1:512831/1‑150 (MQ=255)
                                                   gTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGc                                                                                                   >  1:1463578/1‑150 (MQ=255)
                                                    tCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGca                                                                                                  >  1:1902032/1‑150 (MQ=255)
                                                     ccTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGcac                                                                                                 >  1:3332336/1‑150 (MQ=255)
                                                      cTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGcaca                                                                                                <  2:1269860/150‑1 (MQ=255)
                                                      cTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGcaca                                                                                                >  1:2638923/1‑150 (MQ=255)
                                                       tGTAACGTCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAAACTTAAAAAAACCCCGAGAGATTGATATGAAATTTAAGCATGGTTAAACAAAAAAATACCGGAGCAACAATGATAAAATTTACGGcaca                                                                                                >  2:790314/1‑150 (MQ=255)
                                                         tAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACAttt                                                                                             <  2:4314941/150‑1 (MQ=255)
                                                            cGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTtgtg                                                                                          >  2:2021499/1‑150 (MQ=255)
                                                             ggCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCCAAACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTgtgt                                                                                         >  1:668851/1‑150 (MQ=255)
                                                              gCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGtt                                                                                        <  2:1817545/150‑1 (MQ=255)
                                                               ccTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTg                                                                                       <  2:1778922/150‑1 (MQ=255)
                                                               ccTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTg                                                                                       >  2:4120863/1‑150 (MQ=255)
                                                                ctctACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGa                                                                                      >  2:3059523/1‑150 (MQ=255)
                                                                 tctACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGAt                                                                                     >  1:2847779/1‑150 (MQ=255)
                                                                     cccATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGAc                                                                                 >  1:428623/1‑150 (MQ=255)
                                                                                cAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGAt                                                                      >  2:2619911/1‑150 (MQ=255)
                                                                                   gACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGc                                                                   <  1:4102323/150‑1 (MQ=255)
                                                                                    aCTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCg                                                                  <  1:3565406/150‑1 (MQ=255)
                                                                                      tCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGtt                                                                <  2:689202/150‑1 (MQ=255)
                                                                                       cATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTg                                                               <  2:865628/150‑1 (MQ=255)
                                                                                           aGGATGAGGAGCGGCATCAAAAAAACTGAAACAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAgg                                                           <  2:2233619/150‑1 (MQ=255)
                                                                                           aGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAgg                                                           <  2:2233463/150‑1 (MQ=255)
                                                                                            ggATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGa                                                          >  2:854381/1‑150 (MQ=255)
                                                                                               tGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGAttt                                                       >  1:1691401/1‑150 (MQ=255)
                                                                                                 aggagCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTaa                                                     <  1:874553/150‑1 (MQ=255)
                                                                                                    agCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAAttt                                                  >  1:3888167/1‑150 (MQ=255)
                                                                                                      cGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCa                                                <  1:361851/150‑1 (MQ=255)
                                                                                                      cGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCa                                                <  2:4523009/150‑1 (MQ=255)
                                                                                                       ggCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAg                                               >  2:3332716/1‑150 (MQ=255)
                                                                                                               aaaaaaCTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCTGCCAACtt                                        >  1:1729618/1‑150 (MQ=255)
                                                                                                                   aaCTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTctg                                    <  1:899739/150‑1 (MQ=255)
                                                                                                                     cTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTctgct                                  <  1:578139/150‑1 (MQ=255)
                                                                                                                     cTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTctgct                                  >  1:3997125/1‑150 (MQ=255)
                                                                                                                      tGAAAAAATCCCCGAGAGATTGCTCTGAAAGTTAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTctgctg                                 <  1:1023133/150‑1 (MQ=255)
                                                                                                                       gAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGa                                <  1:2746858/150‑1 (MQ=255)
                                                                                                                       gAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGa                                <  1:277504/150‑1 (MQ=255)
                                                                                                                        aaaaaaTCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTTCTGAt                               <  2:2026941/150‑1 (MQ=255)
                                                                                                                        aaaaaaTCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGAt                               <  2:258758/150‑1 (MQ=255)
                                                                                                                         aaaaaTCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGAtt                              <  2:110700/150‑1 (MQ=255)
                                                                                                                           aaaTCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGATTta                            <  2:1155230/150‑1 (MQ=255)
                                                                                                                           aaaTCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGATTta                            >  2:4528234/1‑150 (MQ=255)
                                                                                                                            aaTCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGATTtat                           <  2:3771197/150‑1 (MQ=255)
                                                                                                                             aTCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGATTtata                          >  1:3729074/1‑150 (MQ=255)
                                                                                                                              tCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGATTTatat                         >  2:3543063/1‑150 (MQ=255)
                                                                                                                               ccccGAGAGATTGCTCTGAAAGTTAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGATTTATAtt                        <  1:2292358/150‑1 (MQ=255)
                                                                                                                                    agagaTTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGATTTATATTGGAAt                   <  2:3836950/150‑1 (MQ=255)
                                                                                                                                         ttGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGATTTATATTGGAATCGtat              <  1:371206/150‑1 (MQ=255)
                                                                                                                                         ttGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGATTTATATTGGAATCGtat              >  2:1015581/1‑150 (MQ=255)
                                                                                                                                          tGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGATTTATATTGGAATCGtata             <  1:219797/150‑1 (MQ=255)
                                                                                                                                          tGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGATTTATATTGGAATCGtata             <  1:218538/150‑1 (MQ=255)
                                                                                                                                          tGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGATTTATATTGGAATCGtata             <  2:2491398/150‑1 (MQ=255)
                                                                                                                                          tGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGATTTATATTGGAATCGtata             <  2:579004/150‑1 (MQ=255)
                                                                                                                                           gCTCTGATAGTGAATCCTGGTTAAACGAATACAGACCGGCTGCACCATGATTAAATTTTCGGCACATTTGTGTTGATGGACCATTATGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGATTTATATTGGAATCGTatat            <  2:2887995/150‑1 (MQ=255)
                                                                                                                                           gCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCTTGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGATTTATATTGGAATCGTatat            <  2:1827301/150‑1 (MQ=255)
                                                                                                                                           gCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGATTTATATTGGAATCGTatat            >  1:148337/1‑150 (MQ=255)
                                                                                                                                             tctGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTTATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGATTTATATTGGAATCGTATATCa          >  1:3942997/1‑150 (MQ=255)
                                                                                                                                               tGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGATTTATATTGGAATCGTATATCaaa        >  2:1830714/1‑150 (MQ=255)
                                                                                                                                                 aaaGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGATTTATATTGGAATCGTATATCAAACt      <  1:4707903/150‑1 (MQ=255)
                                                                                                                                                 aaaGTGAATCCTGGTTAAACCAAGACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGATTTATATTGGAATCGTATATCAAACt      <  1:1052139/150‑1 (MQ=255)
                                                                                                                                                    gTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGATTTATATTGGAATCGTATATCAAACTTTc   <  2:3888167/150‑1 (MQ=255)
                                                                                                                                                    gTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATATGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGATTTATATTGGAATCGTATATCAAACTTTc   <  2:3889649/150‑1 (MQ=255)
                                                                                                                                                     tGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGATTTATATTGGAATCGTATATCAAACTTTcc  >  1:1080845/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
TTTCCCAACGCGACTTATCTTGATCCTGGGAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGATTTATATTGGAATCGTATATCAAACTTTCC  >  CP009273/2880053‑2880350

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