Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,880,202 G→T 9.8% 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 9.8% 755.4 / 44.3 214S126* (TCA→TAA) ygcBR‑loop helicase‑annealase Cas3 needed for Cascade anti‑viral activity
Reads supporting (aligned to +/- strand):  ref base G (86/107);  new base T (10/11);  total (96/118)
Fisher's exact test for biased strand distribution p-value = 8.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 216 total aligned reads displayed.

TTTCCCAACGCGACTTATCTTGATCCTGGGAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGATTTATATTGGAATCGTATATCAAACTTTCCA  >  CP009273/2880053‑2880351
                                                                                                                                                     |                                                                                                                                                     
 ttCCCAACGCGACTTATCTTGATCCTGGGAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGa                                                                                                                                                      <  2:1217992/150‑1 (MQ=255)
 ttCCCAACGCGACTTATCTTGATCCTGGGAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGa                                                                                                                                                      <  1:130708/150‑1 (MQ=255)
    ccAACGCGACTTATCTTGATCCTGGGAATGTAATATAAAACCCTGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTTAATc                                                                                                                                                   <  2:2416809/150‑1 (MQ=255)
      aaCGCGACTTATCTTGATCCTGGGAGTGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCt                                                                                                                                                 <  2:169939/150‑1 (MQ=255)
      aaCGCGACTTATCTTGATCCTGGGAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCt                                                                                                                                                 >  2:2309413/1‑150 (MQ=255)
        cgcgACTTATCTTGATCCTGGGAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTgg                                                                                                                                               >  1:3495651/1‑150 (MQ=255)
         gcgACTTATCTTGATCCTGGGAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGt                                                                                                                                              <  1:2445536/150‑1 (MQ=255)
           gACTTATCTTGATCCTGGGAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTa                                                                                                                                            >  1:1480535/1‑150 (MQ=255)
           gACTTATCTTGATCCTGGGAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTa                                                                                                                                            >  1:182396/1‑150 (MQ=255)
            aCTTATCTTGATCCTGGGAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTTAATCCTGGTTaa                                                                                                                                           >  2:3558822/1‑150 (MQ=255)
            aCTTATCTTGATCCTGGGAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTaa                                                                                                                                           <  2:3803552/150‑1 (MQ=255)
             cTTATCTTGATCCTGGGAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTaaa                                                                                                                                          <  2:1679690/150‑1 (MQ=255)
              ttATCTTGATCCTGGGAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAAc                                                                                                                                         <  1:3104060/150‑1 (MQ=255)
                aTCTTGATCCTGGGAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCa                                                                                                                                       <  2:182648/150‑1 (MQ=255)
                aTCTTGATCCTGGGAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCa                                                                                                                                       >  2:392097/1‑150 (MQ=255)
                aTCTTGATCCTGGGAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCa                                                                                                                                       >  1:4182562/1‑150 (MQ=255)
                 tCTTGATCCTGGGAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCaa                                                                                                                                      <  2:964015/150‑1 (MQ=255)
                     gATCCTGGGAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTTAACCAATAct                                                                                                                                  <  2:1725189/150‑2 (MQ=255)
                      aTCCTGGGAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAg                                                                                                                                 >  2:2555070/1‑150 (MQ=255)
                         cTGGGAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGAcc                                                                                                                              >  2:638988/1‑150 (MQ=255)
                         cTGGGAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGAcc                                                                                                                              <  1:113096/150‑1 (MQ=255)
                          tGGGAATGTAATATATCTCCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCg                                                                                                                             <  1:3660204/150‑1 (MQ=255)
                           gggAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCgg                                                                                                                            >  2:1818707/1‑150 (MQ=255)
                             gAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGACTCCTGGTTAAACCAATACAGACCGGCt                                                                                                                          >  1:2166711/1‑150 (MQ=255)
                             gAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCt                                                                                                                          <  2:4296129/150‑1 (MQ=255)
                                tGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCa                                                                                                                       <  1:75300/150‑1 (MQ=255)
                                      ataAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGa                                                                                                                 >  1:49908/1‑150 (MQ=255)
                                      ataAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGa                                                                                                                 <  1:4265053/150‑1 (MQ=255)
                                       taAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGAt                                                                                                                <  2:2829338/150‑1 (MQ=255)
                                       taAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGAt                                                                                                                >  2:2977751/1‑150 (MQ=255)
                                           aCCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAATTCCGCGAGAGTTTGCTCTGAAAGTGAATCATGGTTTAACCAATACATACCGGCTGCACCATGATTaaa                                                                                                            <  2:3966083/150‑1 (MQ=255)
                                            ccATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAAt                                                                                                           >  1:3314781/1‑150 (MQ=255)
                                             cATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCCTGATTAAAtt                                                                                                          >  1:3294897/1‑150 (MQ=255)
                                              atgatgTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTTAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAAttt                                                                                                         >  1:2696286/1‑150 (MQ=255)
                                              atgatgTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAAttt                                                                                                         <  1:77861/150‑1 (MQ=255)
                                               tgatgTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTa                                                                                                        >  2:3452284/1‑150 (MQ=255)
                                                gatgTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTAc                                                                                                       <  1:3022305/150‑1 (MQ=255)
                                                gatgTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTAc                                                                                                       <  1:435511/150‑1 (MQ=255)
                                                 atgTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACg                                                                                                      <  2:2711663/150‑1 (MQ=255)
                                                      cTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGcaca                                                                                                 <  2:1341315/150‑1 (MQ=255)
                                                         tAACGGCCTCTACCCATGGAAAACAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATTATTAAATTTACGGCACAttt                                                                                              <  2:1881655/150‑1 (MQ=255)
                                                          aaCGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTtg                                                                                             >  2:318980/1‑150 (MQ=255)
                                                           aCGGCCAATACCCATTTAAACCAGGACTCATAATGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTtgt                                                                                            <  2:3410481/150‑1 (MQ=255)
                                                            cGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTTAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTtgtg                                                                                           >  2:3450419/1‑150 (MQ=255)
                                                            cGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTtgtg                                                                                           <  2:1451488/150‑1 (MQ=255)
                                                            cGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTtgtg                                                                                           <  2:4178429/150‑1 (MQ=255)
                                                                ctctACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGa                                                                                       >  1:331500/1‑150 (MQ=255)
                                                                 tctACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGAt                                                                                      <  1:1973638/150‑1 (MQ=255)
                                                                  ctACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTTAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATg                                                                                     <  2:412586/150‑1 (MQ=255)
                                                                     cccATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGAc                                                                                  >  2:3850657/1‑150 (MQ=255)
                                                                       cATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCa                                                                                <  1:1855908/150‑1 (MQ=255)
                                                                       cATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCa                                                                                <  1:3361228/150‑1 (MQ=255)
                                                                        aTGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCAt                                                                               <  1:3113016/150‑1 (MQ=255)
                                                                             aaCCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGt                                                                          >  2:2785339/1‑150 (MQ=255)
                                                                              aCCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACAGTTGTGTTGATGGACCATTAAGTg                                                                         >  2:791953/1‑150 (MQ=255)
                                                                               ccAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGa                                                                        <  1:2340434/150‑1 (MQ=255)
                                                                                cAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTTAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGAt                                                                       <  2:745033/150‑1 (MQ=255)
                                                                                cAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGAt                                                                       >  2:2794915/1‑150 (MQ=255)
                                                                                  ggACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATgg                                                                     >  1:3150249/1‑150 (MQ=255)
                                                                                   gACTCATAAGGATGAGGAGCGGCAACAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTTAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGc                                                                    <  1:432792/150‑1 (MQ=255)
                                                                                     cTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGt                                                                  >  2:2954419/1‑150 (MQ=255)
                                                                                        aTAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGAGTTGc                                                               <  2:2888643/150‑1 (MQ=255)
                                                                                           aGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAgg                                                            <  1:180480/150‑1 (MQ=255)
                                                                                           aGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAgg                                                            <  2:450576/150‑1 (MQ=255)
                                                                                                 aggagCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTaa                                                      <  1:2116755/150‑1 (MQ=255)
                                                                                                  ggagCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAAt                                                     <  1:535091/150‑1 (MQ=255)
                                                                                                    agCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAAttt                                                   <  2:685825/150‑1 (MQ=255)
                                                                                                    agCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATATGTGTTGATTGACCATTAAGTTATGGCGTTGCAGGATTTAAatt                                                   >  1:1265835/1‑150 (MQ=255)
                                                                                                      cGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCa                                                 >  2:1316716/1‑150 (MQ=255)
                                                                                                       ggCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTTAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAg                                                >  2:730081/1‑150 (MQ=255)
                                                                                                       ggCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAg                                                >  1:2219145/1‑150 (MQ=255)
                                                                                                       ggCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAg                                                <  2:2112263/150‑1 (MQ=255)
                                                                                                       ggCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAg                                                <  2:306894/150‑1 (MQ=255)
                                                                                                         cATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGcc                                              <  1:2351406/150‑1 (MQ=255)
                                                                                                          aTCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCa                                             >  1:1796588/1‑150 (MQ=255)
                                                                                                             aaaaaaaCTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACt                                          >  2:110820/1‑150 (MQ=255)
                                                                                                             aaaaaaaCTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACt                                          >  2:2636961/1‑150 (MQ=255)
                                                                                                             aaaaaaaCTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACt                                          <  1:1765709/150‑1 (MQ=255)
                                                                                                             aaaaaaaCTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACt                                          <  1:898410/150‑1 (MQ=255)
                                                                                                              aaaaaaCTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACtt                                         >  1:1157767/1‑150 (MQ=255)
                                                                                                               aaaaaCTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACttt                                        >  1:4159246/1‑150 (MQ=255)
                                                                                                               aaaaaCTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATGTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACtgt                                        >  1:2818830/1‑148 (MQ=255)
                                                                                                                aaaaCTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTc                                       >  2:2094081/1‑150 (MQ=255)
                                                                                                                   aCTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTctgc                                    <  1:2780279/150‑1 (MQ=255)
                                                                                                                   aCTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTctgc                                    <  2:4193581/150‑1 (MQ=255)
                                                                                                                   aCTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATGTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTctgc                                    >  2:4029419/1‑150 (MQ=255)
                                                                                                                    cTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATCAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTctgct                                   <  2:1186118/150‑1 (MQ=255)
                                                                                                                      gAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGa                                 <  2:1282026/150‑1 (MQ=255)
                                                                                                                       aaaaaaTCCCCGAGAGAGTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGAt                                <  1:3742910/150‑1 (MQ=255)
                                                                                                                        aaaaaTCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGAtt                               >  1:652797/1‑150 (MQ=255)
                                                                                                                           aaTCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGATTtat                            <  1:1138541/150‑1 (MQ=255)
                                                                                                                            aTCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGATTtata                           <  2:1021488/150‑1 (MQ=255)
                                                                                                                             tCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGATTTatat                          <  1:162842/150‑1 (MQ=255)
                                                                                                                                ccGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGATTTATATTgg                       >  2:2724009/1‑150 (MQ=255)
                                                                                                                                 cGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGATTTATATTGGa                      >  1:388224/1‑150 (MQ=255)
                                                                                                                                  gagagaTTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGATTTATATTGGaa                     >  1:2259868/1‑150 (MQ=255)
                                                                                                                                        ttGCTCTGAAAGTTAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGATTTATATTGGAATCGtat               <  1:3450419/150‑1 (MQ=255)
                                                                                                                                         tGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGATTTATATTGGAATCGtata              >  1:136841/1‑150 (MQ=255)
                                                                                                                                          gCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGATTTATATTGGAATCGTatat             >  1:1014825/1‑150 (MQ=255)
                                                                                                                                                 aaGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGATTTATATTGGAATCGTATATCAAACtt      <  1:3606649/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTTCCCAACGCGACTTATCTTGATCCTGGGAATGTAATATAAAACCATGATGTCCTGTAACGGCCTCTACCCATGGAAACCAGGACTCATAAGGATGAGGAGCGGCATCAAAAAAACTGAAAAAATCCCCGAGAGATTGCTCTGAAAGTGAATCCTGGTTAAACCAATACAGACCGGCTGCACCATGATTAAATTTACGGCACATTTGTGTTGATGGACCATTAAGTGATGGCGTTGCAGGATTTAATTTCAGCCAACTTTCTGCTGATTTATATTGGAATCGTATATCAAACTTTCCA  >  CP009273/2880053‑2880351

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: