Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,885,766 G→A T485I (ACC→ATC)  cas3 ← CRISPR‑associated endonuclease/helicase Cas3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,885,7660GA100.0% 283.2 / NA 79T485I (ACC→ATC) cas3CRISPR‑associated endonuclease/helicase Cas3
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (48/31);  total (48/31)

GATGAGCTAGCAGATCAAAACGTTGCGCACCATTCACACCTCGCCAGTTAATGAGTGGATATGCGGAGTTATTTTCCACTGGATCTGTATGCAGACCATAAGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGGTTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTCCCAAACCA  >  NC_000913/2885629‑2885900
                                                                                                                                         |                                                                                                                                      
gATGAGCTAGCAGATCAAAACGTTGCGCACCATTCACACCTCGCCAGTTAATGAGTGGATATGCGGAGTTATTTTCCACTGGATCTGTATGCAGACCATAAGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGAtt                                                                                                                                      >  2:518947/1‑140 (MQ=255)
     gCTAGCAGATCAAAACGTTGCGCACCATTCACACCTCGCCAGTTAATGAGTGGATATGCGGAGTTATTTTCCACTGGATCTGTATGCAGACCATAAGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGa                                                                                                                                 >  1:1010329/1‑140 (MQ=255)
     gCTAGCAGATCAAAACGTTGCGCACCATTCACACCTCGCCAGTTAATGAGTGGATATGCGGAGTTATTTTCCACTGGATCTGTATGCAGACCATAAGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGa                                                                                                                                 >  2:309623/1‑140 (MQ=255)
     gCTAGCAGATCAAAACGTTGCGCACCATTCACACCTCGCCAGTTAATGAGTGGATATGCGGAGTTATTTTCCACTGGATCTGTATGCAGACCATAAGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGa                                                                                                                                 >  2:680198/1‑140 (MQ=255)
     gCTAGCAGATCAAAACGTTGCGCACCATTCACACCTCGCCAGTTAATGAGTGGATATGCGGAGTTATTTTCCACTGGATCTGTATGCAGACCATAAGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGa                                                                                                                                 >  2:884201/1‑140 (MQ=255)
     gCTAGCAGATCAAAACGTTGCGCACCATTCACACCTCGCCAGTTAATGAGTGGATATGCGGAGTTATTTTCCACTGGATCTGTATGCAGACCATAAGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGa                                                                                                                                 <  1:298680/140‑1 (MQ=255)
        aGCAGATCAAAACGTTGCGCACCATTCACACCTCGCCAGTTAATGAGTGGATATGCGGAGTTATTTTCCACTGGATCTGTATGCAGACCATAAGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGAaag                                                                                                                              <  1:1062649/140‑1 (MQ=255)
          cAGATCAAAACGTTGCGCACCATTCACACCTCGCCAGTTAATGAGTGGATATGCGGAGTTATTTTCCACTGGATCTGTATGCAGACCATAAGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGAaagaa                                                                                                                            >  1:382804/1‑140 (MQ=255)
                    cGTTGCGCACCATTCACACCTCGCCAGTTAATGAGTGGATATGCGGAGTTATTTTCCACTGGATCTGTATGCAGACCATAAGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGAAAGAAGAATAACAct                                                                                                                  <  2:74699/140‑1 (MQ=255)
                     gTTGCGCACCATTCACACCTCGCCAGTTAATGAGTGGATATGCGGAGTTATTTTCCACTGGATCTGTATGCAGACCATAAGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGAAAGAAGAATAACActc                                                                                                                 <  1:1077381/140‑1 (MQ=255)
                       tGCGCACCATTCACACCTCGCCAGTTAATGAGTGGATATGCGGAGTTATTTTCCACTGGATCTGTATGCAGACCATAAGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGAAAGAAGAATAACActccc                                                                                                               >  2:370272/1‑140 (MQ=255)
                       tGCGCACCATTCACACCTCGCCAGTTAATGAGTGGATATGCGGAGTTATTTTCCACTGGATCTGTATGCAGACCATAAGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGAAAGAAGAATAACActccc                                                                                                               >  1:429366/1‑140 (MQ=255)
                         cgcACCATTCACACCTCGCCAGTTAATGAGTGGATATGCGGAGTTATTTTCCACTGGATCTGTATGCAGACCATAAGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGAAAGAAGAATAACActccctc                                                                                                             <  2:126675/140‑1 (MQ=255)
                         cgcACCATTCACACCTCGCCAGTTAATGAGTGGATATGCGGAGTTATTTTCCACTGGATCTGTATGCAGACCATAAGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGAAAGAAGAATAACActccctc                                                                                                             >  2:511581/1‑140 (MQ=255)
                               aTTCACACCTCGCCAGTTAATGAGTGGATATGCGGAGTTATTTTCCACTGGATCTGTATGCAGACCATAAGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGAAAGAAGAATAACACTCCCTCCCACAt                                                                                                       <  1:518947/140‑1 (MQ=255)
                                     acCTCGCCAGTTAATGAGTGGATATGCGGAGTTATTTTCCACTGGATCTGTATGCAGACCATAAGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCt                                                                                                 <  1:370272/140‑1 (MQ=255)
                                     acCTCGCCAGTTAATGAGTGGATATGCGGAGTTATTTTCCACTGGATCTGTATGCAGACCATAAGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCt                                                                                                 <  1:309623/140‑1 (MQ=255)
                                        tCGCCAGTTAATGAGTGGATATGCGGAGTTATTTTCCACTGGATCTGTATGCAGACCATAAGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAg                                                                                              >  2:334514/1‑140 (MQ=255)
                                          gCCAGTTAATGAGTGGATATGCGGAGTTATTTTCCACTGGATCTGTATGCAGACCATAAGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGcc                                                                                            >  2:1051951/1‑140 (MQ=255)
                                          gCCAGTTAATGAGTGGATATGCGGAGTTATTTTCCACTGGATCTGTATGCAGACCATAAGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGcc                                                                                            >  1:13825/1‑140 (MQ=255)
                                          gCCAGTTAATGAGTGGATATGCGGAGTTATTTTCCACTGGATCTGTATGCAGACCATAAGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGcc                                                                                            >  1:710505/1‑140 (MQ=255)
                                           ccAGTTAATGAGTGGATATGCGGAGTTATTTTCCACTGGATCTGTATGCAGACCATAAGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCt                                                                                           >  2:600043/1‑140 (MQ=255)
                                               ttAATGAGTGGATATGCGGAGTTATTTTCCACTGGATCTGTATGCAGACCATAAGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAg                                                                                       >  1:701317/1‑140 (MQ=255)
                                                   tGAGTGGATATGCGGAGTTATTTTCCACTGGATCTGTATGCAGACCATAAGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACt                                                                                   <  1:1051951/140‑1 (MQ=255)
                                                     aGTGGATATGCGGAGTTATTTTCCACTGGATCTGTATGCAGACCATAAGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGc                                                                                 >  2:91549/1‑140 (MQ=255)
                                                        ggATATGCGGAGTTATTTTCCACTGGATCTGTATGCAGACCATAAGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTc                                                                              >  2:564042/1‑140 (MQ=255)
                                                         gATATGCGGAGTTATTTTCCACTGGATCTGTATGCAGACCATAAGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGAAAGAAGAATGACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTcc                                                                             >  2:254162/1‑140 (MQ=255)
                                                             tGCGGAGTTATTTTCCACTGGATCTGTATGCAGACCATAAGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCa                                                                         <  2:1010329/140‑1 (MQ=255)
                                                              gCGGAGTTATTTTCCACTGGATCTGTATGCAGACCATAAGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGAAAGAAGAATGACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCaa                                                                        <  1:254162/140‑1 (MQ=255)
                                                               cGGAGTTATTTTCCACTGGATCTGTATGCAGACCATAAGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAg                                                                       <  2:554620/140‑1 (MQ=255)
                                                               cGGAGTTATTTTCCACTGGATCTGTATGCAGACCATAAGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAg                                                                       <  1:198419/140‑1 (MQ=255)
                                                               cGGAGTTATTTTCCACTGGATCTGTATGCAGACCATAAGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAg                                                                       <  2:429366/140‑1 (MQ=255)
                                                                 gAGTTATTTTCCACTGGATCTGTATGCAGACCATAAGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCaa                                                                        >  1:786241/1‑137 (MQ=255)
                                                                 gAGTTATTTTCCACTGGATCTGTATGCAGACCATAAGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCaa                                                                        <  2:786241/137‑1 (MQ=255)
                                                                 gAGTTATTTTCCACTGGATCTGTATGCAGACCATAAGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGcc                                                                     >  2:433863/1‑140 (MQ=255)
                                                                 gAGTTATTTTCCACTGGATCTGTATGCAGACCATAAGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGcc                                                                     <  2:710505/140‑1 (MQ=255)
                                                                 gAGTTATTTTCCACTGGATCTGTATGCAGACCATAAGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGcc                                                                     >  1:162634/1‑140 (MQ=255)
                                                                   gTTATTTTCCACTGGATCTGTATGCAGACCATAAGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGt                                                                   >  1:521243/1‑140 (MQ=255)
                                                                    ttATTTTCCACTGGATCTGTATGCAGACCATAAGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGtt                                                                  >  1:932331/1‑140 (MQ=255)
                                                                    ttATTTTCCACTGGATCTGTATGCAGACCATAAGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGtt                                                                  <  1:618110/140‑1 (MQ=255)
                                                                        tttCCACTGGATCTGTATGCAGACCATAAGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCata                                                              <  2:13825/140‑1 (MQ=255)
                                                                              cTGGATCTGTATGCAGACCATAAGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGAAAGAAGTATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGgtg                                                        >  2:957825/1‑140 (MQ=255)
                                                                              cTGGATCTGTATGCAGACCATAAGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGgtg                                                        >  1:282499/1‑140 (MQ=255)
                                                                                 gATCTGTATGCAGACCATAAGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAg                                                                       <  2:967198/122‑1 (MQ=255)
                                                                                 gATCTGTATGCAGACCATAAGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAg                                                                       >  1:967198/1‑122 (MQ=255)
                                                                                        aTGCAGACCATAAGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCAt                                              <  1:433863/140‑1 (MQ=255)
                                                                                           cAGACCATAAGTATCCAGACGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGaa                                           >  1:901489/1‑140 (MQ=255)
                                                                                            aGACCATAAGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAAc                                          >  2:175190/1‑140 (MQ=255)
                                                                                            aGACCATAAGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAAc                                          >  2:831357/1‑140 (MQ=255)
                                                                                            aGACCATAAGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAAc                                          >  1:668166/1‑140 (MQ=255)
                                                                                              aCCATAAGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGt                                                    >  1:811065/1‑128 (MQ=255)
                                                                                              aCCATAAGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGt                                                    <  2:811065/128‑1 (MQ=255)
                                                                                              aCCATAAGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATg                                             >  2:856526/1‑135 (MQ=255)
                                                                                              aCCATAAGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATg                                             <  1:856526/135‑1 (MQ=255)
                                                                                              aCCATAAGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACtt                                        >  1:449975/1‑140 (MQ=255)
                                                                                              aCCATAAGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACtt                                        >  1:857144/1‑140 (MQ=255)
                                                                                              aCCATAAGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACtt                                        >  1:694910/1‑140 (MQ=255)
                                                                                              aCCATAAGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACtt                                        >  1:367945/1‑140 (MQ=255)
                                                                                                 aTAAGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTtcat                                     >  1:180128/1‑140 (MQ=255)
                                                                                                   aaGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTtcatca                                   <  2:878532/140‑1 (MQ=255)
                                                                                                     gTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATcaac                                 <  1:735075/140‑1 (MQ=255)
                                                                                                     gTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATcaac                                 >  2:735075/1‑140 (MQ=255)
                                                                                                        tCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGTCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAAt                              >  2:636579/1‑140 (MQ=255)
                                                                                                        tCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCa                                               >  2:501797/1‑123 (MQ=255)
                                                                                                        tCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCa                                               <  1:501797/123‑1 (MQ=255)
                                                                                                        tCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAAt                              >  1:1081301/1‑140 (MQ=255)
                                                                                                        tCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAAt                              >  2:553072/1‑140 (MQ=255)
                                                                                                        tCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAAt                              >  2:518171/1‑140 (MQ=255)
                                                                                                        tCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAAt                              >  1:758835/1‑140 (MQ=255)
                                                                                                        tCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAAt                              >  2:298766/1‑140 (MQ=255)
                                                                                                            gAAGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTaaa                          <  2:694910/140‑1 (MQ=255)
                                                                                                             aaGTTTCTGTTTTTGTTTCATTGGTAGGATTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTaaaa                         >  1:284450/1‑140 (MQ=255)
                                                                                                                    tGTTTTTGTTTCATTGGTAGGATTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCg                  >  2:1076099/1‑140 (MQ=255)
                                                                                                                       ttttGTTTCATTGGTAGGATTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGAcc               <  2:758835/140‑1 (MQ=255)
                                                                                                                       ttttGTTTCATTGGTAGGATTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGAcc               <  2:857144/140‑1 (MQ=255)
                                                                                                                              tCATTGGTAGGATTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTccc        <  1:600043/140‑1 (MQ=255)
                                                                                                                                aTTGGTAGGATTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTCCCaa      >  2:154410/1‑140 (MQ=255)
                                                                                                                                  tGGTAGGATTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTCCCAAAc    <  1:680198/140‑1 (MQ=255)
                                                                                                                                    gTAGGATTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTCCCAAACCa  <  1:139516/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                      
GATGAGCTAGCAGATCAAAACGTTGCGCACCATTCACACCTCGCCAGTTAATGAGTGGATATGCGGAGTTATTTTCCACTGGATCTGTATGCAGACCATAAGTATCCAGAAGTTTCTGTTTTTGTTTCATTGGTAGGGTTGCGGAAAGAAGAATAACACTCCCTCCCACATCAGCCTGAGCCTTGAGCACTGCCTCCAGCAAGCCGTTCATATAGGTGTCGTAAGCATGAACTTCATCAACAATTAAAACACTTCGACCAATTCCCAAACCA  >  NC_000913/2885629‑2885900

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