Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,525,439 G→A 6.8% intergenic (+16/+150) mrcA → / ← nudE peptidoglycan glycosyltransferase/peptidoglycan DD‑transpeptidase MrcA/ADP‑sugar diphosphatase NudE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,525,4390GA 6.8% 323.8 / 8.3 103intergenic (+16/+150)mrcA/nudEpeptidoglycan glycosyltransferase/peptidoglycan DD‑transpeptidase MrcA/ADP‑sugar diphosphatase NudE
Reads supporting (aligned to +/- strand):  ref base G (45/51);  new base A (5/2);  total (50/53)
Fisher's exact test for biased strand distribution p-value = 2.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 105 total aligned reads displayed.

GCACCGGGCAGTTAGCTAATGGTGGCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCAT  >  NC_000913/3525290‑3525588
                                                                                                                                                     |                                                                                                                                                     
gcaCCGGGCAGTTAGCTAATGGTGGCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCg                                                                                                                                                       >  2:1803291/1‑150 (MQ=255)
  aCCGGGCAGTTAGCTAATGGTGGCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGgc                                                                                                                                                     <  1:1249305/150‑1 (MQ=255)
     gggCAGTTAGCTAATGGTGGCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGcc                                                                                                                                                  >  2:687375/1‑150 (MQ=255)
       gCAGTTAGCTAATGGTGGCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCtt                                                                                                                                                <  1:787003/150‑1 (MQ=255)
            tAGCTAATGGTGGCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCCGAAGCGCCtttttaa                                                                                                                                           <  1:565763/150‑4 (MQ=255)
             aGCTAATGGTGGCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGt                                                                                                                                          >  2:1074307/1‑150 (MQ=255)
               cTAATGGTGGCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCCCAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGttt                                                                                                                                        >  1:1550426/1‑150 (MQ=255)
               cTAATGGTGGCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCCGAAGCGCCtttttaatca                                                                                                                                        <  1:410768/150‑7 (MQ=255)
                 aaTGGTGGCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGc                                                                                                                                      >  1:183899/1‑150 (MQ=255)
                 aaTGGTGGCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGc                                                                                                                                      <  2:1957464/150‑1 (MQ=255)
                   tggtggCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTg                                                                                                                                    <  2:81311/150‑1 (MQ=255)
                   tggtggCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTg                                                                                                                                    <  1:840877/150‑1 (MQ=255)
                    ggtggCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGa                                                                                                                                   >  2:1010090/1‑150 (MQ=255)
                     gtggCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGAc                                                                                                                                  <  1:1004460/150‑1 (MQ=255)
                       ggCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACaa                                                                                                                                >  1:1187871/1‑150 (MQ=255)
                       ggCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACaa                                                                                                                                >  1:1379965/1‑150 (MQ=255)
                         cAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAg                                                                                                                              >  2:1525960/1‑150 (MQ=255)
                         cAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAg                                                                                                                              <  1:1074307/150‑1 (MQ=255)
                             aGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAACTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCa                                                                                                                          <  1:751430/150‑1 (MQ=255)
                              gCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCAc                                                                                                                         <  1:153224/150‑1 (MQ=255)
                              gCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCAc                                                                                                                         >  2:438314/1‑150 (MQ=255)
                               ccGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACt                                                                                                                        <  1:1803291/150‑1 (MQ=255)
                               ccGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACt                                                                                                                        <  1:361699/150‑1 (MQ=255)
                                   taaGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGtt                                                                                                                    <  1:1080883/149‑1 (MQ=255)
                                         tATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGcc                                                                                                              <  2:445722/150‑1 (MQ=255)
                                          aTTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGGACTTGTTTATGCCg                                                                                                             >  1:1055308/1‑150 (MQ=255)
                                            ttCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGa                                                                                                           <  2:95617/150‑1 (MQ=255)
                                              cATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATa                                                                                                         >  2:167781/1‑150 (MQ=255)
                                                 cGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACgg                                                                                                      <  2:478530/150‑1 (MQ=255)
                                                   aaGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCg                                                                                                    <  2:1491766/150‑1 (MQ=255)
                                                     ggTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTg                                                                                                  >  1:1089364/1‑150 (MQ=255)
                                                       tACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGaa                                                                                                >  1:1657359/1‑150 (MQ=255)
                                                         cGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAAcg                                                                                              <  1:1459495/150‑1 (MQ=255)
                                                           cAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAAcgcg                                                                                            >  1:248415/1‑150 (MQ=255)
                                                             gCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGtt                                                                                          >  1:696398/1‑150 (MQ=255)
                                                               cGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTAt                                                                                        >  2:734880/1‑150 (MQ=255)
                                                               cGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGGGAACGCGTTAt                                                                                        <  2:1121181/150‑1 (MQ=255)
                                                                gACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATc                                                                                       >  2:962285/1‑150 (MQ=255)
                                                                  cacaACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCg                                                                                     <  1:734880/150‑1 (MQ=255)
                                                                  cacaACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCg                                                                                     >  1:817750/1‑150 (MQ=255)
                                                                   acaacaGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATAACGCGTTCACGCCGTATCCgg                                                                                    <  1:279530/150‑1 (MQ=255)
                                                                    caacaGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGc                                                                                   >  2:845097/1‑150 (MQ=255)
                                                                    caacaGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGc                                                                                   <  1:120070/150‑1 (MQ=255)
                                                                    caacaGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGc                                                                                   >  1:1297025/1‑150 (MQ=255)
                                                                     aacaGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCt                                                                                  <  1:1525960/150‑1 (MQ=255)
                                                                           cAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACaaa                                                                            <  1:687375/150‑1 (MQ=255)
                                                                                cACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGt                                                                       >  2:1556141/1‑150 (MQ=255)
                                                                                cACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGt                                                                       <  2:696398/150‑1 (MQ=255)
                                                                                 aCGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTTCACTTGTTTATGCCGGATACGTCGTGAACGCGTTATCCTGCTAACAAAATCGTg                                                                      <  2:383582/150‑1 (MQ=255)
                                                                                  cGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGa                                                                     <  2:840387/150‑1 (MQ=255)
                                                                                     ggTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTg                                                                      <  2:1679407/146‑1 (MQ=255)
                                                                                     ggTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTg                                                                      >  1:1679407/1‑146 (MQ=255)
                                                                                      gTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAAc                                                                 >  1:1910043/1‑150 (MQ=255)
                                                                                      gTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAAc                                                                 <  1:962285/150‑1 (MQ=255)
                                                                                        gggAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTc                                                               <  1:1941275/150‑1 (MQ=255)
                                                                                        gggAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTc                                                               <  1:947071/150‑1 (MQ=255)
                                                                                                ccATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGAGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAAtt                                                       >  2:982305/1‑150 (MQ=255)
                                                                                                  aTTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGc                                                     <  1:436614/150‑1 (MQ=255)
                                                                                                  aTTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGc                                                     <  2:1709333/150‑1 (MQ=255)
                                                                                                  aTTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGc                                                     >  2:1553526/1‑150 (MQ=255)
                                                                                                   ttATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCa                                                    <  2:74910/150‑1 (MQ=255)
                                                                                                      tCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGaa                                                 >  2:1457965/1‑150 (MQ=255)
                                                                                                        gATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAAc                                               >  2:1275934/1‑150 (MQ=255)
                                                                                                         aTAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAAcc                                              <  1:964003/150‑1 (MQ=255)
                                                                                                           aaTGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAAcccc                                            <  2:447205/150‑1 (MQ=255)
                                                                                                             tGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCAt                                          >  2:13044/1‑150 (MQ=255)
                                                                                                             tGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCAt                                          >  2:1210386/1‑150 (MQ=255)
                                                                                                              ggCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATa                                         >  1:428056/1‑150 (MQ=255)
                                                                                                                  aGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGcc                                     <  2:1297025/150‑1 (MQ=255)
                                                                                                                   ggCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCt                                    <  2:1657359/150‑1 (MQ=255)
                                                                                                                   ggCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCt                                    <  1:982305/150‑1 (MQ=255)
                                                                                                                    gCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTg                                   >  2:473894/1‑150 (MQ=255)
                                                                                                                     cacaGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGa                                  >  2:1415818/1‑150 (MQ=255)
                                                                                                                     cacaGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGa                                  >  2:1086245/1‑150 (MQ=255)
                                                                                                                     cacaGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGa                                  >  2:447373/1‑150 (MQ=255)
                                                                                                                                 tCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGc                      <  2:136918/150‑1 (MQ=255)
                                                                                                                                  caaactgaaaaGGCGCCGAAGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCa                     >  2:223571/8‑150 (MQ=255)
                                                                                                                                   tGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCAt                    >  1:1433793/1‑150 (MQ=255)
                                                                                                                                    aactgaaaaGGCGCCGAAGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATc                   >  1:1978748/6‑150 (MQ=255)
                                                                                                                                    aactgaaaaGGCGCCGAAGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATc                   >  2:47069/6‑150 (MQ=255)
                                                                                                                                     aTTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGACACCCCATAAGCCTGATAAACATTGTGCATCa                  >  1:977321/1‑150 (MQ=255)
                                                                                                                                      ttAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAg                 <  1:1798959/150‑1 (MQ=255)
                                                                                                                                      ttAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAg                 <  1:438314/150‑1 (MQ=255)
                                                                                                                                        gaaaaGGCGCCGAAGCGCCTTTTTAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGc               >  2:1797097/2‑150 (MQ=255)
                                                                                                                                        gaaaaGGCGCCGAAGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGc               >  1:1903437/2‑150 (MQ=255)
                                                                                                                                        aaaaaGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCAGCCGGc               >  2:778594/1‑150 (MQ=255)
                                                                                                                                         aaaaGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATtgtg                       >  2:1536157/1‑141 (MQ=255)
                                                                                                                                         aaaaGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATtgtg                       <  1:1536157/141‑1 (MQ=255)
                                                                                                                                          aaaGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCaa             >  2:487462/1‑150 (MQ=255)
                                                                                                                                            aGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAAc           <  1:1697114/150‑1 (MQ=255)
                                                                                                                                             ggCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACt          >  2:1958615/1‑150 (MQ=255)
                                                                                                                                             ggCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACt          <  1:1210386/150‑1 (MQ=255)
                                                                                                                                               cgcTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTc        >  2:1962725/1‑150 (MQ=255)
                                                                                                                                               cgcTTCGGCGCCTTTTCAGCTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATGCACAATGTTTATCAGGcttatggg        <  2:1862461/150‑8 (MQ=255)
                                                                                                                                                gcTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCa       >  2:967721/1‑150 (MQ=255)
                                                                                                                                                 cTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCAc      <  2:1785121/150‑1 (MQ=255)
                                                                                                                                                  ttCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACg     <  2:1089364/150‑1 (MQ=255)
                                                                                                                                                  ttCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACg     <  2:1187871/150‑1 (MQ=255)
                                                                                                                                                   tCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGc    <  2:1065424/150‑1 (MQ=255)
                                                                                                                                                     ggCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCAt  >  1:307902/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCACCGGGCAGTTAGCTAATGGTGGCAACAGCCGCGAAGAGTATTTCATCGAAGGTACGCAGCCGACACAACAGGCAGTGCACGAGGTGGGAACGACCATTATCGATAATGGCGAGGCACAGGAATTGTTCTGATTAAAAAGGCGCTTCGGCGCCTTTTCAGTTTGCTGACAAAGTGCACTTGTTTATGCCGGATACGGCGTGAACGCGTTATCCGGCTAACAAAATCGTGAAAACTCAATAAATTGCAGAAACCCCATAAGCCTGATAAACATTGTGCATCAGGCAAACTTCACGCAT  >  NC_000913/3525290‑3525588

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: