Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,766,869 C→T R318C (CGT→TGT)  rnlA → CP4‑57 prophage; RNase LS, toxin of the RnlAB toxin‑antitoxin system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,766,8690CT98.8% 296.1 / ‑4.3 85R318C (CGT→TGT) rnlACP4‑57 prophage; RNase LS, toxin of the RnlAB toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (38/46);  total (39/46)
Fisher's exact test for biased strand distribution p-value = 4.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.53e-01

TACAGCAGCCGGCAGGTTTTGGAACTTACTTTGATAAACCTGCTGCTCATTACATTCTGAAACCGCAATTTGCAGCTACTCTTAGACCGGAACAGATTAACATCATCAGCACAGCCTATACTTTTTTTAATGTGGAACGTCATTCTCTGTTCCACATGGAAACTGTGGTCGATGCCAGCCGTATGATTTCTGATATGGCCCGGTTGATGGGTAAAGCCACTAGAGCGTGGGGAATAATCAAGGACTTATATATTGTTTGAAATCACCGGAATT  >  NC_000913/2766732‑2767004
                                                                                                                                         |                                                                                                                                       
tACAGCAGCCGGCAGGTTTTGGAACTTACTTTGATAAACCTGCTGCTCATTACATTCTGAAACCGCAATTTGCAGCTACTCTTAGACCGGAACAGATTAACATCATCAGCACAGCCTATACTTTTTTTAATGTGGAATGt                                                                                                                                       >  1:864072/1‑140 (MQ=255)
tACAGCAGCCGGCAGGTTTTGGAACTTACTTTGATAAACCTGCTGCTCATTACATTCTGAAACCGCAATTTGCAGCTACTCTTAGACCGGAACAGATTAACATCATCAGCACAGCCTATACTTTTTTTAATGTGGAATGt                                                                                                                                       >  1:389048/1‑140 (MQ=255)
 aCAGCAGCCGGCAGGTTTTGGAACTTACTTTGATAAACCTGCTGCTCATTACATTCTGAAACCGCAATTTGCAGCTACTCTTAGACCGGAACAGATTAACATCATCAGCACAGCCTATACTTTTTTTAATGTGGAATGTc                                                                                                                                      <  2:905469/140‑1 (MQ=255)
 aCAGCAGCCGGCAGGTTTTGGAACTTACTTTGATAAACCTGCTGCTCATTACATTCTGAAACCGCAATTTGCAGCTACTCTTAGACCGGAACAGATTAACATCATCAGCACAGCCTATACTTTTTTTAATGTGGAATGTc                                                                                                                                      <  1:152043/140‑1 (MQ=255)
 aCAGCAGCCGGCAGGTTTTGGAACTTACTTTGATAAACCTGCTGCTCATTACATTCTGAAACCGCAATTTGCAGCTACTCTTAGACCGGAACAGATTAACATCATCAGCACAGCCTATACTTTTTTTAATGTGGAATGTc                                                                                                                                      <  2:516686/140‑1 (MQ=255)
    gcagcCGGCAGGTTTTGGAACTTACTTTGATAAACCTGCTGCTCATTACATTCTGAAACCGCAATTTGCAGCTACTCTTAGACCGGAACAGATTAACATCATCAGCACAGCCTATACTTTTTTTAATGTGGAATGTCAtt                                                                                                                                   <  1:492320/140‑1 (MQ=255)
     cagcCGGCAGGTTTTGGAACTTACTTTGATAAACCTGCTGCTCATTACATTCTGAAACCGCAATTTGCAGCTACTCTTAGACCGGAACAGATTAACATCATCAGCACAGCCTATACTTTTTTTAATGTGGAATGTCATtc                                                                                                                                  >  1:4781/1‑140 (MQ=255)
     cagcCGGCAGGTTTTGGAACTTACTTTGATAAACCTGCTGCTCATTACATTCTGAAACCGCAATTTGCAGCTACTCTTAGACCGGAACAGATTAACATCATCAGCACAGCCTATACTTTTTTTAATGTGGAATGTCATtc                                                                                                                                  >  1:115635/1‑140 (MQ=255)
     cagcCGGCAGGTTTTGGAACTTACTTTGATAAACCTGCTGCTCATTACATTCTGAAACCGCAATTTGCAGCTACTCTTAGACCGGAACAGATTAACATCATCAGCACAGCCTATACTTTTTTTAATGTGGAATGTCATtc                                                                                                                                  >  1:387296/1‑140 (MQ=255)
      agcCGGCAGGTTTTGGAACTTACTTTGATAAACCTGCTGCTCATTACATTCTGAAACCGCAATATGCGGCTACTCTTAGACCGGAACAGATTAACATCATCAGCACAGCCTATACTTTTTTTAATGTGGAATGTCATtct                                                                                                                                 <  2:79885/140‑1 (MQ=255)
       gcCGGCAGGTTTTGGAACTTACTTTGATAAACCTGCTGCTCATTACATTCTGAAACCGCAATTTGCAGCTACTCTTAGACCGGAACAGATTAACATCATCAGCACAGCCTATACTTTTTTTAATGTGGAATGTCATtctc                                                                                                                                >  1:27710/1‑140 (MQ=255)
         cGGCAGGTTTTGGAACTTACTTTGATAAACCTGCTGCTCATTACATTCTGAAACCGCAATTTGCAGCTACTCTTAGACCGGAACAGATTAACATCATCAGCACAGCCTATACTTTTTTTAATGTGGAATGTCATTCTCTg                                                                                                                              >  1:1012391/1‑140 (MQ=255)
              ggTTTTGGAACTTACTTTGATAAACCTGCTGCTCATTACATTCTGAAACCGCAATTTGCAGCTACTCTTAGACCGGAACAGATTAACATCATCAGCACAGCCTATACTTTTTTTAATGTGGAATGTCATTCTCTGTTCca                                                                                                                         <  2:882676/140‑1 (MQ=255)
                   tGGAACTTACTTTGATAAACCTGCTGCTCATTACATTCTGAAACCGCAATTTGCAGCTACTCTTAGACCGGAACAGATTAACATCATCAGCACAGCCTATACTTTTTTTAATGTGGAATGTCATTCTCTGTTCCACATgg                                                                                                                    <  2:387296/140‑1 (MQ=255)
                         ttacttTGATAAACCTGCTGCTCATTACATTCTGAAACCGCAATTTGCAGCTACTCTTAGACCGGAACAGATTAACATCATCAGCACAGCCTATACTTTTTTTAATGTGGAATGTCATTCTCTGTTCCACATGGAAACtg                                                                                                              <  2:241024/140‑1 (MQ=255)
                          tacttTGATAAACCTGCTGCTCATTACATTCTGAAACCGCAATTTGCAGCTACTCTTAGACCGGAACAGATTAACATCATCAGCACAGCCTATACTTTTTTTAATGTGGAATGTCATTCTCTGTTc                                                                                                                           <  2:1154167/126‑1 (MQ=255)
                          tacttTGATAAACCTGCTGCTCATTACATTCTGAAACCGCAATTTGCAGCTACTCTTAGACCGGAACAGATTAACATCATCAGCACAGCCTATACTTTTTTTAATGTGGAATGTCATTCTCTGTTc                                                                                                                           >  1:1154167/1‑126 (MQ=255)
                             tttGATAAACCTGCTGCTCATTACATTCTGAAACCGCAGTTTGCAGCTACTCTTAGACCGGAACAGATTAACATCATCAGCACAGCCTATACTTTTTTTAATGTGGAATGTCATTCTCTGTTCCACATGGAAACTGTGGt                                                                                                          >  2:524119/1‑140 (MQ=255)
                              ttGATAAACCTGCTGCTCATTACATTCTGAAACCGCAATTTGCAGCTACTCTTAGACCGGAACAGATTAACATCATCAGCACAGCCTATACTTTTTTTAATGTGGAATGTCATTCTCTGTTCCACATGGAAACTGTGGTc                                                                                                         <  2:948921/140‑1 (MQ=255)
                                 aTAAACCTGCTGCTCATTACATTCTGAAACCGCAATTTGCAGCTACTCTTAGACCGGAACAGATTAACATCATCAGCACAGCCTATACTTTTTTTAATGTGGAATGTCATTCTATGTTCCACATGGAAACTGTGGTCGAt                                                                                                      <  2:4781/140‑1 (MQ=255)
                                     aCCTGCTGCTCATTACATTCTGAAACCGCAATTTGCAGCTACTCTTAGACCGGAACAGATTAACATCATCAGCACAGCCTATACTTTTTTTAATGTGGAATGTCATTCTCTGTTCCACATGGAAACTGTGGTCGATGCCa                                                                                                  >  2:43702/1‑140 (MQ=255)
                                         gctgctCATTACATTCTGAAACCGCAATTTGCAGCTACTCTTAGACCGGAACAGATTAACATCATCAGCACAGCCTATACTTTTTTTAATGTGGAATGTCATTCTCTGTTCCACATGGAAACTGTGGTCGATGCCAGCCg                                                                                              >  2:149383/1‑140 (MQ=255)
                                         gctgctCATTACATTCTGAAACCGCAATTTGCAGCTACTCTTAGACCGGAACAGATTAACATCATCAGCACAGCCTATACTTTTTTTAATGTGGAATGTCATTCTCTGTTCCACATGGAAACTGTGGTCGATGCCAGCCg                                                                                              >  2:1013253/1‑140 (MQ=255)
                                         gctgctCATTACATTCTGAAACCGCAATTTGCAGCTACTCTTAGACCGGAACAGATTAACATCATCAGCACAGCCTATACTTTTTTTAATGTGGAATGTCATTCTCTGTTCCACATGGAAACTGTGGTCGATGCCAGCCg                                                                                              >  1:419073/1‑140 (MQ=255)
                                           tgctCATTACATTCTGAAACCGCAATTTGCAGCTACTCTTAGACCGGAACAGATTAACATCATCAGCACAGCCTATACTTTTTTTAATGTGGAATGTCATTCTCTGTTCCACATGGAAACTGTGGTCGATGCCAGCCGTa                                                                                            >  1:94934/1‑140 (MQ=255)
                                            gctCATTACATTCTGAAACCGCAATTTGCAGCTACTCTTAGACCGGAACAGATTAACATCATCAGCACAGCCTATACTTTTTTTAATGTGGAATGTCATTCTCTGTTCCACATGGAAACTGTGGTCGATGCCAGCCGTAt                                                                                           <  2:419073/140‑1 (MQ=255)
                                             ctCATTACATTCTGAAACCGCAATTTGCAGCTACTCTTAGACCGGAACAGATTAACATCATCAGCACAGCCTATACTTTTTTTAATGTGGAATGTCATTCTCTGTTCCACATGGAAACTGTGGTCGATGCCAGCCGTATg                                                                                          <  1:1013253/140‑1 (MQ=255)
                                               cATTACATTCTGAAACCGCAATTTGCAGCTACTCTTAGACCGGAACAGATTAACATCATCAGCACAGCCTATACTTTTTTTAATGTGGAATGTCATTCTCTGTTCCACATGGAAACTGTGGTCGATGCCAGc                                                                                                >  1:108432/1‑132 (MQ=255)
                                               cATTACATTCTGAAACCGCAATTTGCAGCTACTCTTAGACCGGAACAGATTAACATCATCAGCACAGCCTATACTTTTTTTAATGTGGAATGTCATTCTCTGTTCCACATGGAAACTGTGGTCGATGCCAGc                                                                                                <  2:108432/132‑1 (MQ=255)
                                                aTTACATTCTGAAACCGCAATTTGCAGCTACTCTTAGACCGGAACAGATTAACATCATCAGCACAGCCTATACTTTTTTTAATGTGGAATGTCATTCTCTGTTCCACATGGAAACTGTGGTCGATGCCAGCCGTa                                                                                            <  1:778597/135‑1 (MQ=255)
                                                aTTACATTCTGAAACCGCAATTTGCAGCTACTCTTAGACCGGAACAGATTAACATCATCAGCACAGCCTATACTTTTTTTAATGTGGAATGTCATTCTCTGTTCCACATGGAAACTGTGGTCGATGCCAGCCGTa                                                                                            >  2:778597/1‑135 (MQ=255)
                                                 ttACATTCTGAAACCGCAATTTGCAGCTACTCTTAGACCGGAACAGATTAACATCATCAGCACAGCCTATACTTTTTTTAATGTGGAATGTCATTCTCTGTTCCACATGGAAACTGTGGTCGATGCCAGCCGTATGAttt                                                                                      <  2:470114/140‑1 (MQ=255)
                                                       tCTGAAACCGCAATTTGCAGCTACTCTTAGACCGGAACAGATTAACATCATCAGCACAGCCTATACTTTTTTTAATGTGGAATGTCATTCTCTGTTCCACATGGAAACTGTGGTCGATGCCAGCCGTATGATTTCTGata                                                                                >  2:193137/1‑140 (MQ=255)
                                                       tCTGAAACCGCAATTTGCAGCTACTCTTAGACCGGAACAGATTAACATCATCAGCACAGCCTATACTTTTTTTAATGTGGAATGTCATTCTCTGTTCCACATGGAAACTGTGGTCGATGCCAGCCGTATGATTTCTGata                                                                                >  2:168775/1‑140 (MQ=255)
                                                         tGAAACCGCAATTTGCAGCTACTCTTAGACCGGAACAGATTAACATCATCAGCACAGCCTATACTTTTTTTAATGTGGAATGTCATTCTCTGTTCCACATGGAAACTGTGGTCGATGCCAGCCGTATGAt                                                                                        >  2:344301/1‑130 (MQ=255)
                                                         tGAAACCGCAATTTGCAGCTACTCTTAGACCGGAACAGATTAACATCATCAGCACAGCCTATACTTTTTTTAATGTGGAATGTCATTCTCTGTTCCACATGGAAACTGTGGTCGATGCCAGCCGTATGAt                                                                                        <  1:344301/130‑1 (MQ=255)
                                                         tGAAACCGCAATTTGCAGCTACTCTTAGACCGGAACAGATTAACATCATCAGCACAGCCTATACTTTTTTTAATGTGGAATGTCATTCTCTGTTCCACATGGAAACTGTGGTCGATGCCAGCCGTATGATTTCTGATATg                                                                              >  2:138238/1‑140 (MQ=255)
                                                         tGAAACCGCAATTTGCAGCTACTCTTAGACCGGAACAGATTAACATCATCAGCACAGCCTATACTTTTTTTAATGTGGAATGTCATTCTCTGTTCCACATGGAAACTGTGGTCGATGCCAGCCGTATGATTTCTGATATg                                                                              >  1:431394/1‑140 (MQ=255)
                                                          gAAACCGCAATTTGCAGCTACTCTTAGACCGGAACAGATTAACATCATCAGCACAGCCTATACTTTTTTTAATGTGGAATGTCATTCTCTGTTCCACATGGAAACTGTGGTCGATGCCAGCCGTATGATTTCTGATATgg                                                                             >  2:674503/1‑140 (MQ=255)
                                                          gAAACCGCAATTTGCAGCTACTCTTAGACCGGAACAGATTAACATCATCAGCACAGCCTATACTTTTTTTAATGTGGAATGTCATTCTCTGTTCCACATGGAAACTGTGGTCGATGCCAGCCGTATGATTTCTGATATgg                                                                             >  1:248556/1‑140 (MQ=255)
                                                          gAAACCGCAATTTGCAGCTACTCTTAGACCGGAACAGATTAACATCATCAGCACAGCCTATACTTTTTTTAATGTGGAATGTCATTCTCTGTTCCACATGGAAACTGTGGTCGATGCCAGCCGTATGATTTCTGATATgg                                                                             >  1:555210/1‑140 (MQ=255)
                                                           aaaCCGCAATTTGCAGCTACTCTTAGACCGGAACAGATTAACATCATCAGCACAGCCTATACTTTTTTTAATGTGGAATGTCATTCTCTGTTCCACATGGAAACTGTGGTCGATGCCAGCCGTATGATTTCTGATATGGc                                                                            >  2:529554/1‑140 (MQ=255)
                                                            aaCCGCAATTTGCAGCTACTCTTAGACCGGAACAGATTAACATCATCAGCACAGCCTATACTTTTTTTAATGTGGAATGTCATTCTCTGTTCCACATGGAAACTGTGGTCGATGCCAGCCGTATGATTTCTGATATGGcc                                                                           >  1:1122334/1‑140 (MQ=255)
                                                            aaCCGCAATTTGCAGCTACTCTTAGACCGGAACAGATTAACATCATCAGCACAGCCTATACTTTTTTTAATGTGGAATGTCATTCTCTGTTCCACATGGAAACTGTGGTCGATGCCAGCCGTATGATTTCTGATATGGcc                                                                           <  1:529554/140‑1 (MQ=255)
                                                            aaCCGCAATTTGCAGCTACTCTTAGACCGGAACAGATTAACATCATCAGCACAGCCTATACTTTTTTTAATGTGGAATGTCATTCTCTGTTCCACATGGAAACTGTGGACGATGCCGGCCGTATGATTTCTGATATGGcc                                                                           <  2:389048/140‑1 (MQ=255)
                                                              ccGCAATTTGCAGCTACTCTTAGACCGGAACAGATTAACATCATCAGCACAGCCTATACTTTTTTTAATGTGGAATGTCATTCTCTGTTCCACATGGAAACTGTGGTCGATGCCAGCCGTATGATTTCTGATATGGCCCg                                                                         >  1:417724/1‑140 (MQ=255)
                                                                   aTTTGCAGCTACTCTTAGACCGGAACAGATTAACATCATCAGCACAGCCTATACTTTTTTTAATGTGGAATGTCATTCTCTGTTCCACATGGAAACTGTGGTCGATGCCAGCCGTATGATTTCTGATATGGCCCGGTTGa                                                                    <  2:529996/140‑1 (MQ=255)
                                                                   aTTTGCAGCTACTCTTAGACCGGAACAGATTAACATCATCAGCACAGCCTATACTTTTTTTAATGTGGAATGTCATTCTCTGTTCCACATGGAAACTGTGGTCGATGCCAGCCGTATGATTTCTGATATGGCCCGGTTGa                                                                    <  2:315111/140‑1 (MQ=255)
                                                                        cAGCTACTCTTAGACCGGAACAGATTAACATCATCAGCACAGCCTATACTTTTTTTAATGTGGAATGTCATTCTCTGTTCCACATGGAAACTGTGGTCGATGCCAGCCGTATGATTTCTGATATGGCCCGGTTGATGGGt                                                               <  1:168775/140‑1 (MQ=255)
                                                                        cAGCTACTCTTAGACCGGAACAGATTAACATCATCAGCACAGCCTATACTTTTTTTAATGTGGAATGTCATTCTCTGTTCCACATGGAAACTGTGGTCGATGCCAGCCGTATGATTTCTGATATGGCCCGGTTGATGGGt                                                               <  2:431394/140‑1 (MQ=255)
                                                                            tACTCTTAGACCGGAACAGATTAACATCATCAGCACAGCCTATACTTTTTTTAATGTGGAATGTCATTCTCTGTTCCACATGGAAACTGTGGTCGATGCCAGCCGTATGATTTCTGATATGGCCCGGTTGATGGGTAAAg                                                           <  1:149383/140‑1 (MQ=255)
                                                                                       cGGAACAGATTAACATCATCAGCACAGCCTATACTTTTTTTAATGTGGAATGTCATTCTCTGTTCCACATGGAAACTGTGGTCGATGCCAGCCGTATGATTTCTGATATGGCCCGGTTGATGGGTAAAGCCACTAGAGCg                                                <  1:138238/140‑1 (MQ=255)
                                                                                        ggAACAGATTAACATCATCAGCACAGCCTATACTTTTTTTAATGTGGAATGTCATTCTCTGTTCCACATGGAAACTGTGGTCGATGCCAGCCGTATGATTTCTGATATGGCCCGGTTGATGGGTAAAGCCACTAGAGCGt                                               <  2:94934/140‑1 (MQ=255)
                                                                                         gAACAGATTAACATCATCAGCACAGCCTATACTTTTTTTAATGTGGAATGTCATTCTCTGTTCCACATGGAAACTGTGGTCGATGCCAGCCGTATGATTTCTGATATGGCCCGGTTGATGGGTAAAGCCACTAGAGCGTg                                              >  1:179257/1‑140 (MQ=255)
                                                                                         gAACAGATTAACATCATCAGCACAGCCTATACTTTTTTTAATGTGGAATGTCATTCTCTGTTCCACATGGAAACTGTGGTCGATGCCAGCCGTATGATTTCTGATATGGCCCGGTTGATGGGTAAAGCCACTAGAGCGTg                                              >  2:876092/1‑140 (MQ=255)
                                                                                              gATTAACATCATCAGCACAGCCTATACTTTTTTTAATGTGGAATGTCATTCTCTGTTCCACATGGAAACTGTGGTCGATGCCAGCCGTATGATTTCTGATATGGCCCGGTTGATGGGTAAAGCCACTAGAGCGTGGGGaa                                         <  2:179257/140‑1 (MQ=255)
                                                                                               aTTAACATCATCAGCACAGCCTATACTTTTTTTAATGTGGAATGTCATTCTCTGTTCCACATGGAAACTGTGGTCGATGCCAGCCGTATGATTTCTGATATGGCCCGGTTGATGGGTAAAGCCACTAGAGCGTGGGGaat                                        >  1:1025341/1‑140 (MQ=255)
                                                                                                    catcatCAGCACAGCCTATACTTTTTTTAATGTGGAATGTCATTCTCTGTTCCACATGGAAACTGTGGTCGATGCCAGCCGTATGATTTCTGATATGGCCCGGTTGATGGGTAAAGCCACTAGAGCGt                                               >  1:613218/1‑128 (MQ=255)
                                                                                                    catcatCAGCACAGCCTATACTTTTTTTAATGTGGAATGTCATTCTCTGTTCCACATGGAAACTGTGGTCGATGCCAGCCGTATGATTTCTGATATGGCCCGGTTGATGGGTAAAGCCACTAGAGCGt                                               <  2:613218/128‑1 (MQ=255)
                                                                                                    catcatCAGCACAGCCTATACTTTTTTTAATGTGGAATGTCATTCTCTGTTCCACATGGAAACTGTGGTCGATGCCAGCCGTATGATTTCTGATATGGCCCGGTTGATGGGTAAAGCCACTAGAGCGTGGGGAATAATCa                                   <  2:274995/140‑1 (MQ=255)
                                                                                                    catcatCAGCACAGCCTATACTTTTTTTAATGTGGAATGTCATTCTCTGTTCCACATGGAAACTGTGGTCGATGCCAGCCGTATGATTTCTGATATGGCCCGGTTGATGGGTAAAGCCACTAGAGCGTGGGGAATAATCa                                   <  2:1012391/140‑1 (MQ=255)
                                                                                                      tcatcaGCACAGCCTATACTTTTTTTAATGTGGAATGTCATTCTCTGTTCCACATGGAAACTGTGGTCGATGCCAGCCGTATGATTTCTGATATGGCCCGGTTGATGGGTAAAGCCACTAGAGCGTGGGGAATAATCAAg                                 <  2:661050/140‑1 (MQ=255)
                                                                                                      tcatcaGCACAGCCTATACTTTTTTTAATGTGGAATGTCATTCTCTGTTCCACATGGAAACTGTGGTCGATGCCAGCCGTATGATTTCTGATATGGCCCGGTTGATGGGTAAAGCCACTAGAGCGTGGGGAATAATCAAg                                 <  1:876092/140‑1 (MQ=255)
                                                                                                      tcatcaGCACAGCCTATACTTTTTTTAATGTGGAATGTCATTCTCTGTTCCACATGGAAACTGTGGTCGATGCCAGCCGTATGATTTCTGATATGGCCCGGTTGATGGGTAAAGCCACTAGAGCGTGGGGAATAATCAAg                                 <  2:454922/140‑1 (MQ=255)
                                                                                                      tcatcaGCACAGCCTATACTTTTTTTAATGTGGAATGTCATTCTCTGTTCCACATGGAAACTGTGGTCGATGCCAGCCGTATGATTTCTGATATGGCCCGGTTGATGGGTAAAGCCACTAGAGCGTGGGGAATAATCAAg                                 <  1:1200378/140‑1 (MQ=255)
                                                                                                      tcatcaGCACAGCCTATACTTTTTTTAATGTGGAATGTCATTCTCTGTTCCACATGGAAACTGTGGTCGATGCCAGCCGTATGATTTCTGATATGGCCCGGTTGATGGGTAAAGCCACTAGAGCGTGGGGAATAATCAAg                                 <  2:1039915/140‑1 (MQ=255)
                                                                                                      tcatcaGCACAGCCTATACTTTTTTTAATGTGGAATGTCATTCTCTGTTCCACATGGAAACTGTGGTCGATGCCAGCCGTATGATTTCTGATATGGCCCGGTTGATGGGTAAAGCCACTAGAGCGTGGGGAATAATCAAg                                 <  2:584616/140‑1 (MQ=255)
                                                                                                        atcaGCACAGCCTATACTTTTTTTAATGTGGAATGTCATTCTCTGTTCCACATGGAAACTGTGGTCGATGCCAGCCGTATGATTTCTGATATGGCCCGGTTGATGGGTAAAGCCACTAGAGCGTGGGGAATAATCAAGGa                               <  1:193137/140‑1 (MQ=255)
                                                                                                        atcaGCACAGCCTATACTTTTTTTAATGTGGAATGTCATTCTCTGTTCCACATGGAAACTGTGGTCGATGCCAGCCGTATGATTTCTGATATGGCCCGGTTGATGGGTAAAGCCACTAGAGCGTGGGGAATAATCAAGGa                               <  2:912047/140‑1 (MQ=255)
                                                                                                         tcaGCACAGCCTATACTTTTTTTAATGTGGAATGTCATTCTCTGTTCCACATGGAAACTGTGGTCGATGCCAGCCGTATGATTTCTGATATGGCCCGGTTGATGGGTAAAGCCACTAGAGCGTGGGGAATAATCAAGGAc                              >  2:739777/1‑140 (MQ=255)
                                                                                                         tcaGCACAGCCTATACTTTTTTTAATGTGGAATGTCATTCTCTGTTCCACATGGAAACTGTGGTCGATGCCAGCCGTATGATTTCTGATATGGCCCGGTTGATGGGTAAAGCCACTAGAGCGTGGGGAATAATCAAGGAc                              <  2:1025341/140‑1 (MQ=255)
                                                                                                               caGCCTATACTTTTTTTAATGTGGAATGTCATTCTCTGTTCCACATGGAAACTGTGGTCGATGCCAGTCGTATGATTTCTGATATGGCCCGGTTGATGGGTAAAGCCACTAGAGCGTGGGGAATAATCAAGGACTtatat                        >  1:845745/1‑140 (MQ=255)
                                                                                                               caGCCTATACTTTTTTTAATGTGGAATGTCATTCTCTGTTCCACATGGAAACTGTGGTCGATGCCAGCCGTATGATTTCTGATATGGCCCGGTTGATGGGTAAAGCCACTAGAGCGTGGGGAATAATCAAGGACTtatat                        <  2:1124032/140‑1 (MQ=255)
                                                                                                                 gCCTATACTTTTTTTAATGTGGAATGTCATTCTCTGTTCCACATGGAAACTGTGGTCGATGCCAGCCGTATGATTTCTGATATGGCCCGGTTGATGGGTAAAGCCACTAGAGCGTGGGGAATAATCAAGGACTTatatat                      <  2:417724/140‑1 (MQ=255)
                                                                                                                  ccTATACTTTTTTTAATGTGGAATGTCATTCTCTGTTCCACATGGAAACTGTGGTCGATGCCAGCCGTATGATTTCTGATATGGCCCGGTTGATGGGTAAAGCCACTAGAGCGTGGGGAATAATCAAg                                 <  2:548858/128‑1 (MQ=255)
                                                                                                                  ccTATACTTTTTTTAATGTGGAATGTCATTCTCTGTTCCACATGGAAACTGTGGTCGATGCCAGCCGTATGATTTCTGATATGGCCCGGTTGATGGGTAAAGCCACTAGAGCGTGGGGAATAATCAAg                                 >  1:548858/1‑128 (MQ=255)
                                                                                                                  ccTATACTTTTTTTAATGTGGAATGTCATTCTCTGTTCCACATGGAAACTGTGGTCGATGCCAGCCGTATGATTTCTGATATGGCCCGGTTGATGGGTAAAGCCACTAGAGCGTGGGGAATAATCAAGGACTTATATAtt                     >  2:259495/1‑140 (MQ=255)
                                                                                                                  ccTATACTTTTTTTAATGTGGAATGTCATTCTCTGTTCCACATGGAAACTGTGGTCGATGCCAGCCGTATGATTTCTGATATGGCCCGGTTGATGGGTAAAGCCACTAGAGCGTGGGGAATAATCAAGGACTTATATAtt                     >  2:325132/1‑140 (MQ=255)
                                                                                                                  ccTATACTTTTTTTAATGTGGAACGTCATTCTCTGTTCCACATGGAAACTGTGGTCGATGCCAGCCGTATGATTTCTGATATGGCCCGGTTGATGGGTAAAGCCACTAGAGCGTGGGGAATAATCAAGGACTTATATAtt                     >  2:264098/1‑140 (MQ=255)
                                                                                                                         tttttttAATGTGGAATGTCATTCTCTGTTCCACATGGAAACTGTGGTCGATGCCAGCCGTATGATTTCTGATATGGCCCGGTTGATGGGTAAAGCCACTAGAGCGTGGGGAATAATCAAGGACTTATATATTGTTTGaa              <  2:27710/140‑1 (MQ=255)
                                                                                                                            ttttAATGTGGAATGTCATTCTCTGTTCCACATGGAAACTGTGGTCGATGCCAGCCGTATGATTTCTGATATGGCCCGGTTGATGGGTAAAGCCACTAGAGCGTGGGGAATAATCAAGGACTtatata                       <  2:566062/128‑1 (MQ=255)
                                                                                                                            ttttAATGTGGAATGTCATTCTCTGTTCCACATGGAAACTGTGGTCGATGCCAGCCGTATGATTTCTGATATGGCCCGGTTGATGGGTAAAGCCACTAGAGCGTGGGGAATAATCAAGGACTtatata                       >  1:566062/1‑128 (MQ=255)
                                                                                                                            ttttAATGTGGAATGTCATTCTCTGTTCCACATGGAAACTGTGGTCGATGCCAGCCGTATGATTTCTGATATGGCCCGGTTGATGGGTAAAGCCACTAGAGCGTGGGGAATAATCAAGGACTTATATATTGTTTGAAATc           >  2:310849/1‑140 (MQ=255)
                                                                                                                              ttAATGTGGAATGTCATTCTCTGTTCCACATGGAAACTGTGGTCGATGCCAGCCGTATGATTTCTGATATGGCCCGGTTGATGGGTAAAGCCACTAGAGCGTGGGGAATAATCAAGGACTTATATATTGTTTGAAATCAc         <  1:325132/140‑1 (MQ=255)
                                                                                                                                     ggAATGTCATTCTCTGTTCCACATGGAAACTGTGGTCGATGCCAGCCGTATGATTTCTGATATGGCCCGGTTGATGGGTAAAGCCACTAGAGCGTGGGGAATAATCAAGGACTTATATATTGTTTGAAATCACCGGAAtt  <  2:248556/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                       
TACAGCAGCCGGCAGGTTTTGGAACTTACTTTGATAAACCTGCTGCTCATTACATTCTGAAACCGCAATTTGCAGCTACTCTTAGACCGGAACAGATTAACATCATCAGCACAGCCTATACTTTTTTTAATGTGGAACGTCATTCTCTGTTCCACATGGAAACTGTGGTCGATGCCAGCCGTATGATTTCTGATATGGCCCGGTTGATGGGTAAAGCCACTAGAGCGTGGGGAATAATCAAGGACTTATATATTGTTTGAAATCACCGGAATT  >  NC_000913/2766732‑2767004

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: