Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,884,577 T→C 6.1% R881R (AGA→AGG) ‡ ygcB ← Cascade complex anti‑viral R‑loop helicase‑annealase Cas3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,884,5770TC 6.1% 530.1 / 17.2 181R881R (AGA→AGG) ‡ygcBCascade complex anti‑viral R‑loop helicase‑annealase Cas3
Reads supporting (aligned to +/- strand):  ref base T (80/90);  new base C (3/8);  total (83/98)
Fisher's exact test for biased strand distribution p-value = 2.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 184 total aligned reads displayed.

ACTGTGAATATAAAATCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATT  >  NC_000913/2884428‑2884726
                                                                                                                                                     |                                                                                                                                                      
    tgAATATAAAATCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTgg                                                                                                                                                    <  2:358112/150‑1 (MQ=255)
     taaTATAAAATCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGt                                                                                                                                                   <  2:678618/149‑1 (MQ=255)
      aaTATAAAATCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTc                                                                                                                                                  <  2:812680/150‑1 (MQ=255)
         ataAAATCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTTCTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATc                                                                                                                                               <  1:391005/150‑1 (MQ=255)
         ataAAATCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATc                                                                                                                                               <  2:942121/150‑1 (MQ=255)
         ataAAATCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCATACTTAAGTAGGGATAAAACGTTATTGCTATAATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGTGATGACTCTGGTCATc                                                                                                                                               >  1:291158/1‑150 (MQ=255)
           aaaaTCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACCAGAGTCATccc                                                                                                                                             <  2:549718/150‑1 (MQ=255)
           aaaaTCTCATACCGGGAAATTAAAACAAGATGTACATTGAGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGGCTTATTATAGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATccc                                                                                                                                             <  2:493453/150‑1 (MQ=255)
            aaaTCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCt                                                                                                                                            >  2:71978/1‑150 (MQ=255)
            aaaTCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCt                                                                                                                                            >  1:1662289/1‑150 (MQ=255)
            aaaTCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCt                                                                                                                                            <  1:1228832/150‑1 (MQ=255)
            aaaTCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACCAGAGTCATCCCt                                                                                                                                            >  1:1569287/1‑150 (MQ=255)
             aaTCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCtt                                                                                                                                           <  1:600137/150‑1 (MQ=255)
             aaTCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACCAGAGTCATCCCtg                                                                                                                                           <  1:1395693/150‑2 (MQ=255)
             aaTCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACCAGAGTCATCCCtg                                                                                                                                           <  2:1009731/150‑2 (MQ=255)
               tctcATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCa                                                                                                                                         <  2:728594/150‑1 (MQ=255)
                 tcATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACCAGAGTCATCCCTGCaaa                                                                                                                                       <  2:1322024/150‑3 (MQ=255)
                      ccGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTg                                                                                                                                  <  1:1076682/150‑1 (MQ=255)
                      ccGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATATGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTg                                                                                                                                  <  2:535859/150‑1 (MQ=255)
                      ccGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATATGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTg                                                                                                                                  <  2:755919/150‑1 (MQ=255)
                      ccGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAAACATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGTTCATCACTTCATCCCCTg                                                                                                                                  >  2:209330/1‑150 (MQ=255)
                         ggAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTACGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGtat                                                                                                                               <  1:933506/150‑1 (MQ=255)
                         ggAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGtat                                                                                                                               <  1:740855/150‑1 (MQ=255)
                          gAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGCGATTTGCAGGGATGAATCTGGTCATCCCTTCATCACCTGtata                                                                                                                              >  2:1021945/1‑150 (MQ=255)
                           aaaTTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAg                                                                                                                             <  1:1021513/150‑1 (MQ=255)
                             aTTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGt                                                                                                                           >  1:1198184/1‑150 (MQ=255)
                              ttAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGta                                                                                                                          >  2:1385983/1‑150 (MQ=255)
                                aaaaGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCAGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCATTCATCCACTGTATAGGtaat                                                                                                                        >  2:1552915/1‑150 (MQ=255)
                                     aagaTGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACaa                                                                                                                   >  2:1460956/1‑150 (MQ=255)
                                      agaTGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAAt                                                                                                                  >  1:1718466/1‑150 (MQ=255)
                                         tGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACt                                                                                                               >  1:1606604/1‑150 (MQ=255)
                                         tGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACt                                                                                                               >  2:139194/1‑150 (MQ=255)
                                              aTTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTAc                                                                                                          <  1:709459/150‑1 (MQ=255)
                                              aTTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTAc                                                                                                          >  2:1042614/1‑150 (MQ=255)
                                                tgtgCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTAccc                                                                                                        <  2:1207077/150‑1 (MQ=255)
                                                    cACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGcc                                                                                                   >  2:600529/1‑150 (MQ=255)
                                                    cACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGcc                                                                                                   >  2:1602256/1‑150 (MQ=255)
                                                      ccTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAg                                                                                                 <  1:784650/150‑1 (MQ=255)
                                                        atCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGAc                                                                                               >  2:266516/2‑150 (MQ=255)
                                                         tCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGAcc                                                                                              <  1:871966/150‑1 (MQ=255)
                                                          cccTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGAccc                                                                                             >  1:109794/1‑150 (MQ=255)
                                                          cccTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGAccc                                                                                             >  1:1028015/1‑150 (MQ=255)
                                                              aCTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACccatcc                                                                                         <  1:723235/150‑1 (MQ=255)
                                                               cTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACccatcca                                                                                        <  1:690305/150‑1 (MQ=255)
                                                               cTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACccatcca                                                                                        <  2:1000830/150‑1 (MQ=255)
                                                                 tAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCcatccatc                                                                                      >  2:1240192/1‑150 (MQ=255)
                                                                   aGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGTTATTAACAAACAGTTTCCCATGCAAGACCAaaccaaccc                                                                                    >  2:847310/1‑149 (MQ=255)
                                                                   aGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCAtccatcca                                                                                    <  1:1551849/150‑1 (MQ=255)
                                                                        ggATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTc                                                                               <  1:344874/150‑1 (MQ=255)
                                                                        ggATAAAACGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTc                                                                               <  2:1604530/150‑1 (MQ=255)
                                                                          aTAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTg                                                                             >  2:1258568/1‑150 (MQ=255)
                                                                            aaaCCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGtt                                                                           >  1:655254/1‑150 (MQ=255)
                                                                             aaCCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCctgtttctgttt                                                                          >  1:1391594/1‑150 (MQ=255)
                                                                             aaCCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGttt                                                                          >  1:500869/1‑150 (MQ=255)
                                                                             aaCCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTATGttt                                                                          >  2:339229/1‑150 (MQ=255)
                                                                              aCCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTc                                                                         >  2:916775/1‑150 (MQ=255)
                                                                               ccGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTcc                                                                        >  2:1624772/1‑150 (MQ=255)
                                                                                       ggTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGc                                                                >  2:473650/1‑150 (MQ=255)
                                                                                       ggTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGc                                                                >  1:1183504/1‑150 (MQ=255)
                                                                                        gTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGcc                                                               >  1:414306/1‑150 (MQ=255)
                                                                                         tCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCa                                                              >  2:422109/1‑150 (MQ=255)
                                                                                          cTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCaa                                                             >  1:1199842/1‑150 (MQ=255)
                                                                                            tattatCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAg                                                           >  1:971343/1‑150 (MQ=255)
                                                                                                atCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAAc                                                       <  1:769289/150‑1 (MQ=255)
                                                                                                atCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAAc                                                       <  1:840107/150‑1 (MQ=255)
                                                                                                 tCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAAcc                                                      >  1:1524646/1‑150 (MQ=255)
                                                                                                 tCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAAcc                                                      <  2:1086689/150‑1 (MQ=255)
                                                                                                 tCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAAcc                                                      >  1:1570845/1‑150 (MQ=255)
                                                                                                      aTTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCt                                                 >  1:1293692/1‑150 (MQ=255)
                                                                                                       ttGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCtt                                                <  1:219031/150‑1 (MQ=255)
                                                                                                          aTAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAGCAATACTGTTACCCTGCCAGACACATCCATCCAGATTCTGTTTCACTTCCAGCCAAAGTAACCCATCTTCAt                                             >  1:1423323/1‑150 (MQ=255)
                                                                                                           tAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATc                                            >  2:628473/1‑150 (MQ=255)
                                                                                                            aaCAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATct                                           <  1:1385983/150‑1 (MQ=255)
                                                                                                             aCAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATcta                                          >  1:427580/1‑150 (MQ=255)
                                                                                                              cAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATctac                                         <  1:600529/150‑1 (MQ=255)
                                                                                                               aaTCATTCCCGAAGTTATTTGGGATTTGCAGGGATTACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATctact                                        <  2:55940/150‑1 (MQ=255)
                                                                                                               aaTCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATctact                                        <  2:1198184/150‑1 (MQ=255)
                                                                                                                    ttCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTc                                   <  1:339096/150‑1 (MQ=255)
                                                                                                                     tCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCa                                  <  1:1602256/150‑1 (MQ=255)
                                                                                                                     tCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCa                                  >  1:637128/1‑150 (MQ=255)
                                                                                                                     tCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCa                                  >  1:643245/1‑150 (MQ=255)
                                                                                                                      cccGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAg                                 <  2:1570845/150‑1 (MQ=255)
                                                                                                                          aaGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGaaaa                             <  1:422109/150‑1 (MQ=255)
                                                                                                                            gTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAAc                           <  2:1661778/150‑1 (MQ=255)
                                                                                                                            gTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAAc                           <  1:473650/150‑1 (MQ=255)
                                                                                                                              tATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTa                         <  1:101893/150‑1 (MQ=255)
                                                                                                                                  tGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCCTCCCTCCAGATTCTGTTTCCCTTCCAGCCAAAGTCACCCCTCTTCCTCTACTACTTCAGAAAAACTACggt                     >  2:276591/1‑148 (MQ=255)
                                                                                                                                   gggATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGttt                    <  1:1460956/150‑1 (MQ=255)
                                                                                                                                        ttGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAgg               <  2:109794/150‑1 (MQ=255)
                                                                                                                                        ttGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAgg               <  1:266516/150‑1 (MQ=255)
                                                                                                                                         tGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGt              <  1:274466/150‑1 (MQ=255)
                                                                                                                                          gCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGCTTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTg             >  1:493192/1‑150 (MQ=255)
                                                                                                                                           cAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGAAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGa            >  2:225702/1‑150 (MQ=255)
                                                                                                                                            aGGGATGACCAGAGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGaa           >  2:1686625/1‑150 (MQ=255)
                                                                                                                                              ggATGACTCTGGTCATCTCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAgg         >  2:386518/1‑150 (MQ=255)
                                                                                                                                              ggATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAgg         <  1:1258568/150‑1 (MQ=255)
                                                                                                                                              ggATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTATTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAgg         <  2:354899/150‑1 (MQ=255)
                                                                                                                                                  gACTCTGGTCATCCCTTCATCCCCTGTTTAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTAc     <  2:1076266/150‑1 (MQ=255)
                                                                                                                                                  gACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTAc     <  1:758939/150‑1 (MQ=255)
                                                                                                                                                    ctctGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACAt   >  1:573395/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
ACTGTGAATATAAAATCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATT  >  NC_000913/2884428‑2884726

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