Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,477,345 A→T 5.1% A25A (GCT→GCA fkpA ← peptidyl‑prolyl cis‑trans isomerase FkpA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,477,3450AT 5.1% 302.3 / 3.1 92A25A (GCT→GCAfkpApeptidyl‑prolyl cis‑trans isomerase FkpA
Reads supporting (aligned to +/- strand):  ref base A (40/47);  new base T (2/3);  total (42/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAGTTTGATGCCCAGTTTTTCTTGTTCTTTTAGAGAGTTTTCCATGTAACGACCCAGCGAGGCACCCAGTGCATAAGCTGATTTCTGATCGTCATTTTTGAACGCTGCTTTGCTGTCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGATCGGGGCATCTCACCCCAGGGTTAACTATTATCAGAAGGGTACTATAAAGCGTTGTCGAAC  >  NC_000913/3477197‑3477494
                                                                                                                                                    |                                                                                                                                                     
ccAGTTTGATGCCCAGTTTTTCTTGTTCTTTTAGAGAGTTTTCCATGTAACGCCCCAGCGAGGCACCCAGTGCATAAGCTGATTTCTGATCGTCATTTTTGAACGCTGCTTTGCTGTCAGCAGCTGTAGCAGGTTTTGCAGCTTcagcag                                                                                                                                                      >  2:1602182/1‑150 (MQ=255)
    tttGATGCCCAGTTTTTCTTGTTCTTTTAGAGAGTTTTCCATGTAACGACCCAGCTAGGAACCCAGTGCATAAGCTGATTTCTGATCGTCATTTTTGAACGCTGCTTTGCTGTCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCaaa                                                                                                                                                  <  2:660877/150‑1 (MQ=255)
      tGATGCCCAGTTTTTCTTGTTCTTTTAGAGAGTTTTCCATGTAACGACCCAGCGAGGCACCCAGTGCATAAGCTGATTTCTGATCGTCATTTTTGAACGCTGCTTTGCTGTCAGCAGCTGTAGCAGGTTTTGCTGCTGAAGCTGCaaaac                                                                                                                                                <  1:442080/150‑2 (MQ=255)
      tGATGCCCAGTTTTTCTTGTTCTTTTAGAGAGTTTTCCATGTAACGACCCAGCGAGGCACCCAGTGCATAAGCTGATTTCTGATCGTCATTTTTGAACGCTGCTTTGCTGTCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAg                                                                                                                                                <  1:1294336/150‑1 (MQ=255)
       gATGCCCAGTTTTTCTTGTTCTTTTAGAGAGTTTTCCATGTAACGACCCAGCGAGGCACCCAGTGCATAAGCTGATTTCTGATCGTCATTTTTGAACGCTGCTTTGCTGTCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCaaaa                                                                                                                                                 <  2:1567728/148‑1 (MQ=255)
       gATGCCCAGTTTTTCTTGTTCTTTTAGAGAGTTTTCCATGTAACGACCCAGCGAGGCACCCAGTGCATAAGCTGATTTCTGATCGTCATTTTTGAACGCTGCTTTGCTGTCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCaaaa                                                                                                                                                 >  1:1567728/1‑148 (MQ=255)
        ctGCCCAGCTTTTCTTGTTCTTTTAGAGAGTTTTCCATGTAACGACCCAGCGAGGCACCCAGTGCATAAGCTGATTTCTGATCGTCATTTTTGAACGCTGCTTTGCTGTCAGCAGCTGTAGCAGGTTTTGCTGCTGAAGCTGCaaaacct                                                                                                                                              <  1:984113/149‑4 (MQ=255)
         tGCCCAGTTTTTCTTGTTCTTTTAGAGAGTTTTCCATGTAACGACCCAGCGAGGCACCCAGTGCATAAGCTGATTTCTGATCGTCATTTTTGAACGCTGCTTTGCTGTCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGa                                                                                                                                             >  2:62398/1‑150 (MQ=255)
          gCCCAGTTTTTCTTGTTCTTTTAGAGAGTTTTCCATGTAACGACCCAGCGAGGCACCCAGTGCATAAGCTGATTTCTGATCGTCATTTTTGAACGCTGCTTTGCTGTCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGAt                                                                                                                                            >  2:873982/1‑150 (MQ=255)
           cccAGTTTTTCTTGTTCTTTTAGAGAGTTTTCCATGTAACGACCCAGCGAGGCACCCAGTGCATAAGCTGATTTCTGATCGTCATTTTTGAACGCTGCTTTGCTGTCAGCAGCTGTAGCAGGTTTTGCTGCTGAAGCTGCaaaacctgct                                                                                                                                           <  2:363066/150‑7 (MQ=255)
           cccAGTTTTTCTTGTTCTTTTAGAGAGTTTTCCATGTAACGACCCAGCGAGGCACCCAGTGCATAAGCTGATTTCTGATCGTCATTTTTGAACGCTGCTTTGCTGTCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGAtt                                                                                                                                           >  2:1265285/1‑150 (MQ=255)
           cccAGTTTTTCTTGTTCTTTTAGAGAGTTTTCCATGTAACGACCCAGCGAGGCACCCAGTGCATAAGCTGATTTCTGATCGTCATTTTTGAACGCTGCTTTGCTGTCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGAtt                                                                                                                                           >  1:1054855/1‑150 (MQ=255)
            ccAGTTTTTCTTGTTCTTTTAGAGAGTTTTCCATGTAACGACCCAGCGAGGCACCCAGTGCATAAGCTGATTTCTGATCGTCATTTTTGAACGCTGCTTTGCTGTCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTg                                                                                                                                          <  1:1085165/150‑1 (MQ=255)
                    tCTTGTTCTTTTAGAGAGTTTTCCATGTAACGACCCAGCGAGGCACCCAGTGCATAAGCTGATTTCTGATCGTCATTTTTGAACGCTGCTTTGCTGTCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGtgcatgc                                                                                                                                  >  1:228190/1‑150 (MQ=255)
                       tGTTCTTTTAGAGAGTTTTCCATGTAACGACCCAGCGAGGCACCCAGTGCATAAGCTGATTTCTGATCGTCATTTTTGAACGCTGCTTTGCTGTCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAgg                                                                                                                               <  1:1150443/150‑1 (MQ=255)
                          tCTTTTAGAGAGTTTTCCATGTAACGACCCAGCGAGGCACCCAGTGCATAAGCTGATTTCTGATCGTCATTTTTGAACGCTGCTTTGCTGTCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCa                                                                                                                            >  1:209984/1‑150 (MQ=255)
                           cTTTTAGAGAGTTTTCCATGTAACGACCCAGCGAGGCACCCAGTGCATAAGCTGATTTCTGATCGTCATTTTTGAACGCTGCTTTGCTGTCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCaa                                                                                                                           >  2:1563578/1‑150 (MQ=255)
                              ttAGAGAGTTTTCCATGTAACGACCCAGCGAGGCACCCAGTGCATAAGCTGATTTCTGATCGTCATTTTTGAACGCTGCTTTGCTGTCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACgg                                                                                                                        <  2:636859/150‑1 (MQ=255)
                              ttAGAGAGTTTTCCATGTAACGACCCAGCGAGGCACCCAGTGCATAAGCTGATTTCTGATCGTCATTTTTGAACGCTGCTTTGCTGTCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACgg                                                                                                                        <  2:757803/150‑1 (MQ=255)
                               tAGAGAGTTTTCCATGTAACGACCCAGCGAGGCACCCAGTGCATAAGCTGATTTCTGATCGTCATTTTTGAACGCTGCTTTGCTGTCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGc                                                                                                                       >  2:141031/1‑150 (MQ=255)
                                 gagagTTTTCCATGTAACGACCCAGCGAGGCACCCAGTGCATAAGCTGATTTCTGATCGTCATTTTTGAACGCTGCTTTGCAGTCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCa                                                                                                                     >  1:1120329/1‑150 (MQ=255)
                                 gagagTTTTCCATGTAACGACCCAGCGAGGCACCCAGTGCATAAGCCGATTTCTGATCGTCATTTTTGAACGCTGCTTTGCTGTCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCa                                                                                                                     <  2:1136059/150‑1 (MQ=255)
                                  agagTTTTCCATGTAACGACCCAGCGAGGCACCCAGTGCATAAGCTGATTTCTGATCGTCATTTTTGAACGCTGCTTTGCTGTCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCAt                                                                                                                    >  1:711703/1‑150 (MQ=255)
                                  agagTTTTCCATGTAACGACCCAGCGAGGCACCCAGTGCATAAGCTGATTTCTGATCGTCATTTTTGAACGCTGCTTTGCTGTCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCAt                                                                                                                    <  1:1047908/150‑1 (MQ=255)
                                  agagTTTTCCATGTAACGACCCAGCGAGGCACCCAGTGCATAAGCTGATTTCTGATCGTCATTTTTGAACGCTGCTTTGCTGTCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCAt                                                                                                                    <  2:1064750/150‑1 (MQ=255)
                                    agTTTTCCATGTAACGACCCAGCGAGGCACCCAGTGCATAAGCTGATTTCTGATCGTCATTTTTGAACGCTGCTTTGCTGTCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATtg                                                                                                                  <  1:603123/150‑1 (MQ=255)
                                     gTTTTCCATGTAACGACCCAGCGAGGCACCCAGTGCATAAGCTGATTTCTGATCGTCATTTTTGAACGCTGCTTTGCTGTCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATtgt                                                                                                                 <  2:1265697/150‑1 (MQ=255)
                                           cATGTAACGACCCAGCGAGGCACCCAGTGCATAAGCTGATTTCTGATCGTCATTTTTGAACGCTGCTTTGCTGTCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGc                                                                                                           >  1:756937/1‑150 (MQ=255)
                                           cATGTAACGACCCAGCGAGGCACCCAGTGCATAAGCTGATTTCTGATCGTCATTTTTGAACGCTGCTTTGCTGTCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGc                                                                                                           >  2:1041285/1‑150 (MQ=255)
                                           cATGTAACGACCCAGCGAGGCACCCAGTGCATAAGCTGATTTCTGATCGTCATTTTTGAACGCTGCTTTGCTGTCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGc                                                                                                           >  2:582223/1‑150 (MQ=255)
                                             tGTAACGACCCAGCTAGGCACCCAGTGCATAAGCTGATTTCTGATCGTCATTTTTGAACGCTGCTTTGCTGTCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGccc                                                                                                         >  2:439580/1‑149 (MQ=255)
                                               tAACGACCCAGCGAGGCACCCAGTGCATAAGCTGATTTCTGATCGTCATTTTTGAACGCTGCTTTGCTGTCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCcagc                                                                                                       >  1:1059127/1‑150 (MQ=255)
                                               tAACGACCCAGCGAGGCACCCAGTGCATAAGCTGATTTCTGATCGTCATTTTTGAACGCTGCTTTGCTGTCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCcagc                                                                                                       >  2:253103/1‑150 (MQ=255)
                                                 aCGACCCAGCGAGGCACCCAGTGCATAAGCTGATTTCTGATCGTCATTTTTGAACGCTGCTTTGCTGTCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCcagcag                                                                                                     <  1:956027/150‑1 (MQ=255)
                                                    aCCCAGCGAGGCACCCAGTGCATAAGCTGATTTCTGATCGTCATTTTTGAAAGCTGCTTTGCTGTCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGt                                                                                                  >  2:658394/1‑150 (MQ=255)
                                                         gCGAGGCACCCAGTGCATAAGCTGATTTCTGATCGTCATTTTTGAACGCTGCTTTGCTGTCAGCAGCTGTAGCAGGGTTTGCAGCTTCAGGAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACtt                                                                                             <  1:1074416/150‑1 (MQ=255)
                                                           gAGGCACCCAGTGCATAAGCTGATTTCTGATCGTCATTTTTGAACGCTGCTTTGCTGTCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTa                                                                                           <  2:924395/150‑1 (MQ=255)
                                                           gAGGCACCCAGTGCATAAGCTGATTTCTGATCGTCATTTTTGAACGCTGCTTTGCTGTCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTa                                                                                           <  2:522733/150‑1 (MQ=255)
                                                              gCACCCAGTGCATAAGCTGATTTCTGATCGTCATTTTTGAACGCTGCTTTGCTGTCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAAc                                                                                        <  1:607706/150‑1 (MQ=255)
                                                                aCCCAGTGAATAAGCTGATTTCTGATCGTCATTTTTGAACGCTGCTTTGCTGTCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAg                                                                                      >  2:787461/1‑150 (MQ=255)
                                                                   cAGTGCATAAGCTGATTTCTGATCGTCATTTTTGAACGCTGCTTTGCTGTCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGa                                                                                   >  1:901995/1‑150 (MQ=255)
                                                                      tGCATAAGCTGATTTCTGATCGTCATTTTTGAACGCTGCTTTGCTGTCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGAttt                                                                                <  2:1299341/150‑1 (MQ=255)
                                                                      tGCATAAGCTGATTTCTGATCGTCATTTTTGAACGCGGCTTTGCTGTCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGAttt                                                                                <  2:1305044/150‑1 (MQ=255)
                                                                       gCATAAGCTGATTTCTGATCGTCATTTTTGAACGCTGCTTTGCTGTCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGAttta                                                                               >  2:1432817/1‑149 (MQ=255)
                                                                         aTAAGCTGATTTCTGATCGTCATTTTTGAACGCTGCTTTGCTGTCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCAt                                                                             >  2:190060/1‑150 (MQ=255)
                                                                          tAAGCTGATTTCTGATCGTCATTTTTGAACGCTGCTTTGCTGTCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATc                                                                            <  1:1265285/150‑1 (MQ=255)
                                                                            aGCTGATTTCTGATCGTCATTTTTGAACGCTGCTTTGCTGTCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCa                                                                          >  1:1505176/1‑150 (MQ=255)
                                                                            aGCTGATTTCTGATCGTCATTTTTGAACGCTGCTTTGCTGTCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCa                                                                          >  2:284149/1‑150 (MQ=255)
                                                                            aGCTGATTTCTGATCGTCATTTTTGAACGCTGCTTTGCTGTCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCa                                                                          >  1:427880/1‑150 (MQ=255)
                                                                             gCTGATTTCTGATCGTCATTTTTGAACGCTGCTTTGCTGTCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCat                                                                         >  2:1446339/1‑150 (MQ=255)
                                                                             gCTGATTTCTGATCGTCATTTTTGAACGCTGCTTTGCTGTCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCat                                                                         >  1:120986/1‑150 (MQ=255)
                                                                                    tCTGATCGTCATTTTTGAACGCTGCTTTGCTGTCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCa                                                                  <  1:1445061/150‑1 (MQ=255)
                                                                                      tGATCGTCATTTTTGAACGCTGCTTTGCTGTCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAgg                                                                >  2:1488147/1‑150 (MQ=255)
                                                                                      tGATCGTCATTTTTGAACGCTGCTTTGCTGTCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAgg                                                                >  2:1412295/1‑150 (MQ=255)
                                                                                       gATCGTCATTTTTGAACGCTGCTTTGCTGTCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGa                                                               >  1:1457256/1‑150 (MQ=255)
                                                                                          cGTCATTTTTGAACGCTGCTTTGCTGTCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGATCg                                                            <  2:737812/150‑1 (MQ=255)
                                                                                           gTCATTTTTGAACGCTGCTTTGCTGTCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGATCgg                                                           >  2:1092452/1‑150 (MQ=255)
                                                                                           gTCATTTTTGAACGCTGCTTTGCTGTCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGATCgg                                                           <  1:1092452/150‑1 (MQ=255)
                                                                                             cATTTTTGAACGCTGCTTTGCTGTCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGATCgggg                                                         <  2:1059127/150‑1 (MQ=255)
                                                                                             cATTTTTGAACGCTGCTTTGCAGTCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGATCgggg                                                         <  2:1120329/150‑1 (MQ=255)
                                                                                                 tttGAACGCTGCTTTGCTGTCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGATCGGGGCAtc                                                     <  2:1054855/150‑1 (MQ=255)
                                                                                                  ttGAACGCTGCTTTGCTGTCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGATCGGGGCAtct                                                    >  2:1239496/1‑150 (MQ=255)
                                                                                                  ttGAACGCTGCTTTGCTGTCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGATCGGGGCAtct                                                    <  2:1041442/150‑1 (MQ=255)
                                                                                                         ctgctTTGCTGTCAGCAGCTGTAGCAGGTTTTGCAGAATCAGCAGCAAAACTGATTGGTGCATGCAGGGCAACGGCGATTATGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGCATCGGTGCATCTCACACCa                                             >  2:105689/1‑150 (MQ=255)
                                                                                                          tgctTTGCTGTCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGATCGGGGCATCTCACCCCAg                                            <  1:1446339/150‑1 (MQ=255)
                                                                                                            ctTTGCTGTCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGATCGGGGCATCTCACCCCAggg                                          <  2:618166/150‑1 (MQ=255)
                                                                                                              ttGCTGTCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGATCGGGGCATCTCACCCCAGGGtt                                        <  1:36494/150‑1 (MQ=255)
                                                                                                              ttGCTGTCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGATCGGGGCATCTCACCCCAGGGtt                                        <  1:141031/150‑1 (MQ=255)
                                                                                                               tGCTGTCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGGCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGATCGGGGCATCTCACCCCAGGGTTa                                       <  1:660390/150‑1 (MQ=255)
                                                                                                                 cTGTCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGATCGGGGCATCTCACCCCAGGGTTAAc                                     <  1:787461/150‑1 (MQ=255)
                                                                                                                 cTGTCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGATCGGGGCATCTCACCCCAGGGTTAAc                                     <  2:502799/150‑1 (MQ=255)
                                                                                                                 cTGTCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGATCGGGGCATCTCACCCCAGGGTTAAc                                     <  1:658394/150‑1 (MQ=255)
                                                                                                                  tGTCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGATCGGGGCATCTCACCCCAGGGTTAACt                                    <  2:228190/150‑1 (MQ=255)
                                                                                                                   gTCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGATCGGGGCATCTCACCCCAGGGTTAACta                                   >  1:1446519/1‑150 (MQ=255)
                                                                                                                    tCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGATCGGGGCATCTCACCCAAGGGTTAACtat                                  <  1:198110/150‑1 (MQ=255)
                                                                                                                      agcagcTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGATCGGGGCATCTCACCCCAGGGTTAACtatta                                <  2:756937/150‑1 (MQ=255)
                                                                                                                      agcagcTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGATCGGGGCATCTCACCCCAGGGTTAACtatta                                >  2:76422/1‑150 (MQ=255)
                                                                                                                       gcagcTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGATCGGGGCATCTCACCCCAGGGTTAACtattat                               >  1:1603549/1‑150 (MQ=255)
                                                                                                                        cagcTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGATCGGGGCATCTCACCCCAGGGTTAACTATTATc                              <  2:1603549/150‑1 (MQ=255)
                                                                                                                        cagcTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGATCGGGGCATCTCACCCCAGGGTTAACTATTATc                              <  1:1041285/150‑1 (MQ=255)
                                                                                                                              tAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGATCGGGGCATCTCACCCCAGGGTTAACTATTATCAGAAgg                        >  1:545732/1‑150 (MQ=255)
                                                                                                                               aatcacttttGCTGCTGAAGCTGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGATCGGGGCATCTCACCCCAGGGTTAACTATTATCAGAAggg                       >  1:456408/7‑150 (MQ=255)
                                                                                                                               aGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGATCGGGGCATCTCACCCCAGGGTTAACTATTATCAGAAggg                       <  1:1310589/150‑1 (MQ=255)
                                                                                                                               aGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGATCGGGGCATCTCACCCCAGGGTTAACTATTATCAGAAggg                       >  1:1141540/1‑150 (MQ=255)
                                                                                                                               aGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGATCGGGGCATCTCACCCCAGGGTTAACTATTATCAGAAggg                       >  1:637964/1‑150 (MQ=255)
                                                                                                                                gCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGATCGGGGCATCTCACCCCAGGGTTa                                       >  1:919840/1‑133 (MQ=255)
                                                                                                                                gCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGATCGGGGCATCTCACCCCAGGGTTa                                       <  2:919840/133‑1 (MQ=255)
                                                                                                                                  aGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGAAACGGCCATTGTGGGCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGATCGGGGCATCTCACCCCAGGGTTAACTATTATCAGAAGGGTAc                    <  2:483858/150‑1 (MQ=255)
                                                                                                                                    gTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGATCGGGGCATCTCACCCCAGGGTTAACTATTATCAGAAGGGTACta                  <  2:711703/150‑1 (MQ=255)
                                                                                                                                    cttttGCTGCTGAAGCTGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGATCGGGGCATCTCACCCCAGGGTTAACTATTATCAGAAGGGTACta                  >  1:954889/2‑150 (MQ=255)
                                                                                                                                              ttCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGATCGGGGCATCTCACCCCAGGGTTAACTATTATCAGTAGGGTACTATAAAGCGTTg        <  2:1422365/150‑1 (MQ=255)
                                                                                                                                               tCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAAGGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGATCGGGGCATCTCACCCCAGGGTTAACTATTATCAGAAGGGTACTATAAAGCGTTGt       <  1:873982/150‑1 (MQ=255)
                                                                                                                                               tCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGATCGGGGCATCTCACCCCAGGGTTAACTATTATCAGAAGGGTACTATAAAGCGTTGt       <  2:901995/150‑1 (MQ=255)
                                                                                                                                                 agcagcAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGATCGGGGCATCTCACCCCAGGGTTAACTATTATCAGAAGGGTACTATAAAGCGTTGTCg     >  1:1410869/1‑150 (MQ=255)
                                                                                                                                                 agcagcAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGGCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGATCGGGGCATCTCACCCCAGGGTTAACTATTATCAGAAGGGTACTATAAAGCGTTGTCg     <  1:1602182/150‑1 (MQ=255)
                                                                                                                                                   cagcaAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGATCGGGGCATCTCACCCCAGGGTTAACTATTATCAGAAGGGTACTATAAAGCGTTGTCGaa   <  1:637704/150‑1 (MQ=255)
                                                                                                                                                    agcaAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGATCGGGGCATCTCACCCCAGGGTTAACTATTATCAGAAGGGTACTATAAAGCGTTGTCGAAc  <  1:582223/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CCAGTTTGATGCCCAGTTTTTCTTGTTCTTTTAGAGAGTTTTCCATGTAACGACCCAGCGAGGCACCCAGTGCATAAGCTGATTTCTGATCGTCATTTTTGAACGCTGCTTTGCTGTCAGCAGCTGTAGCAGGTTTTGCAGCTTCAGCAGCAAAAGTGATTGGTGCATGCAGGGCAACGGCCATTGTGGTCGCCAGCAGCGTTACTTTAAACAGTGATTTCATCCATATCTCCAGGATCGGGGCATCTCACCCCAGGGTTAACTATTATCAGAAGGGTACTATAAAGCGTTGTCGAAC  >  NC_000913/3477197‑3477494

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: