Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,263,306 C→G 7.9% R309P (CGT→CCT) ‡ qorA ← putative quinone oxidoreductase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,263,3060CG 7.9% 280.7 / 10.2 85R309P (CGT→CCT) ‡qorAputative quinone oxidoreductase 1
Reads supporting (aligned to +/- strand):  ref base C (36/42);  new base G (2/5);  total (38/47)
Fisher's exact test for biased strand distribution p-value = 4.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCCTCTAACAGATTCATCGTGCTGTACCCTACATACAGCCGAACTATAAAAAGAAAGGGCTTCCCAGGTGGAAGCCCTATTTCTTTTATGGAATCAGCAGGCTGGAACCTTGCGTCGCCCGGCTTTCCAGAATCTCATGCGCACGCTGCGCATCCTTCAGCGGATATTTCTGCTGCTCGGCGACATCGACCTTAATCACACCGCTGGCAATCAAAGAGAACAGTTCATTACTGGCCTCGGTTAATTCCTCCCGCGTGGTGATATAGCCTTGCAGGGAAGGGCGTGTCAC  >  NC_000913/4263162‑4263451
                                                                                                                                                |                                                                                                                                                 
cgcCTCTAACAGATTCATCGTGCTGTACCCTACATACAGCCGAACTATAAAAAGAAAGGGCTTCCCAGGTGGAAGCCCTATTTCTTTTATGGAATCAGCAGGCTGGAACCTTGCGTCGCCCGGCTTTCCAGAATCTCATGCGCACGCTgc                                                                                                                                              >  1:909255/1‑150 (MQ=255)
cgcCTCTAACAGATTCATCGTGCTGTACCCTACATACAGCCGAACTATAAAAAGAAAGGGCTTCCCAGGTGGAAGCCCTATTTCTTTTATGGAATCAGCAGGCTGGAACCTTGCGTCGCCCGGCTTTCCAGAATCTCATGCGCACGCTgc                                                                                                                                              <  2:502596/150‑1 (MQ=255)
 gcCTCTAACAGATTCATCGTGCTGTACCCTACATACAGCCGAACTATAAAAAGAAAGGGCTTCCCAGGTGGAAGCCCTATTTCTTTTATGGAATCAGCAGGCTGGAACCTTGCGTCGCCCGGCTTTCCAGAATCTCATGCGCACGCTgcg                                                                                                                                             <  1:293076/150‑1 (MQ=255)
 gcCTCTAACAGATTCATCGTGCTGTACCCTACATACAGCCGAACTATAAAAAGAAAGGGCTTCCCAGGTGGAAGCCCTATTTCTTTTATGGAATCAGCAGGCTGGAACCTTGCGTCGCCCGGCTTTCCAGAATCTCATGCGCACGCTgcg                                                                                                                                             <  2:242629/150‑1 (MQ=255)
 gcCTCTAACAGATTCATCGTGCTGTACCCTACATACAGCCGAACTATAAAAAGAAAGGGCTTCCCAGGTGGAAGCCCTATTTCTTTTATGGAATCAGCAGGCTGGAACCTTGCGTCGCCCGGCTTTCCAGAATCTCATGCGCACGCTgcg                                                                                                                                             <  2:256946/150‑1 (MQ=255)
 gcCTCTAACAGATTCATCGTGCTGTACCCTACATACAGCCGAACTATAAAAAGAAAGGGCTTCCCAGGTGGAAGCCCTATTTCTTTTATGGAATCAGCAGGCTGGAACCTTGCGTCGCCCGGCTTTCCAGAATCTCATGCGCACGCTgcg                                                                                                                                             <  2:262833/150‑1 (MQ=255)
  ccTCTAACAGATTCATCGTGCTGTACCCTACATACAGCCGAACTATAAAAAGAAAGGGCTTCCCAGGTGGAAGCCCTATTTCTTTTATGGAATCAGCAGGCTGGAACCTTGCGTCGCCCGGCTTTCCAGAATCTCATGCGCACGCTgcgc                                                                                                                                            <  1:351527/150‑1 (MQ=255)
     ctAACAGATTCATCGTGCTGTACCCTACATACAGCCGAACTATAAAAAGAAAGGGCTTCCCAGGTGGAAGCCCTATTTCTTTTATGGAATCAGCAGGCTGGAACCTTGCGTCGCCCGGCTTTCCAGAATCTCATGCGCAGCGTGCGCAtg                                                                                                                                         <  1:1274796/150‑2 (MQ=255)
       aaCAGATTCATCGTGCTGTACCCTACATACAGCCGAACTATAAAAAGAAAGGGCTTCCCAGGTGGAAGCCCTATTTCTTTTATGGAATCAGCAGGCTGGAACCTTGCGTCGCCCGGCTTTCCAGAATCTCATGCGCAGCGTGCGCAtgag                                                                                                                                       <  1:1103246/150‑4 (MQ=255)
       aaCAGATTCATCGTGCTGTACCCTACATACAGCCGAACTATAAAAAGAAAGGGCTTCCCAGGTGGAAGCCCTATTTCTTTTATGGAATCAGCAGGCTGGAACCTTGCGTCGCCCGGCTTTCCAGAATCTCATGCGCAGCGTGCGCAtgag                                                                                                                                       <  2:655473/150‑4 (MQ=255)
       aaCAGATTCATCGTGCTGTACCCTACATACAGCCGAACTATAAAAAGAAAGGGCTTCCCAGGTGGAAGCCCTATTTCTTTTATGGAATCAGCAGGCTGGAACCTTGCGTCGCCCGGCTTTCCAGAATCTCATGCGCACGCTGCGCATCCt                                                                                                                                       <  2:187221/150‑1 (MQ=255)
        aCAGATTCATCGTGCTGTACCCTACATACAGCCGAACTATAAAAAGAAAGGGCTTCCCAGGTGGAAGCCCTATTTCTTTTATGGAATCAGCAGGCTGGAACCTTGCGTCGCCCGGCTTTCCAGAATCTCATGCGCACGCTGCGCATCCtt                                                                                                                                      >  1:1507886/1‑150 (MQ=255)
                     gCTGTACCCTACATACAGCCGAACTATAAAAAGAAAGGGCTTCCCAGGTGGAAGCCCTATTTCTTTTATGGAATCAGCAGGCTGGAACCTTGCGTCGCCCGGCTTTCCAGAATCTCATGCGCACGCTGCGCATCCTTCAGCGGATATTTc                                                                                                                         >  1:889309/1‑150 (MQ=255)
                       tGTACCCTACATACAGCCGAACTATAAAAAGAAAGGGCTTCCCAGGTGGAAGCCCTATTTCTTTTATGGAATCAGCAGGCTGGAACCTTGCGTCGCCCGGCTTTCCAGAATCTCATGCGCACGCTGCGCATCCTTCAGCGGATATTTctg                                                                                                                       <  1:484364/150‑1 (MQ=255)
                       tGTACCCTACATACAGCCGAACTATAAAAAGAAAGGGCTTCCCAGGTGGAAGCCCTATTTCTTTTATGGAATCAGCAGGCTGGAACCTTGCGTCGCCCGGCTTTCCAGAATCTCATGCGCACGCTGCGCATCCTTCAGCGGATATTTctg                                                                                                                       <  2:995853/150‑1 (MQ=255)
                        ttACCCTACATACAGCCTAACTATAAAAAGAAAGGGCTTCCCAGGTGGAAGCCCTATTTCTTTTATGGAATCAGCAGGCTGGAACCTTGCGTCGCCCGGCTTTCCAGAATCTCATGCGCACGCTGCGCATCCTTCAGCGGATATTTctgc                                                                                                                      <  1:561096/149‑1 (MQ=255)
                        gTACCCTACATACAGCCGAACTATAAAAAGAAAGGGCTTCCCAGGTGGAAGCCCTATTTCTTTTATGGAATCAGCAGGCTGGAACCTTGCGTCGCCCGGCTTTCCAGAATCTCATGCGCACGCTGCGCATCCTTCAGCGGATATTTctgc                                                                                                                      >  2:1275218/1‑150 (MQ=255)
                        gTACCCTACATACAGCCGAACTATAAAAAGAAAGGGCTTCCCAGGTGGAAGCCCTATTTCTTTTATGGAATCAGCAGGCTGGAACCTTGCGTCGCCCGGCTTTCCAGAATCTCATGCGCACGCTGCGCATCCTTCAGCGGATATTTctgc                                                                                                                      <  1:1034577/150‑1 (MQ=255)
                        gTACCCTACATACAGCCGAACTATAAAAAGAAAGGGCTTCCCAGGTGGAAGCCCTATTTATTTTATGGAATCAGCAGGCTGGAACCTTGCGTCGCCCGGCTTTCCAGAATCTCATGCGCACGCTGCGCATCCTTCAGCGGATATTTctgc                                                                                                                      >  2:428779/1‑150 (MQ=255)
                           cccTACATACAGCCGAACTATAAAAAGAAAGGGCTTCCCAGGTGGAAGCCCTATTTCTTTTATGGAATCAGCAGGCTGGAACCTTGCGTCGCCCGGCTTTCCAGAATCTCATGCGCACGCTGCGCATCCTTCAGCGGATATTTCtgctgc                                                                                                                   >  1:977293/1‑150 (MQ=255)
                              tacatacaGCCGAACTATAAAAAGAAAGGGCTTCCCAGGTGGAAGCCCTATTTCTTTTATGGAATCAGCAGGCTGGAACCTTGCGTCGCCCGGCTTTCCAGAATCTCATGCGCACGCTGCGCATCCTTCAGCGGATATTTCTGCTGCTCg                                                                                                                >  2:68987/1‑150 (MQ=255)
                              tacatacaGCCGAACTATAAAAAGAAAGGGCTTCCCAGGTGGAAGCCCTATTTCTTTTATGGAATCAGCAGGCTGGAACCTTGCGTCGCCCGGCTTTCCAGAATCTCATGCGCACGCTGCGCATCCTTCAGCGGATATTTCTGCTGCTCg                                                                                                                >  1:1303524/1‑150 (MQ=255)
                                 atacaGCCGAACTATAAAAAGAAAGGGCTTCCCAGGTGGAAGCCCTATTTCTTTTATGGAATCAGCAGGCTGGAACCTTGCGTCGCCCGGCTTTCCAGAATCTCATGCGCACGCTGCGCATCCTTCAGCGGATATTTCTGCTGCTCGGCg                                                                                                             >  1:1551058/1‑150 (MQ=255)
                                 atacaGCCGAACTATAAAAAGAAAGGGCTTCCCAGGTGGAAGCCCTATTTCTTTTATGGAATCAGCAGGCTGGAACCTTGCGTCGCCCGGCTTTCCAGAATCTCATGCGCACGCTGCGCATCCTTCAGCGGATATTTCTGCTGCTCGGCg                                                                                                             <  1:946915/150‑1 (MQ=255)
                                 atacaGCCGAACTATAAAAAGAAAGGGCTTCCCAGGTGGAAGCCCTATTTCTTTTATGGAATCAGCAGGCTGGAACCTTGCGTCGCCCGGCTTTCCAGAATCTCATGCGCACGCTGCGCATCCTTCAGCAGAAATATccgctgaaggatg                                                                                                             <  1:1596465/150‑8 (MQ=255)
                                 atacaGCCGAACTATAAAAAGAAAGGGCTTCCCAGGTGGAAGCCCTATTTCTTTTATGGAATCAGCAGGCTGGAACCTTGCGTCGCCCGGCTTTCCAGAATCTCATGCGCACGCTGCGCATCCTTCAGCAGAAATATccgctgaaggatg                                                                                                             <  1:761855/150‑8 (MQ=255)
                                  tacaGCCGAACTATAAAAAGAAAGGGCTTCCCAGGTGGAAGCCCTATTTCTTTTATGGAATCAGCAGGCTGGAACCTTGCGTCGCCCGGCTTTCCAGAATCTCATGCGCACGCTGCGCATCCTTCAGCGGATATTTCTGCTGCTCGGCGa                                                                                                            >  1:414214/1‑150 (MQ=255)
                                    caGCCGAACTATAAAAAGAAAGGGCTTCCCAGGTGGAAGCCCTATTTCTTTTATGGAATCAGCAGGCTGGAACCTTGCGTCGCCCGGCTTTCCAGAATCTCATGCGCACGCTGCGCATCCTTCAGCGGATATTTCTGCTGCTCGGTGACa                                                                                                          >  1:893100/1‑150 (MQ=255)
                                      gCCGAACTATAAAAAGAAAGGGCTTCCCAGGTGGAAGCCCTATTTCTTTTATGGAATCAGCAGGCTGGAACCTTGCGTCGCCCGGCTTTCCAGAATCTCATGCGCACGCTGCGCATCCTTCAGCGGATATTTCTGCTGCTCGGCGACATc                                                                                                        >  1:1498654/1‑150 (MQ=255)
                                        cGAACTATAAAAAGAAAGGGCTTCCCAGGTGGAAGCCCTATTTCTTTTATGGAATCAGCAGGCTGGAACCTTGCGTCGCCCGGCTTTCCAGAATCTCATGCGCACGCTGCGCATCCTTCAGCGGATATTTCTGCTGCTCGGCGACATCGa                                                                                                      <  2:909255/150‑1 (MQ=255)
                                         gAACTATAAAAAGAAAGGGCTTCCCAGGTGGAAGCCCTATTTCTTTTATGGAATCAGCAGGCTGGAACCTTGCGTCGCCCGGCTTTCCAGAATCTCATGCGCACGCTGCGCATCCTTCAGCGGATATTTCTGCTGCTCGGCGACATCGAc                                                                                                     <  1:1289025/150‑1 (MQ=255)
                                         gAACTATAAAAAGAAAGGGCTTCCCAGGTGGAAGCCCTATTTCTTTTATGGAATCAGCAGGCTGGAACCTTGCGTCGCCCGGCTTTCCAGAATCTCATGCGCACGCTGCGCATCCTTCAGCGGATATTTCTGCTGCTCGGCGACATCGAc                                                                                                     >  1:6475/1‑150 (MQ=255)
                                           aCTATAAAAAGAAAGGGCTTCCCAGGTGGAAGCCCTATTTCTTTTATGGAATCAGCAGGCTGGAACCTTGCGTCGCCCGGCTTTCCAGAATCTCATGCGCACGCTGCGCATCCTTCAGCGGATATTTCTGCTGCTCGGCGACATCGACCt                                                                                                   >  2:51273/1‑150 (MQ=255)
                                              ataAAAAGAAAGGGCTTCCCAGGTGGAAGCCCTATTTCTTTTATGGAATCAGCAGGCTGGAACCTTGCGTCGCCCGGCTTTCCAGAATCTCATGCGCACGCTGCGCATCCTTCAGCGGATATTTCTGCTGCTCGGCGACATCGACCTTaa                                                                                                >  2:732723/1‑150 (MQ=255)
                                               taAAAAGAAAGGGCTTCCCAGGTGGAAGCCCTATTTCTTTTATGGAATCAGCAGGCTGGAACCTTGCGTCGCCCGGCTTTCCAGAATCTCATGCGCACGCTGCGCATCCTTCAGCGGATATTTCTGCTGCTCGGCGACATCGAcc                                                                                                    <  1:324551/145‑1 (MQ=255)
                                               taAAAAGAAAGGGCTTCCCAGGTGGAAGCCCTATTTCTTTTATGGAATCAGCAGGCTGGAACCTTGCGTCGCCCGGCTTTCCAGAATCTCATGCGCACGCTGCGCATCCTTCAGCGGATATTTCTGCTGCTCGGCGACATCGAcc                                                                                                    >  2:324551/1‑145 (MQ=255)
                                               taAAAAGAAAGGGCTTCCCAGGTGGAAGCCCTATTTCTTTTATGGAATCAGCAGGCTGGAACCTTGCGTCGCCCGGCTTTCCAGAATCTCATGCGCACGCTGCGCATCCTTCAGCGGATATTTCTGCTGCTCGGCGACATCGACCTTAAt                                                                                               >  1:1527925/1‑150 (MQ=255)
                                               taAAAAGAAAGGGCTTCCCAGGTGGAAGCCCTATTTCTTTTATGGAATCAGCAGGCTGGAACCTTGCGTCGCCCGGCTTTCCAGAATCTCATGCGCACGCTGCGCATCCTTCAGCGGATATTTCTGCTGCTCGGCGACATCGACCTTAAt                                                                                               >  1:151406/1‑150 (MQ=255)
                                               taAAAAGAAAGGGCTTCCCAGGTGGAAGCCCTATTTCTTTTATGGAATCAGCAGGCTGGAACCTTGCGTCGCCCGGCGTTCCAGAATCTCATGCGCACGCTGCGCATCCTTCAGCGGATATTTCTGCTGCTCGGCGACATCGACCTTAAt                                                                                               >  1:1627573/1‑150 (MQ=255)
                                                aaaaaGAAAGGGCTTCCCAGGTGGAAGCCCTATTTCTTTTATGGAATCAGCAGGCTGGAACCTTGCGTCGCCCGGCTTTCCAGAATCTCATGCGCACGCTGCGCATCCTTCAGCGGATATTTCTGCTGCTCGGCGACATCGACCTTAATc                                                                                              <  1:1468795/150‑1 (MQ=255)
                                                aaaaaGAAAGGGCTTCCCAGGTGGAAGCCCTATTTCTTTTATGGAATCAGCAGGCTGGAACCTTGCGTCGCCCGGCTTTCCAGAATCTCATGCGCACGCTGCGCATCCTTCAGCGGATATTTCTGCTGCTCGGCGACATCGACCTTAATc                                                                                              <  2:1087015/150‑1 (MQ=255)
                                                 aaaaGAAAGGGCTTCCCAGGTGGAAGCCCTATTTCTTTTATGGAATCAGCAGGCTGGAACCTTGCGTCGCCCGGCTTTCCAGAATCTCATGCGCACGCTGCGCATCCTTCAGCGGATATTTCTGCTGCTCGGCGACATCGACCTTAATca                                                                                             >  1:998151/1‑150 (MQ=255)
                                                   aagaaagGGCTTCCCAGGTGGAAGCCCTATTTCTTTTATGGAATCAGCAGGCTGGAACCTTGCGTCGCCCGGCTTTCCAGAATCTCATGCGCACGCTGCGCATCCTTCAGCGGATATTTCTGCTGCTCGGCGACATCGACCTTAATcaca                                                                                           <  2:957425/150‑1 (MQ=255)
                                                    agaaagGGCTTCCCAGGTGGAAGCCCTATTTCTTTTATGGAATCAGCAGGCTGGAACCTTGCGTCGCCCGGCTTTCCAGAATCTCATGCGCACGCTGCGCATCCTTCAGCGGATATTTCTGCTGCTCGGCGACATCGACCTTAATCacac                                                                                          <  1:1045426/150‑1 (MQ=255)
                                                         gggCTTCCCAGGTGGAAGCCCTATTTCTTTTATGGAATCAGCAGGCTGGAACCTTGCGTCGCCCGGCTTTCCAGAATCTCATGCGCACGCTGCGCATCCTTCAGCGGATATTTCTGCTGCTCGGCGACATCGACCTTAATCACACCGCTg                                                                                     <  2:414214/150‑1 (MQ=255)
                                                         gggCTTCCCAGGTGGAAGCCCTATTTCTTTTATGGAATCAGCAGGCTGGAACCTTGCGTCGCCCGGCTTTCCAGAATCTCATGCGCACGCTGCGCATCCTTCAGCGGATATTTCTGCTGCTCGGCGACATCGACCTTAATCACACCGCTg                                                                                     >  1:514415/1‑150 (MQ=255)
                                                         gggCTTCCCAGGTGGAAGCCCTATTTCTTTTATGGAATCAGCAGGCTGGAACCTTGCGTCACCCGGCTTTCCAGAATCTCATGCGCAGCGTGCGCAtg                                                                                                                                         <  2:1105418/98‑2 (MQ=255)
                                                         gggCTTCCCAGGTGGAAGCCCTATTTCTTTTATGGAATCAGCAGGCTGGAACCTTGCGTCACCCGGCTTTCCAGAATCTCATGCGCAGCGTGCGCAtg                                                                                                                                         >  1:1105418/1‑97 (MQ=255)
                                                          ggCTTCCCAGGTGGAAGCCCTATTTCTTTTATGGAATCAGCAGGCTGGAACCTTGCGTCGCCCGGCTTTCCAGAATCTCATGCGCACGCTGCGCATCCTTCAGCGGATATTTCTGCTGCTCGGCGACATCGACCTTAATCACACCGCTgg                                                                                    <  2:6475/150‑1 (MQ=255)
                                                           gCTTCCCAGGTGGAAGCCCTATTTCTTTTATGGAATCAGCAGGCTGGAACCTTGCGTCGCCCGGCTTTCCAGAATCTCATGCGCACGCTGCGCATCCTTCAGCGGATATTTCTGCTGCTCGGCGACATCGACCTTAATCACACCGCTGGc                                                                                   >  1:1606756/1‑150 (MQ=255)
                                                                ccAGGTGGAAGCCCTATTTCTTTTATGGAATCAGCAGGCTGGAACCTTGCGTCGCCCGGCTTTCCAGAATCTCATGCGCACGCTGCGCATCCTTCAGCGGATATTTCTGCTGCTCGGCGACATCGACCTTAATCACACCGCTGGCAATCa                                                                              >  2:508073/1‑150 (MQ=255)
                                                                 cAGGTGGAAGCCCTATTTCTTTTATGGAATCAGCAGGCTGGAACCTTGCGTCGCCCGGCTTTCCAGAATCTCATGCGCACGCTGCGCATCCTTCAGCGGATATTTCTGCTGCTCGGCGACATCGACCTTAATCacac                                                                                          >  2:535234/1‑137 (MQ=255)
                                                                 cAGGTGGAAGCCCTATTTCTTTTATGGAATCAGCAGGCTGGAACCTTGCGTCGCCCGGCTTTCCAGAATCTCATGCGCACGCTGCGCATCCTTCAGCGGATATTTCTGCTGCTCGGCGACATCGACCTTAATCacac                                                                                          <  1:535234/137‑1 (MQ=255)
                                                                   ggTGGAAGCCCTATTTCTTTTATGGAATCAGCAGGCTGGAACCTTGCGTCGCCCGGCTTTCCAGAATCTCATGCGCACGCTGCGCATCCTTCAGCGGATATTTCTGCTGCTCGGCGACATCGACCTTAATCACACCGCTGGCAATCAAag                                                                           <  2:111112/150‑1 (MQ=255)
                                                                   ctgggAAGCCCTATTTCTTTTATGGAATCAGCAGGCTGGAACCTTGCGTCGCCCGGCTTTCCAGAATCTCATGCGCACGCTGCGCATCCTTCAGCGGATATTTCTGCTGCTCGGCGACATCGACCTTAATCACACCGCTGGCAATCAAag                                                                           <  2:797929/147‑1 (MQ=255)
                                                                       gAAGCCCTATTTCTTTTATGGAATCAGCAGGCTGGAACCTGGCGTCGCCCGGCTTTCCAGAATCTCATGCGCACGCTGCGCATCCTTCAGCGGATATTTCTGCTGCTCGGCGACATCGACCTTAATCACACCGCTGGCAATCAAAGAGaa                                                                       <  1:135492/150‑1 (MQ=255)
                                                                            ccTATTTCTTTTATGGAATCAGCAGGCTGGAACCTTGCGTCGCCCGGCTTTCCAGAATCTCATGCGCACGCTGCGCATCCTTCAGCGGATATTTCTGCTGCTCGGCGACATCGACCTTAATCACACCGCTGGCAATCAAAGAGAACAGtt                                                                  <  2:998151/150‑1 (MQ=255)
                                                                            ccTATTTCTTTTATGGAATCAGCAGGCTGGAACCTTGCGTCGCCCGGCTTTCCAGAATCTCATGCGCACGCTGCGCATCCTTCAGCGGATATTTCTGCTGCTCGGCGACATCGACCTTAATCACACCGCTGGCAATCAAAGAGAACAGtt                                                                  <  2:262931/150‑1 (MQ=255)
                                                                               aTTTCTTTTATGGAATCAGCAGGCTGGAACCTTGCGTCGCCCGGCTTTCCAGAATCTCATGCGCACGCTGAGCATCCTTCACCGGATATTTCTGCTGCTCGGCGACATCGCCCTTAATCACACCGCTGGCAATCAAAGAGAACAGTTCAt                                                               >  2:1062744/1‑150 (MQ=255)
                                                                                   cTTTTATGGAATCAGCAGGCTGGAACCTTGCGTCGCCCGGCTTTCCAGAATCTCATGCGCACGCTGCGCATCCTTCAGCGGATATTTCTGCTGCTCGGCGACATCGACCTTAATCACACCGCTGGCAATCAAAGAGAACAGTTCATTACt                                                           >  1:1206713/1‑150 (MQ=255)
                                                                                       tATGGAATCAGCAGGCTGGAACCTTGCGTCGCCCGGCTTTCCAGAATCTCATGCGCACGCTGCGCATCCTTCAGCGGATATTTCTGCTGCTCGGCGACATCGACCTTAATCACACCGCTGGCAATCAAAGAGAACAGTTCATTACTGGcc                                                       >  2:657340/1‑150 (MQ=255)
                                                                                        aTGGAATCAGCAGGCTGGAACCTTGCGTCGCCCGGCTTTCCAGAATCTCATGCGCACGCTGCGCATCCTTCAGCGGATATTTCTGCTGCTCGGCGACATCGACCTTAATCACACCGCTGGCAATCAAAGAGAGCAGTTCATTTCTGGCCt                                                      >  1:1154885/1‑150 (MQ=255)
                                                                                                          aaCCTTGCGTCGCCCGGCTTTCCAGAATCTCATGCGCACGCTGCGCATCCTTCAGCGGATATTTCTGCTGCTCGGCGACATCGACCTTAATCACACCGCTGGCAATCAAAGAGAACAGTTCATTACTGGCCTCGGTTAATTCCTCCcgcg                                    <  2:1206713/150‑1 (MQ=255)
                                                                                                          aaCCTTGCGTCGCCCGGCTTTCCAGAATCTCATGCGCACGCTGCGCATCCTTCAGCGGATATTTCTGCTGCTCGGCGACATCGACCTTAATCACACCGCTGGCAATCAAAGAGAACAGTTCATTACTGGCCTCGGTTAATTCCTCCcgcg                                    <  2:1308711/150‑1 (MQ=255)
                                                                                                             cTTGCGTCGCCCGGCTTTCCAGAATCTCATGCGCACGCTGCGCATCCTTCAGCGGATATTTCTGCTGCTCGGCGACATCGACCTTAATCACACCGCTGGCAATCAAAGAGAACAGTTCATTACTGGCCTCGGTTAATTCCTCCCGCgtgg                                 <  2:1562500/150‑1 (MQ=255)
                                                                                                             cTTGCGTCGCCCGGCTTTCCAGAATCTCATGCGCACGCTGCGCATCCTTCAGCGGATATTTCTGCTGCTCGGCGACATCGACCTTAATCACACCGCTGGCAATCAAAGAGAACAGTTCATTACTGGCCTCGGTTAATTCCTCCCGCgtgg                                 <  1:362125/150‑1 (MQ=255)
                                                                                                             cTTGCGTCGCCCGGCTTTCCAGAATCTCATGCGCACGCTGCGCATCCTTCAGCGGATATTTCTGCTGCTCGGCGACATCGACCTTAATCACACCGCTGGCAATCAAAGAGAACAGTTCATTACTGGCCTCGGTTAATTCCTCCCGCgtgg                                 <  2:1250026/150‑1 (MQ=255)
                                                                                                              ttGCGTCGCCCGGCTTTCCAGAATCTCATGCGCACGCTGCGCATCCTTCAGCGGATATTTCTGCTGCTCGGCGACATCGACCTTAATCACACCGCTGGCAATCAAAGAGAACAGTTCATTACTGGCCTCGGTTAATTCCTCCCGCgtggt                                <  2:1226887/150‑1 (MQ=255)
                                                                                                              ttGCGTCGCCCGGCTTTCCAGAATCTCATGCGCACGCTGCGCATCCTTCAGCGGATATTTCTGCTGCTCGGCGACATCGACCTTAATCACACCGCTGGCAATCAAAGAGAACAGTTCATTACTGGCCTCGGTTAATTCCTCCCGCgtggt                                <  2:1498654/150‑1 (MQ=255)
                                                                                                               tGCGTCGCCCGGCTTTCCAGAATCTCATGCGCACGCTGCGCATCCTTCAGCGGATATTTCTGCTGCTCGGCGACATCGACCTTAATCACACCGCTGGCAATCAAAGAGAACAGTTCATTACTGGCCTCGGTTAATTCCTCCCGCgtggtg                               <  2:1606756/150‑1 (MQ=255)
                                                                                                                  gTCGCCCGGCTTTCCAGAATCTCATGCGCACGCTGCGCATCCTTCAGCGGATATTTCTGCTGCTCGGCGACATCGACCTTAATCACACCGCTGGCAATCAAAGAGAACAGTTCATTACTGGCCTCGGTTAATTCCTCCCGCGTGGTGata                            <  1:1379293/150‑1 (MQ=255)
                                                                                                                     gCCCGGCTTTCCAGAATCTCATGCGCACGCTGCGCATCCTTCAGCGGATATTTCTGCTGCTCGGCGACATCGACCTTAATCACACCGCTGGCAATCAAAGAGAACAGTTCATTACTGGCCTCGGTTAATTCCTCCCGCGTGGTGATATAg                         >  2:138902/1‑150 (MQ=255)
                                                                                                                            tttCCAGATTCTCATTCGCTGGCTGCGCTTCCTTCAGCGTATATTTCTGCTTCTCGGCGAGATCGACCTTAATCACACCGCTGGCAATCAAAGAGTAGAGTTGATTACTGGCTTCGGTTAATTCGTCCCGCGTGGTGATATAGCCTTGct                  <  2:343124/150‑2 (MQ=255)
                                                                                                                            tttCCAGAATCTCATGCGCACGCTGCGCATCCTTCAGCGGATATTTCTGCTGCTCGGCGACATCGACCTTAATCACACCGCTGGCAATCAAAGAGAACAGTTCATTACTGGCCTCGGTTAATTCCTCCCGCGTGGTGATATAGCCTTGCa                  <  1:678902/150‑1 (MQ=255)
                                                                                                                            tttCCAGAATCTCATGCGCACGCTGCGCATCCTTCAGCGGATATTTCTGCTGCTCGGCGACATCGACCTTAATCACACCGCTGGCAATCAAAGAGAACAGTTCATTACTGGCCTCGGTTAATTCCTCCCGCGTGGTGATATAGCCTTGCa                  >  1:333201/1‑150 (MQ=255)
                                                                                                                            tttCCAGAATCTCATGCGCACGCTGCGCATCCTTCAGCGGATATTTCTGCTGCTCGGCGACATCGACCTTAATCACACCGCTGGCAATCAAAGAGAACAGTTCATTACTGGCCTCGGTTAATTCCTCCCGCGTGGTGATATAGCCTTGCa                  <  1:558314/150‑1 (MQ=255)
                                                                                                                               ccAGAATCTCATGCGCACGCTGCGCATCCTTCAGCGGATATTTCTGCTGCTCTGCGACATCGACCTTAATCACACCGCTGGCAATCAAAGAGAACAGTTCATTACTGGCCTCGGTTAATTCCTCCCGCGTGGTGATATAGCCTTGCAggg               <  2:1458288/150‑1 (MQ=255)
                                                                                                                                 aGAATCTCATGCGCACGCTGCGCATCCTTCAGCGGATATTTCTGCTGCTCGGCGACATCGACCTTAATCACACCGCTGGCAATCAAAGAGAACAGTTCATTACTGGCCTCGGTTAATTCCTCCCGCGTGGTGATATAGCCTTGCAGGGaa             >  2:166760/1‑150 (MQ=255)
                                                                                                                                      ctcATGCGCACGCTGCGCATCCTTCAGCGGATATTTCTGCTGCTCGGCGACATCGACCTTAATCACACCGCTGGCAATCAAAGAGAACAGTTCATTACTGGCCTCGGTTAATTCCTCCCGCGTGGTGATATAGCCTTGCAGGGAAGGGCg        >  1:1213226/1‑150 (MQ=255)
                                                                                                                                      aggaTGCGCAGCGTGCGCATCCTTCAGCGGATATTTCTGCTGCTCGGCGACATCGACCTTAATCACACCGCTGGTAATCAAAGAGAACAGTTCATTACTGGCCTCGGTTAATTCCTCCCGCGTGGTGATATAGCCTTGCAGGGAAGGGCg        >  2:325228/4‑150 (MQ=255)
                                                                                                                                        cATGCGCACGCTGCGCATCCTTCAGCGGATATTTCTGCTGCTCGGCGACATCGACCTTAATCACACCGCTGGCAATCAAAGAGAACAGTTCATTACTGGCCTCGGTTAATTCCTCCCGCGTGGTGATATAGCCTTGCAGGGAAGGGCgtg      >  2:458786/1‑150 (MQ=255)
                                                                                                                                         aTGCGCACGCTGCGCATCCTTCAGCGGATATTTCTGCTGCTCGGCGACATCGACCTTAATCACACCGCTGGCAATCAAAGAGAACAGTTCATTACTGGCCTCGGTTAATTCCTCCCGCGTGGTGATATAGCCTTGCAGGGAAGGGCgtgt     <  1:732723/150‑1 (MQ=255)
                                                                                                                                          tGCGCACGCTGCGCATCCTTCAGCGGATATTTCTGCTGCTCGGCGACATCGACCTTAATCACACCGCTGGCAATCAAAGAGAACAGTTCATTACTGGCCTCGGTTAATTCCTCCCGCGTGGTGATATAGCCTTGCAGGGAAGGGCGTGTc    >  2:675949/1‑150 (MQ=255)
                                                                                                                                           gcgcACGCTGCGCATCCTTCAGCGGATATTTCTGCTGCTCGGCGACATCGACCTTAATCACACCGCTGGCAATCAAAGAGAACAGTTCATTACTGGCCTCGGTTAATTCCTCCCGCGTGGTGATATAGCCTTGCAGGGAAGGGCGTGTca   >  2:1595631/1‑150 (MQ=255)
                                                                                                                                            cccACGCTGCGCATCCTTCAGCGGATATTTCTGCTGCTCGGCGACATCGACCTTAATCACACCGCTGGCAATCAAAGAGAACAGTTCATTACTGGCCTCGGTTAATTCCTCCCGCGTGGTGATATAGCCTTGCAGGGAAGGGCGTGTcac  >  1:422212/3‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
CGCCTCTAACAGATTCATCGTGCTGTACCCTACATACAGCCGAACTATAAAAAGAAAGGGCTTCCCAGGTGGAAGCCCTATTTCTTTTATGGAATCAGCAGGCTGGAACCTTGCGTCGCCCGGCTTTCCAGAATCTCATGCGCACGCTGCGCATCCTTCAGCGGATATTTCTGCTGCTCGGCGACATCGACCTTAATCACACCGCTGGCAATCAAAGAGAACAGTTCATTACTGGCCTCGGTTAATTCCTCCCGCGTGGTGATATAGCCTTGCAGGGAAGGGCGTGTCAC  >  NC_000913/4263162‑4263451

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: