Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,884,579 T→G 5.3% R881R (AGA→CGA) ‡ 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,5790TG 5.3% 227.9 / 4.7 76R881R (AGA→CGA) ‡ygcBCascade complex anti‑viral R‑loop helicase‑annealase Cas3
Reads supporting (aligned to +/- strand):  ref base T (38/34);  new base G (2/2);  total (40/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.25e-01

TGAATATAAAATCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACAT  >  NC_000913/2884432‑2884725
                                                                                                                                                   |                                                                                                                                                  
tgAATATAAAATCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTg                                                                                                                                                   >  1:152889/1‑149 (MQ=255)
tgAATATAAAATCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTg                                                                                                                                                   >  3:19372/1‑149 (MQ=255)
      taAAATCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCAt                                                                                                                                              <  5:268140/148‑1 (MQ=255)
       aaaaTCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATcc                                                                                                                                            >  7:270784/1‑149 (MQ=255)
          aTCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCtg                                                                                                                                         <  1:88748/149‑2 (MQ=255)
            ctcATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCa                                                                                                                                       >  7:83859/1‑149 (MQ=255)
               aTACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCa                                                                                                                                       <  4:36318/146‑1 (MQ=255)
                tACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATccc                                                                                                                                   <  2:135785/149‑1 (MQ=255)
                  ccGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCt                                                                                                                                 <  4:193277/149‑1 (MQ=255)
                  ccGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATAGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCt                                                                                                                                 >  4:317663/1‑149 (MQ=255)
                     ggAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGta                                                                                                                              >  2:343235/1‑149 (MQ=255)
                       aaaTTAAAAGAAGATGTACATTGTGCACCTTCCCTAATTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGtata                                                                                                                            <  4:187744/149‑1 (MQ=255)
                        aaTTAAAAGAAGTTGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAACCCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCACGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCAGCCCTTCATCCCCTGTATAg                                                                                                                           <  6:108613/149‑1 (MQ=255)
                             aaaGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCTTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGtaat                                                                                                                      <  1:208002/149‑1 (MQ=255)
                             aaaGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGtaat                                                                                                                      <  4:111851/149‑1 (MQ=255)
                             aaaGAAGATGTACATTGTGCAACTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGtaat                                                                                                                      <  4:65565/149‑1 (MQ=255)
                                gaagaTGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGtaataa                                                                                                                    >  5:176881/1‑148 (MQ=255)
                                      gTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCACCTGTATAGGTAATAACAATACt                                                                                                             >  4:283209/1‑149 (MQ=255)
                                       tACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTg                                                                                                            <  4:19372/149‑1 (MQ=255)
                                         cATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGtt                                                                                                          <  4:127094/149‑1 (MQ=255)
                                           ttGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTAc                                                                                                        <  7:299650/149‑1 (MQ=255)
                                               gCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCt                                                                                                     <  6:176881/148‑1 (MQ=255)
                                                cACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACCAGAGTCATCCCTGCaaa                                                                                                                                     <  1:333037/115‑3 (MQ=255)
                                                cACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACCAGAGTCATCCCTGCaaa                                                                                                                                     >  2:333037/1‑113 (MQ=255)
                                                 aCCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGcc                                                                                                  >  3:247615/1‑149 (MQ=255)
                                                  ccTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCa                                                                                                 <  3:317663/149‑1 (MQ=255)
                                                        cTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACcca                                                                                           <  1:343235/149‑1 (MQ=255)
                                                            ttAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACccatcca                                                                                       <  8:281831/149‑1 (MQ=255)
                                                               aGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCatccatcc                                                                                    <  2:152889/149‑1 (MQ=255)
                                                                 tAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCACGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACATTACTGTTACCCTGCCAGACGCATCCATCCAg                                                                                  >  4:342784/1‑149 (MQ=255)
                                                                    ggATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGCCCCATCCATCCAGAtt                                                                               >  1:48486/1‑149 (MQ=255)
                                                                    ggATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGCCCCATCCATCCAGAtt                                                                               >  8:59242/1‑149 (MQ=255)
                                                                       tAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTg                                                                            >  4:340780/1‑149 (MQ=255)
                                                                        aaaCCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGt                                                                           >  4:251925/1‑149 (MQ=255)
                                                                         aaCCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGtt                                                                          <  4:130377/149‑1 (MQ=255)
                                                                          aCCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGttt                                                                         >  7:152614/1‑149 (MQ=255)
                                                                                 ttGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTcc                                                                  <  6:208128/149‑1 (MQ=255)
                                                                                 ttGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACCAGAGTCATCCCtg                                                                                                                                         <  6:218635/78‑2 (MQ=38)
                                                                                 ttGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACCAGAGTCATCCCtg                                                                                                                                         >  5:218635/1‑77 (MQ=38)
                                                                                    gTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCa                                                                                                 >  3:348045/1‑115 (MQ=255)
                                                                                    gTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCa                                                                                                 <  4:348045/115‑1 (MQ=255)
                                                                                    gTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTCCCAGACCCATCCATCCAGATTCTGtt                                                                          <  2:62248/138‑1 (MQ=255)
                                                                                    gTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTACCAGACCCATCCATCCAGATTCTGtt                                                                          >  1:62248/1‑138 (MQ=255)
                                                                                       ttattaTCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCa                                                             >  2:300138/1‑148 (MQ=255)
                                                                                             tCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTaa                                                       <  8:83859/148‑1 (MQ=255)
                                                                                             tCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATACATCCAGATACTGTTTCCCTTCCAGCCAAAGTAAc                                                      >  4:323360/1‑149 (MQ=255)
                                                                                                     gATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCtt                                               >  5:311638/1‑148 (MQ=255)
                                                                                                      aTAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCa                                             >  1:99960/1‑149 (MQ=255)
                                                                                                         aCAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATc                                           <  6:311638/148‑1 (MQ=255)
                                                                                                          cAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATcta                                         >  4:148308/1‑149 (MQ=255)
                                                                                                           aaTCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATctac                                        >  3:240713/1‑149 (MQ=255)
                                                                                                           aaTCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATctac                                        >  1:124542/1‑149 (MQ=255)
                                                                                                           aaTCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCATGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATctac                                        >  4:42835/1‑149 (MQ=255)
                                                                                                             tCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATctacta                                      >  6:215917/1‑149 (MQ=255)
                                                                                                               aTTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTactact                                    >  5:10389/1‑149 (MQ=255)
                                                                                                                ttCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACtt                                   >  6:68160/1‑149 (MQ=255)
                                                                                                                ttCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACtt                                   >  7:232009/1‑149 (MQ=255)
                                                                                                                 tCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACtt                                   >  6:315068/1‑148 (MQ=255)
                                                                                                                  cccGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACtt                                   >  7:34106/1‑147 (MQ=255)
                                                                                                                   ccGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATccc                                                                                                                                   >  3:106015/1‑50 (MQ=255)
                                                                                                                   ccGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATccc                                                                                                                                   <  4:106015/50‑1 (MQ=255)
                                                                                                                   ccGAAGTTAATTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCa                                 <  2:262343/148‑1 (MQ=255)
                                                                                                                    cGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAAccc                                                    <  4:85953/128‑1 (MQ=255)
                                                                                                                    cGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAAccc                                                    >  3:85953/1‑128 (MQ=255)
                                                                                                                     gAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATctac                                        >  4:17266/1‑139 (MQ=255)
                                                                                                                     gAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATctac                                        <  3:17266/139‑1 (MQ=255)
                                                                                                                         ttATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATCCTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAAc                          <  4:247615/149‑1 (MQ=255)
                                                                                                                            tttGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTAc                       >  4:101760/1‑149 (MQ=255)
                                                                                                                            tttGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTAc                       <  2:118444/149‑1 (MQ=255)
                                                                                                                                   tttGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCa                <  4:240713/149‑1 (MQ=255)
                                                                                                                                       cAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCAt                                            <  7:16059/117‑1 (MQ=255)
                                                                                                                                       cAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCAt                                            >  8:16059/1‑117 (MQ=255)
                                                                                                                                          ggATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAg         <  3:101760/149‑1 (MQ=255)
                                                                                                                                           gATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCGCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAgg        <  3:323360/149‑1 (MQ=255)
                                                                                                                                               aCTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACGGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCAATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTa     >  2:228358/1‑148 (MQ=255)
                                                                                                                                               aCTCTGGTCATCCCGTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACACATCCAGCCAGATTCTGGTTCCCTTCCAGCCAAAGTACCCCATCTTCATCTACTACTTCAGAAAAACTACGGTTCCAGGTGAAGGGTa     >  2:237330/1‑148 (MQ=255)
                                                                                                                                                 tctGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCa                >  5:350561/1‑135 (MQ=255)
                                                                                                                                                 tctGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCa                >  7:296586/1‑135 (MQ=255)
                                                                                                                                                 tctGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCa                >  7:291352/1‑135 (MQ=255)
                                                                                                                                                 tctGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCa                <  8:291352/135‑1 (MQ=255)
                                                                                                                                                 tctGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCa                <  8:296586/135‑1 (MQ=255)
                                                                                                                                                 tctGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCa                <  6:350561/135‑1 (MQ=255)
                                                                                                                                                 tctGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACAt  >  4:254776/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TGAATATAAAATCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACAT  >  NC_000913/2884432‑2884725

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