Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,884,093 (C)6→7 coding (46/1509 nt) casA ← type I‑E CRISPR system Cascade subunit CasA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,884,0871.C94.2% 229.5 / 5.5 69coding (52/1509 nt)casAtype I‑E CRISPR system Cascade subunit CasA
Reads supporting (aligned to +/- strand):  ref base . (2/2);  new base C (41/24);  total (43/26)
Fisher's exact test for biased strand distribution p-value = 6.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCCGGGGCGATAATTTGCCCAATGCAAACCAGCAGTGCTAAAGCGGCCAGTTCCATATCGTCACGGGGCAAACTTAATCGCCACTGATCTCTACTGCAGTATAGCGATTGCAGATTTATGATTTGGACTTTCCCCCCGTTTCGCGGGCGTACAGGGATCCAGTTATCAATAAGCAAATTCATTTGTTCTCCTTCATATGCTCCGACATTTCTCCTGCATTTCTATACGTCGGCACACTTCCGTTATCGGTAAAAACAAAGAAAA  >  NC_000913/2883956‑2884220
                                                                                                                                    |                                                                                                                                     
tGCCGGGGCGATAATTTGCCCAATGCAAACCAGCAGTGCTAAAGCGGCCAGTTCCATATCGTCACGGGGCAAACTTAATCGCCACTGATCTCTACTGCAGTATAGCGATTGCAGATTTATGATTTGGACTTTcccccccg                                                                                                                               >  2:905892/1‑138 (MQ=255)
  ccGGGGCGATAATTTGCCCAATGCAAACCAGCAGTGCTAAAGCGGCCAGTTCCATATCGTCACGGGGCAAACTTAATCGCCACTGATCTCTACTGCAGTATAGCGATTGCAGATTTATGATTTGGACTTTcccccccgtt                                                                                                                             >  1:1098418/1‑136 (MQ=255)
   cGGGGCGATAATTTGCCCAATGCAAACCAGCAGTGCTAAAGCGGCCAGTTCCATATCGTCACGGGGCAAACTTAATCGCCACTGATCTCTACTGCAGTATAGCGATTGCAGATTTATGATTTGGACTTTcccccccgttt                                                                                                                            <  2:237781/140‑6 (MQ=255)
     gggCGATAATTTGCCCAATGCAAACCAGCAGTGCTAAAGCGGCCAGTTCCATATCGTCACGGGGCAAACTTAATCGCCACTGATCTCTACTGCAGTATAGCGATTGCAGATTTATGATTTGGACTTTCCCCCCCGTTTcg                                                                                                                           >  1:643936/1‑140 (MQ=255)
     gggCGATAATTTGCCCAATGCAAACCAGCAGTGCTAAAGCGGCCAGTTCCATATCGTCACGGGGCAAACTTAATCGCCACTGATCTCTACTGCAGTATAGCGATTGCAGATTTATGATTTGGACTTTCCCCCCCGTTTcg                                                                                                                           <  1:522202/140‑1 (MQ=255)
      ggCGATAATTTGCCCAATGCAAACCAGCAGTGCTAAAGCGGCCAGTTCCATATCGTCACGGGGCAAACTTAATCGCCACTGATCTCTACTGCAGTATAGCGATTGCAGATTTATGATTTGGACTTTCCCCCCCGTTTcgc                                                                                                                          >  1:459395/1‑140 (MQ=255)
        cGATAATTTGCCCAATGCAAACCAGCAGTGCTAAAGCGGCCAGTTCCATATCGTCACGGGGCAAACTTAATCGCCACTGATCTCTACTGCAGTATAGCGATTGCAGATTTATGATTTGGACTTTCCCCCCCGTTTCGCgg                                                                                                                        <  2:324887/140‑1 (MQ=255)
        cGATAATTTGCCCAATGCAAACCAGCAGTGCTAAAGCGGCCAGTTCCATATCGTCACGGGGCAAACTTAATCGCCACTGATCTCTACTGCAGTATAGCGATTGCAGATTTATGATTTGGACTTTCCCCCCCGTTTCGCgg                                                                                                                        <  1:279017/140‑1 (MQ=255)
        cGATAATTTGCCCAATGCAAACCAGCAGTGCTAAAGCGGCCAGTTCCATATCGTCACGGGGCAAACTTAATCGCCACTGATCTCTACTGCAGTATAGCGATTGCAGATTTATGATTTGGACTTTCCCCCCCGTTTCGCgg                                                                                                                        <  1:374219/140‑1 (MQ=255)
        cGATAATTTGCCCAATGCAAACCAGCAGTGCTAAAGCGGCCAGTTCCATATCGTCACGGGGCAAACTTAATCGCCACTGATCTCTACTGCAGTATAGCGATTGCAGATTTATGATTTGGACTTTCCCCCCCGTTTCGCgg                                                                                                                        <  1:1073557/140‑1 (MQ=255)
          aTAATTTGCCCAATGCAAACCAGCAGTGCTAAAGCGGCCAGTTCCATATCGTCACGGGGCAATCTTAATCGCCACTGATCTCTACTGCAGTATAGCGATTGCAGATTTATGATTTGGACTTTCCCCCCGTTTCGCGGGCg                                                                                                                     >  2:1084165/1‑140 (MQ=255)
           tAATTTGCCCAATGCAAACCAGCAGTGCTAAAGCGGCCAGTTCCATATCGTCACGGGGCAAACTTAATCGCCACTGATCTCTACTGCAGTATAGCGATTGCAGATTTATGATTTGGACTTTCCCCCCCGTTTCGCGGGCg                                                                                                                     <  1:973750/140‑1 (MQ=255)
                 gCCCAATGCAAACCAGCAGTGCTAAAGCGGCCAGTTCCATATCGTCACGGGGCAAACTTAATCGCCACTGATCTCTACTGCAGTATAGCGATTGCAGATTTATGATTTGGACTTTCCCCCCCGTTTCGCGGGCGTACAgg                                                                                                               >  2:122676/1‑140 (MQ=255)
                     aaTGCAAACCAGCAGTGCTAAAGCGGCCAGTTCCATATCGTCACGGGGCAAACTTAATCGCCACTGATCTCTACTGCAGTATAGCGATTGCAGATTTATGATTTGGACTTTCCCCCCCGTTTCGCGGGCGTACAGGGATc                                                                                                           <  1:122676/140‑1 (MQ=255)
                     aaTGCAAACCAGCAGTGCTAAAGCGGCCAGTTCCATATCGTCACGGGGCAAACTTAATCGCCACTGATCTCTACTGCAGTATAGCGATTGCAGATTTATGATTTGGACTTTCCCCCCCGTTTCGCGGGCGTACAGGGATc                                                                                                           <  1:682805/140‑1 (MQ=255)
                                 cagTGCTAAAGCGGCCAGTTCCATATCGTCACGGGGCAAACTTAATCGCCACTGATCTCTACTGCAGTATAGCGATTGCAGATTTATGATTTGGACTTTCCCCCCCGTTTCGCGGGCGTACAggg                                                                                                              <  2:549132/125‑1 (MQ=255)
                                 cagTGCTAAAGCGGCCAGTTCCATATCGTCACGGGGCAAACTTAATCGCCACTGATCTCTACTGCAGTATAGCGATTGCAGATTTATGATTTGGACTTTCCCCCCCGTTTCGCGGGCGTACAggg                                                                                                              >  1:549132/1‑125 (MQ=255)
                                 cagTGCTAAAGCGGCCAGTTCCATATCGTCACGGGGCAAACTTAATCGCCACTGATCTCTACTGCAGTATAGCGATTGCAGATTTATGATTTGGACTTTCCCCCCCGTTTCGCGGGCGTACAGGGATCCAGTTATCAATa                                                                                               >  1:1129574/1‑140 (MQ=255)
                                            cGGCCAGTTCCATATCGTCACGGGGCAAACTTAATCGCCACTGATCTCTACTGCAGTATAGCGATTGCAGATTTATGATTTGGACTTTCCCCCCCGTTTCGCGGGCGTACAGGGATCCAGTTATCAATAAGCAAATTCAt                                                                                    <  2:916446/140‑1 (MQ=255)
                                              gCCAGTTCCATATCGTCACGGGGCAAACTTAATCGCCACTGATCTCTACTGCAGTATAGCGATTGCAGATTTATGATTTGGACTTTCCCCCCCGTTTCGCGGGCGTACAGGGATCCAGTTATCAATAAGCaaa                                                                                         >  2:716651/1‑133 (MQ=255)
                                              gCCAGTTCCATATCGTCACGGGGCAAACTTAATCGCCACTGATCTCTACTGCAGTATAGCGATTGCAGATTTATGATTTGGACTTTCCCCCCCGTTTCGCGGGCGTACAGGGATCCAGTTATCAATAAGCaaa                                                                                         <  1:716651/133‑1 (MQ=255)
                                              gCCAGTTCCATATCGTCACGGGGCAAACTTAATCGCCACTGATCTCTACTGCAGTATAGCGATTGCAGATTTATGATTTGGACTTTCCCCCCCGTTTCGCGGGCGTACAGGGATCCAGTTATCAATAAGCAAATTCAttt                                                                                  >  2:258603/1‑140 (MQ=255)
                                                 ggTTCCATATCGTCACGGGGCAAACTTAATCGCCACTGATCTCTACTGCAGTATAGCGATTGCAGATTTATGATTTGGACTTTCCCCCCCGTTTCGCGGGCGTACAGGGATCCAGTTATCAATAAGCAAATTCATTTGtt                                                                               <  1:753469/139‑1 (MQ=255)
                                                     ccATATCGTCACGGGGCAAACTTAATCGCCACTGATCTCTACTGCAGTATAGCGATTGCAGATTTATGATTTGGACTTTCCCCCCCGTTTCGCGGGCGTACAGGGATCCAGTTATCAATAAGCAAATTCATTTGTTCTcc                                                                           <  2:98612/140‑1 (MQ=255)
                                                        tatCGTCACGGGGCAAACTTAATCGCCACTGATCTCTACTGCAGTATAGCGATTGCAGATTTATGATTTGGACTTTCCCCCCCGTTTCGCGGGCGTACAGGGATCCAGTTATCAATAAGCAAATTCATTTGTTCTCCTTc                                                                        >  1:719363/1‑140 (MQ=255)
                                                        tatCGTCACGGGGCAAACTTAATCGCCACTGATCTCTACTGCAGTATAGCGATTGCAGATTTATGATTTGGACTTTCCCCCCCGTTTCGCGGGCGTACAGGGATCCAGTTATCAATAAGCAAATTCATTTGTTCTCCTTc                                                                        >  1:1212143/1‑140 (MQ=255)
                                                          tCGTCACGGGGCAAACTTAATCGCCACTGATCTCTACTGCAGTATAGCGATTGCAGATTTATGATTTGGACTTCCCCCCCGTTTCGCGGGCGTACAGGGATCCAGTTATCAATAAGCAAATTCATTTGTTCTCCTTCata                                                                     <  2:396473/140‑1 (MQ=255)
                                                               aCGGGGCAAACTTAATCGCCACTGATCTCTACTGCAGTATAGCGATTGCAGATTTATGATTTGGACTTTCCCCCCCGTTTCGCGGGCGTACAGGGATCCAGTTATCAATAAGCAAATTCATTTGTTCTCCTTCATATGCt                                                                 >  2:779271/1‑140 (MQ=255)
                                                                        aCTTAATCGCCACTGATCTCTACTGCAGTATAGCGATTGCAGATTTATGATTTGGACTTTCCCCCCCGTTTCGCGGGCGTACAGGGATCCAGTTATCAATAAGCAAATTCATTTGTTCTCCTTCATATGCTCCGACAttt                                                        >  1:305073/1‑140 (MQ=255)
                                                                              tCGCCACTGATCTCTACTGCAGTATAGCGATTGCAGATTTATGATTTGGACTTTCCCCCCCGTTTCGCGGGCGTACAGGGATCCAGTTATCAATAAGCAAATTCATTTGTTCTCCTTCATATGCTCCGGCATTTCTCCTg                                                  <  2:52247/140‑1 (MQ=255)
                                                                              tCGCCACTGATCTCTACTGCAGTATAGCGATTGCAGATTTATGATTTGGACTTTCCCCCCCGTTTCGCGGGCGTACAGGGATCCAGTTATCAATAAGCAAATTCATTTGTTCTCCTTCATATGCTCCGACATTTCTCCTg                                                  <  1:79513/140‑1 (MQ=255)
                                                                              tCGCCACTGATCTCTACTGCAGTATAGCGATTGCAGATTTATGATTTGGACTTTCCCCCCCGTTTCGCGGGCGTACAGGGATCCAGTTATCAATAAGCAAATTCATTTGTTCTCCTTCATATGCTCCGACATTTCTCCTg                                                  <  2:512729/140‑1 (MQ=255)
                                                                               cGCCACTGATCTCTACTGCAGTATAGCGATTGCAGATTTATGATTTGGACTTTCCCCCCCGTTTCGCGGGCGTACAGGGATCCAGTTATCAATAAGCAAATTCATTTGTTCTCCTTCATATGCTCCGACATTTCTCCTGc                                                 <  2:1098418/140‑1 (MQ=255)
                                                                               cGCCACTGATCTCTACTGCAGTATAGCGATTGCAGATTTATGATTTGGACTTTCCCCCCCGTTTCGCGGGCGTACAGGGATCCAGTTATCAATAAGCAAATTCATTTGTTCTCCTTCATATGCTCCGACATTTCTCCTGc                                                 >  2:750143/1‑140 (MQ=255)
                                                                                 ccACTGATCTCTACTGCAGTATAGCGATTGCAGATTTATGATTTGGACTTTCCCCCCCGTTTCGCGGGCGTACAGGGATCCAGTTATCAATAAGCAAATTCATTTGTTCTCCTTCATATGCTCCGACATTTCTCCTGCAt                                               <  2:719363/140‑1 (MQ=255)
                                                                                    cTGATCTCTACTGCAGTATAGCGATTGCAGATATATGATTTGGACTTTCCCCCCCGTTTCGCGGGCGTACAGGGATCCAGTTATCAATAAGCAAATTCATTTGTTCTCCTTCATATGCTCCGACATTTCTCCTGCATTTc                                            >  2:86113/1‑140 (MQ=255)
                                                                                         ctctACTGCAGTATAGCGATTGCAGATTTATGATTTGGACTTTCCCCCCCGTTTCGCGGGCGTACAGGGATCCAGTTATCAATAAGCAAATTCATTTGTTCTCCTTCATATGCTCCGACATTTCTCCTGCATTTCTATAc                                       >  1:819664/1‑140 (MQ=255)
                                                                                            tACTGCAGTATAGCGATTGCAGATTTATGATTTGGACTTTCCCCCCCGTTTCGCGGGCGTACAGGGATCCAGTTATCAATAAGCAAATTCATTTGTTCTCCTTCATATGCTCCGACATTTCTCCTGCATTTCTATACGTc                                    >  1:1078543/1‑140 (MQ=255)
                                                                                             aCTGCAGTATAGCGATTGCAGATTTATGATTTGGACTTTCCCCCCGTTTCGCGGGCGTACAGGGATCCAGTTATCAATAAGCAAATTCATTTGTTCTCCTTCATATGCTCCGACATTTCTCCTGCATTTCTATACGTCgg                                  <  1:1084165/140‑1 (MQ=255)
                                                                                              cTGCAGTATAGCGATTGCAGATTTATGATTTGGACTTTCCCCCCCGTTTCGCGGGCGTACAGGGATCCAGTTATCAATAAGCAAATTCATTTGTTCTCCTTCATATGCTCCGACATTTCTCCTGCATTTCTATACGTCgg                                  <  2:459395/140‑1 (MQ=255)
                                                                                              cTGCAGTATAGCGATTGCAGATTTATGATTTGGACTTTCCCCCCCGTTTCGCGGGCGTACAGGGATCCAGTTATCAATAAGCAAATTCATTTGTTCTCCTTCATATGCTCCGACATTTCTCCTGCATTTCTATACGTCgg                                  <  2:1078194/140‑1 (MQ=255)
                                                                                                 cAGTATAGCGATTGCAGATTTATGATTTGGACTTTCCCCCCCGTTTCGCGGGCGTACAGGGATCCAGTTATCAATAAGCAAATTCATTTGTTCTCCTTCATATGCTCCGACATTTCTCCTGCATTTCTATACGTCGGcac                               >  2:957655/1‑140 (MQ=255)
                                                                                                 cAGTATAGCGATTGCAGATTTATGATTTGGACTTTCCCCCCCGTTTCGCGGGCGTACAGGGATCCAGTTATCAATAAGCAAATTCATTTGTTCTCCTTCATATGCTCCGACATTTCTCCTGCATTTCTATACGTCGGcac                               >  2:506012/1‑140 (MQ=255)
                                                                                                 cAGTATAGCGATTGCAGATTTATGATTTGGACTTTCCCCCCCGTTTCGCGGGCGTACAGGGATCCAGTTATCAATAAGCAAATTCATTTGTTCTCCTTCATATGCTCCGACATTTCTCCTGCATTTCTATACGTCGGca                                >  2:947935/1‑139 (MQ=255)
                                                                                                 cAGTATAGCGATTGCAGATTTATGATTTGGACTTTCCCCCCCGTTTCGCGGGCGTACAGGGATCCAGTTATCAATAAGCAAATTCATTTATTCTCCTTCATATGCTCCGACATTTCTCCTGCATTTCTATACGTCGGca                                <  1:947935/139‑1 (MQ=255)
                                                                                                   gTATAGCGATTGCAGATTTATGATTTGGACTTTCCCCCCCGTTTCGCGGGCGTACAGGGATCCAGTTATCAATAAGCAAATTCATTTGTTCTCCTTCATATGCTCCGACATTTCTCCTGCATTTCTATACGTCGGCacac                             >  2:784071/1‑140 (MQ=255)
                                                                                                   gTATAGCGATTGCAGATTTATGATTTGGACTTTCCCCCCCGTTTCGCGGGCGTACAGGGATCCAGTTATCAATAAGCAAATTCATTTGTTCTCCTTCATATGCTCCGACATTTCTCCTGCATTTCTATACGTCGGCacac                             >  1:189182/1‑140 (MQ=255)
                                                                                                   gTATAGCGATTGCAGATTTATGATTTGGACTTTCCCCCCCGTTTCGCGGGCGTACAGGGATCCAGTTATCAATAAGCAAATTCATTTGTTCTCCTTCATATGCTCCGACATTTCTCCTGCATTTCTATACGTCGGCacac                             >  2:972105/1‑140 (MQ=255)
                                                                                                    tataGCGATTGCAGATTTATGATTTGGACTTTCCCCCCCGTTTCGCGGGCGTACAGGGATCCAGTTATCAATAAGCAAATTCATTTGTTCTCCTTCATATGCTCCGACATTTCTCCTGCATTTCTATACGTCGGCACACt                            >  2:416775/1‑140 (MQ=255)
                                                                                                      taGCGATTGCAGATTTATGATTTGGACTTTCCCCCCCGTTTCGCGGGCGTACAGGGATCCAGTTATCAATAAGCAAATTCATTTGTTCTCCTTCATATGCTCCGACATTTCTCCTGCATTTCTATACGTCGGCACACTTc                          >  2:963717/1‑140 (MQ=255)
                                                                                                      taGCGATTGCAGATTTATGATTTGGACTTTCCCCCCCGTTTCGCGGGCGTACAGGGATCCAGTTATCAATAAGCAAATTCATTTGTTCTCCTTCATATGCTCCGACATTTCTCCTGCATTTCTATACGTCGGCACACTTc                          >  2:263111/1‑140 (MQ=255)
                                                                                                         cGATTGCAGATTTATGATTTGGACTTTCCCCCCCGTTTCGCGGGCGTACAGGGATCCAGCTATCAATAAGCAAATTCATTTGTTCTCCTTCATATGCTCCGACATTTCTCCTGCATTTCTATACGTCGGCACACTTCCGt                       >  2:965318/1‑140 (MQ=255)
                                                                                                           aTTGCAGATTTATGATTTGGACTTTCCCCCCCGTTTCGCGGGCGTACAGGGATCCAGTTATCAATAAGCAAATTCATTTGTTCTCCTTCATATGCTCCGACATTTCTCCTGCATTTCTATACGTCGGCacac                             <  1:908164/132‑1 (MQ=255)
                                                                                                           aTTGCAGATTTATGATTTGGACTTTCCCCCCCGTTTCGCGGGCGTACAGGGATCCAGTTATCAATAAGCAAATTCATTTGTTCTCCTTCATATGCTCCGACATTTCTCCTGCATTTCTATACGTCGGCacac                             >  2:908164/1‑132 (MQ=255)
                                                                                                           aTTGCAGATTTATGATTTGGACTTTCCCCCCCGTTTCGCGGGCGTACAGGGATCCAGTTATCAATAAGCAAATTCATTTGTTCTCCTTCATATGCTCCGACATTTCTCCTGCATTTCTATACGTCGGCACACTTCCGTTa                     >  2:282103/1‑140 (MQ=255)
                                                                                                           aTTGCAGATTTATGATTTGGACTTTCCCCCCGTTTCGCGGGCGTACAGGGATCCAGTTATCAATAAGCAAATTCATTTGTTCTCCTTCATATGCTCCGACATTTCTCCTGCATTTCTATACGTCGGCACACTTCCGTTAt                    >  1:573347/1‑140 (MQ=255)
                                                                                                            ttGCAGATTTATGATTTGGACTTTCCCCCCCGTTTCGCGGGCGTACAGGGATCCAGTTATCAATAAGCAAATTCATTTGTTCTCCTTCATATGCTCCGACATTTCTCCTGCATTTCTATACGTCGGCACACTTCCGTTAt                    >  1:360836/1‑140 (MQ=255)
                                                                                                              gCAGATTTATGATTTGGACTTTCCCCCCCGTTTCGCGGGCGTACAGGGATCCAGTTATCAATAAGCAAATTCATTTGTTCTCCTTCATATGCTCCGACATTTCTCCTGCATTTCTATACGTCGGcac                               <  2:337870/127‑1 (MQ=255)
                                                                                                              gCAGATTTATGATTTGGACTTTCCCCCCCGTTTCGCGGGCGTACAGGGATCCAGTTATCAATAAGCAAATTCATTTGTTCTCCTTCATATGCTCCGACATTTCTCCTGCATTTCTATACGTCGGcac                               >  1:337870/1‑127 (MQ=255)
                                                                                                              gCAGATTTATGATTTGGACTTTCCCCCCCGTTTCGCGGGCGTACAGGGATCCAGTTATCAATAAGCAAATTCATTTGTTCTCCTTCATATGCTCCGACATTTCTCCTGCATTTCTATACGTCGGCACACTTCCGTTATCg                  >  1:618601/1‑140 (MQ=255)
                                                                                                              gCAGATTTATGATTTGGACTTTCCCCCCCGTTTCGCGGGCGTACAGGGATCCAGTTATCAATAAGCAAATTCATTTGTTCTCCTTCATATGCTCCGACATTTCTCCTGCATTTCTATACGTCGGCACACTTCCGTTATCg                  >  1:504362/1‑140 (MQ=255)
                                                                                                              gCAGATTTATGATTTGGACTTTCCCCCCCGTTTCGCGGGCGTACAGGGATCCAGTTATCAATAAGCAAATTCATTTGTTCTCCTTCATATGCTCCGACATTTCTCCTGCATTTCTATACGTCGGCACACTTCCGTTATCg                  >  1:934628/1‑140 (MQ=255)
                                                                                                              gCAGATTTATGATTTGGACTTTCCCCCCCGTTTCGCGGGCGTACAGGGATCCAGTTATCAATAAGCAAATTCATTTGTTCTCCTTCATATGCTCCGACATTTCTCCTGCATTTCTATACGTCGGCACACTTCCGTTATCg                  >  1:491438/1‑140 (MQ=255)
                                                                                                              gCAGATTTATGATTTGGACTTTCCCCCCCGTTTCGCGGGCGTACAGGGATCCAGTTATCAATAAGCAAATTCATTTGTTCTCCTTCATATGCTCCGACATTTCTCCTGCATTTCTATACGTCGGCACACTTCCGTTATCg                  >  2:106702/1‑140 (MQ=255)
                                                                                                              gCAGATTTATGATTTGGACTTTCCCCCCCGTTTCGCGGGCGTACAGGGATCCAGTTATCAATAAGCAAATTCATTTGTTCTCCTTCATATGCTCCGACATTTCTCCTGCATTTCTATACGTCGGCACACTTCCGTTATCg                  >  1:213358/1‑140 (MQ=255)
                                                                                                                  aTTTATGATTTGGACTTTCCCCCCCGTTTCGCGGGCGTACAGGGATCCAGTTATCAATAAGCAAATTCATTTGTTCTCCTTCATATGCTCCGACATTTCTCCTGCATTTCTATACGTCGGCACACTTCCGTTATCGGTaa              >  1:1078925/1‑140 (MQ=255)
                                                                                                                    ttATGATTTGGACTTTCCCCCCCGTTTCGCGGGCGTACAGGGATCCAGTTATCAATAAGCAAATTCATTTGTTCTCCTTCATATGCTCCGACATTTCTCCTGCATTTCTATACGTCGGCACACTTCCGTTATCGGTaaaa            >  2:414251/1‑140 (MQ=255)
                                                                                                                         aTTTGGACTTTCCCCCCCGTTTCGCGGGCGTACAGGGATCCAGTTATCAATAAGCAAATTCATTTGTTCTCCTTCATATGCTCCGACATTTCTCCTGCATTTCTATACGTCGGCACACTTCCGTTATCGGTAAAAACaaa       >  2:1115194/1‑140 (MQ=255)
                                                                                                                          tttGGACTTTCCCCCCCGTTTCGCGGGCGTACAGGGATCCAGTTATCAATAAGCAAATTCATTTGTTCTCCTTCATATGCTCCGACATTTCTCCTGCATTTCTATACGTCGGCACACTTCCGTTATCGGTAAAAACAAAg      <  2:643936/140‑1 (MQ=255)
                                                                                                                            tGGACTTTCCCCCCCGTTTCGCGGGCGTACAGGGATCCAGTTATCAATAAGCAAATTCATTTGTTCTCCTTCATATGCTCCGACATTTCTCCTGCATTTCTATACGTCGGCACACTTCCGTTATCGGTAAAAACAAAGaa    >  2:496980/1‑140 (MQ=255)
                                                                                                                              gACTTTCCCCCCCGTTTCGCGGGCGTACAGGGATCCAGTTATCAATAAGCAAATTCATTTGTTCTCCTTCATATGCTCCGACATTTCTCCTGCATTTCTATACGTCGGCACACTTCCGTTATCGGTAAAAACAAAGaaaa  >  1:372774/1‑140 (MQ=255)
                                                                                                                              gACTTTCCCCCCCGTTTCGCGGGCGTACAGGGATCCAGTTATCAATAAGCAAATTCATTTGTTCTCCTTCATATGCTCCGACATTTCTCCTGCATTTCTATACGTCGGCACACTTCCGTTATCGGTAAAAACAAAGaaaa  >  2:416894/1‑140 (MQ=255)
                                                                                                                                    |                                                                                                                                     
TGCCGGGGCGATAATTTGCCCAATGCAAACCAGCAGTGCTAAAGCGGCCAGTTCCATATCGTCACGGGGCAAACTTAATCGCCACTGATCTCTACTGCAGTATAGCGATTGCAGATTTATGATTTGGACTTTCCCCCCGTTTCGCGGGCGTACAGGGATCCAGTTATCAATAAGCAAATTCATTTGTTCTCCTTCATATGCTCCGACATTTCTCCTGCATTTCTATACGTCGGCACACTTCCGTTATCGGTAAAAACAAAGAAAA  >  NC_000913/2883956‑2884220

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