Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 556,018 G→A 5.2% intergenic (+22/+14) cysS → / ← ybcI cysteine‑‑tRNA ligase/conserved inner membrane protein YbcI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913556,0180GA 5.2% 245.2 / 2.1 77intergenic (+22/+14)cysS/ybcIcysteine‑‑tRNA ligase/conserved inner membrane protein YbcI
Reads supporting (aligned to +/- strand):  ref base G (34/38);  new base A (2/2);  total (36/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

TTCAACAGCGTCTGGATGCCCGTAAAGCGAAAGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCAAACGGCGCGACTTTT  >  NC_000913/555869‑556161
                                                                                                                                                     |                                                                                                                                               
ttCAACAGCGTCTGGATGCCCGTAAAGCGAAAGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAg                                                                                                                                                 <  1:377008/150‑1 (MQ=255)
  cAACAGCGTCTGGATGCCCGTAAAGCGAAAGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGtt                                                                                                                                               >  2:949483/1‑150 (MQ=255)
    aCAGCGTCTGGATGCCCGTAAAGCGAAAGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGtttt                                                                                                                                             >  1:185307/1‑150 (MQ=255)
      aGCGTCTGGATGCCCGTAAAGCGAAAGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGTAATATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTTCGCTATTTCCGGATGCGAGTTTTCg                                                                                                                                           <  2:1403219/150‑1 (MQ=255)
       gCGTCTGGATGCCCGTAAAGCGAAAGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGc                                                                                                                                          <  1:1186301/150‑1 (MQ=255)
       gCGTCTGGATGCCCGTAAAGCGAAAGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGc                                                                                                                                          <  2:870285/150‑1 (MQ=255)
        cGTCTGGATGCCCGTAAAGCGAAAGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCa                                                                                                                                         <  1:1197092/150‑1 (MQ=255)
           cTGGATGCCCGTAAAGCGAAAGACTGGGCGGCGTCGGATGCGGCGCGTGATAGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATcc                                                                                                                                      >  2:1345953/1‑150 (MQ=255)
           cTGGATGCCCGTAAAGCGAAAGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAAAACTCGCATcc                                                                                                                                      <  1:1186993/150‑1 (MQ=255)
           cTGGATGCCCGTAAAGCGAAAGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCATAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATcc                                                                                                                                      >  2:1209700/1‑150 (MQ=255)
             ggATGCCCGTAAAGCGAAAGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCgg                                                                                                                                    >  1:1364207/1‑150 (MQ=255)
                   ccGTAAAGCGAAAGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCg                                                                                                                              <  1:382036/150‑1 (MQ=255)
                         aGCGAAAGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGCTGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCCTCCGGTTATCGTCTgcg                                                                                                                        >  1:984000/1‑150 (MQ=255)
                         aGCGAAAGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTgcg                                                                                                                        <  2:185307/150‑1 (MQ=255)
                          gCGAAAGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTgcgc                                                                                                                       >  1:1238189/1‑150 (MQ=255)
                           cGAAAGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGcc                                                                                                                      >  1:792068/1‑150 (MQ=255)
                             aaaGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGccac                                                                                                                    <  1:1369035/150‑1 (MQ=255)
                              aaGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGccacc                                                                                                                   <  2:799685/150‑1 (MQ=255)
                                         gcggcggATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGc                                                                                                                                          <  2:673457/116‑1 (MQ=255)
                                         gcggcggATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGc                                                                                                                                          >  1:673457/1‑116 (MQ=255)
                                          cggcggATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCTTAAGTAATTTAGCTTTTGCCGTATGCGATTTTTCTCATCCGGTTATCGTCTGGGCCACCACTACTTTCGCa                                                                                                       <  2:1490121/150‑1 (MQ=255)
                                          cggcggATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCa                                                                                                       <  1:494637/150‑1 (MQ=255)
                                          cggcggATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCa                                                                                                       <  2:908250/150‑1 (MQ=255)
                                          cggcggATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCa                                                                                                       <  2:792068/150‑1 (MQ=255)
                                           ggcggATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCCCCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCAt                                                                                                      >  1:579879/1‑150 (MQ=255)
                                            gcggATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATc                                                                                                     <  2:640212/150‑1 (MQ=255)
                                            gcggATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATc                                                                                                     <  2:608681/150‑1 (MQ=255)
                                              ggATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAg                                                                                                   <  1:765660/150‑1 (MQ=255)
                                                aggcggcgCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTa                                                                                                 <  2:610014/148‑1 (MQ=255)
                                                   cggcgCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCa                                                                                              <  1:579387/150‑1 (MQ=255)
                                                     gcgcGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATc                                                                                            >  1:345035/1‑150 (MQ=255)
                                                        cgTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGTGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCGGTAGCATcccc                                                                                         <  2:774733/150‑1 (MQ=255)
                                                        cgTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATcccc                                                                                         <  1:618134/150‑1 (MQ=255)
                                                          tGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCgg                                                                                       <  2:1134196/150‑1 (MQ=255)
                                                           gATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCCCCACAACATTCCCATCAGTAGCATCCCCGGc                                                                                      >  2:325509/1‑150 (MQ=255)
                                                           gATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGc                                                                                      >  2:727301/1‑150 (MQ=255)
                                                               gTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAAcc                                                                                  <  2:345035/150‑1 (MQ=255)
                                                                   tAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCacac                                                                              >  2:991491/1‑150 (MQ=255)
                                                                   tAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCacac                                                                              <  1:612811/150‑1 (MQ=255)
                                                                     aCGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACAccc                                                                            >  2:1198132/1‑150 (MQ=255)
                                                                          aTGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATc                                                                       <  1:1443698/150‑1 (MQ=255)
                                                                              ggATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAAtt                                                                   >  1:691458/1‑150 (MQ=255)
                                                                              ggATCGTCCTCGAATATGGCCCGCAAGGGACAAACTGGCGACGTAAGTAATTGCGCTATTGACGGATGCAAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAAtt                                                                   >  2:399595/1‑150 (MQ=255)
                                                                                       tGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAaataa                                                          <  2:221418/150‑1 (MQ=255)
                                                                                        ggAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCa                                                                 >  1:878182/1‑142 (MQ=255)
                                                                                        ggAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCa                                                                 <  2:878182/142‑1 (MQ=255)
                                                                                            gATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCAcc                                                     >  2:1494172/1‑150 (MQ=255)
                                                                                             aTGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCACCCACACCCACATCAATTCAGAAATAATCACCt                                                    >  1:1122671/1‑150 (MQ=255)
                                                                                              tGGCCCGCAAGGGACCACCTGGCGTCATAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTg                                                   <  1:1209700/150‑1 (MQ=255)
                                                                                               ggCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGa                                                  >  2:1052879/1‑150 (MQ=255)
                                                                                               ggCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGa                                                  >  2:891199/1‑150 (MQ=255)
                                                                                                 cccGCAAGGGACCACCTGGCGTCGTTAGTAATTGCGCTAGTGCCGGATGCGAGGTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACGCACATCAATTCAGAAATAATCACCTGATg                                                <  2:1238450/150‑1 (MQ=255)
                                                                                                   cGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGcc                                              <  1:432750/150‑1 (MQ=255)
                                                                                                          gACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACgg                                       >  2:436162/1‑150 (MQ=255)
                                                                                                              accTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCgtg                                   >  1:552043/1‑150 (MQ=255)
                                                                                                                 tGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCgtggtg                                <  1:1085455/150‑1 (MQ=255)
                                                                                                                 tGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCgtggtg                                <  1:991491/150‑1 (MQ=255)
                                                                                                                 tGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCgtggtg                                <  2:80148/150‑1 (MQ=255)
                                                                                                                         taagtaaTTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACgaga                        >  1:365910/1‑150 (MQ=255)
                                                                                                                         taagtaaTTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACgaga                        <  1:1456236/150‑1 (MQ=255)
                                                                                                                          aagtaaTTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGAc                       >  2:643103/1‑150 (MQ=255)
                                                                                                                            gtaaTTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACaa                     <  1:162844/150‑1 (MQ=255)
                                                                                                                               aTTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGc                  <  2:401199/150‑1 (MQ=255)
                                                                                                                                  gcgcTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGACACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCaaa               >  1:1425294/1‑150 (MQ=255)
                                                                                                                                        taaccGGATGCGAAAACTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCAAACGgcgc         >  1:129171/4‑150 (MQ=255)
                                                                                                                                         tGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCAAACGGcgcg        <  1:1345953/150‑1 (MQ=255)
                                                                                                                                         tGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCAAACGGcgcg        >  1:1371799/1‑150 (MQ=255)
                                                                                                                                         aaccGGATGCGAAAACTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCAAACGGcgcg        <  2:129171/148‑1 (MQ=255)
                                                                                                                                          gCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCAAACGGcgcgc       >  1:1436161/1‑149 (MQ=255)
                                                                                                                                          gCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCAAACGGCGCGa       >  1:894589/1‑150 (MQ=255)
                                                                                                                                          accGGATGCGAAAACTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCAAACGGCGCGa       >  1:708927/2‑150 (MQ=255)
                                                                                                                                             ggATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCAAACGGCGCGACtt    >  2:454857/1‑150 (MQ=255)
                                                                                                                                             ggATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCAAACGGCGCGACtt    >  2:469199/1‑150 (MQ=255)
                                                                                                                                              gATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCAAACGGCGCGACttt   <  2:104462/150‑1 (MQ=255)
                                                                                                                                               aTGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCAAACGGCGCGACtttt  >  2:454148/1‑150 (MQ=255)
                                                                                                                                               aTGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGCCGTGGTGTAACGAGACAATGCAAACGGCGCGACtttt  <  2:746639/150‑1 (MQ=255)
                                                                                                                                                tGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCAAACGGCGCGACtttt  >  1:1277614/1‑149 (MQ=255)
                                                                                                                                                tGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCAAACGGCGCGACtttt  <  2:1277614/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
TTCAACAGCGTCTGGATGCCCGTAAAGCGAAAGACTGGGCGGCGGCGGATGCGGCGCGTGATCGTCTTAACGAGATGGGGATCGTGCTGGAAGATGGCCCGCAAGGGACCACCTGGCGTCGTAAGTAATTGCGCTATTGCCGGATGCGAGTTTTCGCATCCGGTTATCGTCTGCGCCACCACAACATTCCCATCAGTAGCATCCCCGGCAACCACACCCACATCAATTCAGAAATAATCACCTGATGCCCGTACGGCGTGGTGTAACGAGACAATGCAAACGGCGCGACTTTT  >  NC_000913/555869‑556161

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