Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,424,327 C→T 5.2% S61S (TCG→TCA hisP ← histidine ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,424,3270CT 5.2% 227.8 / 2.3 77S61S (TCG→TCAhisPhistidine ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base C (43/30);  new base T (2/2);  total (45/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.55e-01

CACCGTCATATGGCTCCAGAGATTGAAGTGCTGGAATACCATCGTCAGGCGTGTGCGCAGTAAGCGCAGTTGATTTTTATCGGCGACTTTGAGTTGACCGTCTTTGTCGCGCACCAGATTGATCGTCTGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTACTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACATCTCCGGCATTCGCTTGCAGTGATACCCCTTTCAGCACTTCATGTTCGCCGTAGCGT  >  NC_000913/2424183‑2424471
                                                                                                                                                |                                                                                                                                                
cACCGTCATATGGCTCCAGAGATTGAAGTGCTGGAATACCATCGTCAGGCGTGTGCGCAGTAAGCGCAGTTGATTTTTATCGGCGACTTTGAGTTGACCGTCTTTGTCGCGCACCAGATTGATCGTCTGGCCATTGACCACGATCGAccc                                                                                                                                             >  1:1215831/1‑150 (MQ=255)
 aCCGTCATATGGCTCCAGAGATTGAAGTGCTGGAATACCATCGTCAGGCGTGTGCGCAGTAAGCGCAGTTGATTTTTATCGGCGACTTTGAGTTGACCGTCTTTGTCGCGCACCAGATTGATCGTCTGGCCATTGACCACGATCGAccca                                                                                                                                            >  1:584359/1‑149 (MQ=255)
  ccGTCATATGGCTCCAGAGATTGAAGTGCTGGAATACCATCGTCAGGCGTGTGCGCAGTAAGCGCAGTTGATTTTTATCGGCGACTTTGAGTTGACCGTCTTTGTCGCGCACCAGATTGATCGTCTGGCCATTGACCACGATCGACCCtt                                                                                                                                           >  2:7150/1‑150 (MQ=255)
    gTCATATGGCTCCAGAGATTGAAGTGCTGGAATACCATCGTCAGGCGTGTGCGCAGTAAGCGCAGTTGATTTTTATCGGCGACATTGAGTTGACCGTCTGTGTCGCGCACCAGATTGATCGTCTGGCCATTGACAACGATCGACCCTTCa                                                                                                                                         >  1:616815/1‑150 (MQ=255)
      cATATGGCTCCAGAGATTGAAGTGCTGGAATACCATCGTCAGGCGTGTGCGCAGTAAGCGCAGTTGATTTTTATCGGCGACTTTGAGTTGACCGTCTTTGTCGCGCACCAGATTGATCGTCTGGCCATTGACCACGATCGACCCTTCAcg                                                                                                                                       >  2:671980/1‑149 (MQ=255)
      aatatGGCTCCAGAGATTGAAGTGCTTGAATACCATCGTCAGGCGTGTGCGCAGTAAGCGCAGTTGATTTTTATCGGCGACTTTGAGTTGACCGTCTTTGTCGCGCACCAGATTGATCGTCTGGCCATTGACCACGATCGACCCTTCACt                                                                                                                                       <  1:960225/149‑1 (MQ=255)
             cTCCAGAGATTGAAGTGCTGGAATACCATCGTCAGGCGTGTGCGCAGTAAGCGCAGTTGATTTTTATCGGCGACTTTGAGTTGACCGTCTTTGTCGCGCACCAGATTGATCGTCTGGCCATTGACCACGATCGACCCTTCACTCGGtttt                                                                                                                                >  1:420934/1‑150 (MQ=255)
                cAGAGATTGATGTGCTGGAATACCATCGTCAGGCGTGTGCGCAGTAAGCGCAGTTGATTTTTATAGGCGACTTTGAGTTGACCGTCTTTGTCGCGCACCAGATTGATCGTCTGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCg                                                                                                                             <  1:124332/150‑1 (MQ=255)
                cAGAGATTGAAGTGCTGGAATACCATCGTCAGGCGTGTGCGCAGTAAGCGCAGTTGATTTTTATCGGCGACTTTGAGTTTACCGTCTTTGTCGCGCAACAGATTGATCGTCTGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCg                                                                                                                             >  2:698045/1‑150 (MQ=255)
                cAGAGATTGAAGTGCTGGAATACCATCGTCAGGCGTGTGCGCAGTAAGAGCAGTTGATTTTTATAGGCGACTATGAGTTGACCGTCTTTGTCGCGCACCAGATAGATCGTCTGGCCATTTACCACGATCGACACTTCACTCGGTTTTTCg                                                                                                                             >  1:243094/1‑150 (MQ=255)
                  gcgaTTGAAGTGCTGGAATACCATCGTCAGGCGTGTGCGCAGTAAGCGCAGTTGATTTTTATCGGCGACTTTGAGTTAACCGTCTTTGTCGCGCACCAGATTGATCGTCTGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAg                                                                                                                           <  2:10364/148‑1 (MQ=255)
                        gAAGTGCTGGAATACCATCGTCAGGCGTGTGCGCAGTAAGCGCAGTTGATTTTTATCGGCGACTTTGAGTTGACCGTCTTTGTCGCGCACCAGATTGATCGTCTGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGtt                                                                                                                     >  1:677839/1‑150 (MQ=255)
                        gAAGTGCTGGAATACCATCGTCAGGCGTGTGCGCAGTAAGCGCAGTTGATTTTTATCGGCGACTTGGAGTTGACCGTCTTTGTCGCGAACCAGATTGATCGTCTGGCCATTGACCACGATCGACCCTTCACTCGGGTTTTCGAGGAAGtt                                                                                                                     <  1:330320/150‑1 (MQ=255)
                         aaGTGCTGGAATACCATCGTCAGGCGTGTGCGCAGTAAGCGCAGTTGATTTTTATCGGCGACTTTGAGTTGACCGTCTTTGTCGCGCACCAGATTGATCGTCTGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTa                                                                                                                    >  1:955643/1‑150 (MQ=255)
                             gCTGGAATCCCATCGGCAGGCGTGTGCGCAGTAAGCGCAGTTGATTTATATCGGCGACTTTGAGTTGACCGTCTTTGTCACGCACCAGATTGATCGTCTGGCCATTGACCACGATCGACCCTTCACTAGGATTTTCGACGAACTTAATGc                                                                                                                >  2:362502/1‑150 (MQ=255)
                                  aaTACCATCGTCAGGCGTATGCGCAGTAAGCGCAGTTGATTTTTATCGGCGACTTTGAGTTGACCGTCTTTGTCGCGCACCAGATTGATCGTCTGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAgcgc                                                                                                           <  2:690001/150‑1 (MQ=255)
                                     aCCATCGTCAGGCGTGTGCGCAGTAAGCGCAGTTGAGTTTTATCGGCGACTTTGAGTTGACCGTCTTTGTCGCGCCCCTGATTGATCGTCTGGCCATTGCCCACGATCGAACCTTCACTCGGTTTTTCGCGGATGTTAATGGAgcgctgt                                                                                                        >  2:1126429/1‑147 (MQ=255)
                                           gTCAGGCGTGTGCGCAGTAAGCGCAGTTGATTTTTATCGGCGACTTTGAGTTGACCGTCTTTGTCGCGCACCAGATTGATCGTCTGGCCATTGACCACGATCGACCCTTCACTCGGGTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTa                                                                                                  >  1:52697/1‑150 (MQ=255)
                                             catgCGTGTTCGCAGTAAGCGCAGTTGATTTTTATCGGCGACTTTGAGTTGACCGTCTTTGTCTCGCACCAGATTGATCGTCTGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTACt                                                                                                <  1:233389/147‑1 (MQ=255)
                                             cAGGCGTGTGCGCAGTAAGCGCAGTTGATTTTTATCGGCGACTTTGAGTTGACCGTCTTTGTCGCGCACCAGATTGATCGTCTGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCTGCGCAGAAAGGTACt                                                                                                >  1:1088785/1‑150 (MQ=255)
                                                 cGTGTGCGCAGTAAGCGCAGTTGATTTTTATCGGCGACTTTGAGTTGACCGTCTTTGTCGCGCACCAGATTGATCGTCTGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTACTTTTc                                                                                            >  2:1031988/1‑150 (MQ=255)
                                                   tgtgCGCAGTAAGCGCAGTTGATTTTTATCGGCGACTTTGAGTTGACCGTCTTTGTCGCGCACCAGATTGATCGTCTGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTACTTTTccc                                                                                          >  1:1063854/1‑150 (MQ=255)
                                                    gtgCGCAGTAAGCGCAGTTGATTTTTATCGGCGACTTTGAGTTGACCGTCTTTGTCGCGCACCAGATTGATCGTCTGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTACTTTTcccc                                                                                         <  2:1171816/150‑1 (MQ=255)
                                                     tgCGCAGTAAGCGCAGTTGATTTTTATCGGCGACTTTGAGTTGACCGTCTTTGTCGCGCACCAGATTGATCGTCTGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGATGTTAATGCAGCGCAGAAAGGTACTTTTCCCCg                                                                                        >  1:1127804/1‑150 (MQ=255)
                                                         cAGTAAGCGCAGTTGATTTTTATCGGCGACTTTGAGTTGACCGTCTTTGTCGCGCACCAGATTGATCGTCTGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTACTTTTCCCCGATcc                                                                                    >  2:1104128/1‑150 (MQ=255)
                                                         cAGTAAGCGCAGTTGATTTTTATCGGCGACTTTGAGTTGACCGTCTTTGTCGCGCACCAGATTGATCGTCTGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTACTTTTCCCCGATcc                                                                                    <  2:600077/150‑1 (MQ=255)
                                                         cAGTAAGCGCAGTTGATTTTTATCGGCGACTTTGAGTTGACCGTCTTTGTCGCGCACCAGATTGATCGTCTGGCCATTGACCACGATCGACCCTTCACTCGGTGTTTGGAGGAAGTTAATGCAGCGCAGAAAGGTACTTTTCCCCGATcc                                                                                    >  1:1082863/1‑150 (MQ=255)
                                                          attAAGCGCAGTTGATTTTTATCGGCGACTTTGAGTTGACCGTCTTTGTCGCGCACCAGATTGATCGTCTGGCCATTGACCACTATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTACTTTTCCCCGATccc                                                                                   <  2:245927/148‑1 (MQ=255)
                                                           gTAAGCGCAGTTGATTTTTATCGGCGACTTTGAGTTGACCGTCTTTGTGGCGCACAATATTGATCGTCTGGCCATTGACCACGATCGAGCCTTCACTCGGATTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTACTTTTCCCCGTTCCcg                                                                                  >  1:188812/1‑150 (MQ=255)
                                                           gTAAGCGCAGTTGATTTTTATCGGCGACTTTGAGTTGACCGTCTTTGTCGCGCACCAGATTGCTCGTCTGGCCATTGACCACGATCGACCCTTCACTCGGGTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTACTTTTCCCCGATCCcg                                                                                  >  1:86849/1‑150 (MQ=255)
                                                             aaGCGCAGTTGATTTTTATCGGCGACTTTGAGTTGACCGTCTTTGTCGCGCACCAGATTGATCGTCTGGCCATTGACCACGATCGCCCCTTCACTCGGGTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTACTTTTCCCCGATCCcgac                                                                                <  2:758858/150‑1 (MQ=255)
                                                             aaGCGCAGTTGATTTTTATCGGCGACTTTGAGTTGACCGTCTTTGTCGCGCACCAGATTGATCGTCTGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTACTTTTCCCCGATCCcgac                                                                                <  2:175226/150‑1 (MQ=255)
                                                                     ttGATTTTTATCGGCGACTTTGAGTTGACCGTCTTTGTCGCGCACCAGATTGATCGTCTGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTACTTTTCCCCGATCCCGACGATCCgat                                                                        <  1:1220815/150‑1 (MQ=255)
                                                                      tGATTTTTATCGGCGACTTTGAGTTGACCGTCTTTGTCGCGCACCAGATTGATCGTCTGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTACTTTTCCCCGATCCCGACGATCCgatg                                                                       >  1:1184706/1‑150 (MQ=255)
                                                                         tttttATCGGCGACTTTGAGTTGACCGTCTTTGTCGCGCACCAGATTGATCGTCTGGCCATTGACAACGATCGCCCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTACTTTTCACCGATCCCGACGCTCCgaagatt                                                                    >  2:913958/1‑149 (MQ=255)
                                                                             tATCGGCGAGATTGAGTTCACCGTCTTTGTCGCGCACCAGATTGAGCGTCTGGCCATTGACCACGATCGACCCATCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTACTTTTCCCCGAACCCGACGATCCGATGATGCtta                                                                <  2:574669/150‑1 (MQ=255)
                                                                              aTCGGCGACTTTGAGTTGACCGTCTTTGTCGCGCACCAGATTGATCGTCTGGCCATTGACCACGATTGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTACTTTTCCCCGATCCCGACGATCCGATGATGCt                                                                  >  1:1233838/1‑147 (MQ=255)
                                                                              aTCGGCGACTTTGAGTTGACCGTCTTTGTCGCGCACCAGATTGATCGTCTGGCCATTGACCACGATTGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTACTTTTCCCCGATCCCGACGATCCGATGATGCt                                                                  <  2:1233838/147‑1 (MQ=255)
                                                                              aTCGGCGACTTTGAGTTGACCGTCTTTGTCGCGCACCAGATTGATCGTCTGGCCATTGACCACGATTGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTACTTTTCCCCGATCCCGACGATCCGATGATGCt                                                                  >  1:418599/1‑147 (MQ=255)
                                                                              aTCGGCGACTTTGAGTTGACCGTCTTTGTCGCGCACCAGATTGATCGTCTGGCCATTGACCACGATTGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTACTTTTCCCCGATCCCGACGATCCGATGATGCt                                                                  <  2:418599/147‑1 (MQ=255)
                                                                               tCGGCGACTTTGAGTTGACCGTCTTTGTCGCGCACCAGATTGATCGTCTGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGAGCAGAAAGGTACTTTTCCCCGATCCCGACGATCCGATGATGCttatt                                                              >  1:501857/1‑150 (MQ=255)
                                                                               tCGGCGACTTTGAGTTGACCGTCTTTGTCGCGCACCAGATTGATCGGCTGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTACTTTTCCCCGATCCCGACGATCCGATGATGCttatt                                                              <  2:1073754/150‑1 (MQ=255)
                                                                                  gCGACTTTGAGTTGTCCGTCTTTGTCGCGCAGCAGATAGATCGTCTGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTACATTTCCCCGATCCCGACGATCCGATGATGCTTATTACa                                                           <  2:361035/150‑1 (MQ=255)
                                                                                  gCGACTTTGAGTTGACCGTCTTTGTCGCGCACCAGATTGATCGTCTGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTACTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACa                                                           <  2:349905/150‑1 (MQ=255)
                                                                                       tttGAGTTGACCGTCTTTGTCGCGCACCAGATTGATCGTCTGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTACTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACATCTcc                                                      <  2:318933/150‑1 (MQ=255)
                                                                                               gACCGTCTTTGTCGCGCACCAGATTGATCGTCTGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAATGGTACTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACATCTCCGGCATTCg                                              >  1:95899/1‑150 (MQ=255)
                                                                                                 ccGTCTTTGTCGCGCACCAGATTGATCGTCTGGCCATTGACCACGATCGACCATTCACTCGGTTTTTAGAGGAAGTTAATGCAGCGCAGAAAGGTACTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACATCTACGGCATTCGCt                                            >  1:628552/1‑150 (MQ=255)
                                                                                                  cGTCTTTGTCGCGCACCAGATTGATCGTCTGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTACTTTTCCCCGATCCAGACGATCCGATGATGCTTATTACATCTCCGGCATTCGCtt                                           >  1:466288/1‑150 (MQ=255)
                                                                                                    tCTTTGTCGCGCACCAGATTGATCGGCTGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTACTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACATCTCCGGCATTCGCTTGc                                         <  1:782211/150‑1 (MQ=255)
                                                                                                      tttGTCGCGCACCAGATTGATCGTCTGGCCATTGCCCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCATCTCAGAAAGGTACTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACATCTCCGGCATTCGCTTGCAg                                       <  1:224721/150‑1 (MQ=255)
                                                                                                      tttGTCGCGCACCAGATTGATCGTCTGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTACTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACATCTCCGGCATTCGCTTGCAg                                       <  1:1031988/150‑1 (MQ=255)
                                                                                                      gttGTCGCGCACCAGATTGATCGTCTGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTACTATTCCCCGATCCCGACGATCCGATGATGCTTCTTACATCTCCGGCATTCGCTTGCAg                                       <  1:814700/149‑1 (MQ=255)
                                                                                                       ttGTCGCGCACCAGATTGATCGTCTGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTACTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACATCTCCGGCATTCGCTTGCAGt                                      >  2:96745/1‑150 (MQ=255)
                                                                                                        tGTCGCGCACCAGATTGATCGTCTGGCCATTGACCACGATCGACCCTTCCCTCGGTCTTTCGAGGAAGTTAATGCAGCGCAGAACGGTACTTTTCCCCGATCCCGACTATCCGATGATGCTTATTACATCTCCGGCATTCGCTTGCAGTg                                     <  1:864461/150‑1 (MQ=255)
                                                                                                        ggTCGCGCACCAGATTGATCGTCTGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGGACTTTTCCCCGATCCCGCCGATCCGATGATGCTTATTACATCTCCGGCATTCGCTGGCAGTg                                     <  1:917736/149‑1 (MQ=255)
                                                                                                          tCGCGCACCAGATTGATCGTCTGGCCATTGACCACGATCGCCCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTACTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACATCTCCGGCATTCGCTTGCAGTGAt                                   >  2:263131/1‑150 (MQ=255)
                                                                                                          tCGCGCACCAGATTGATCGTCTGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTACTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACATCTCCGGCATTCGCTTGCAGTGAt                                   >  2:1223924/1‑150 (MQ=255)
                                                                                                           cgcgCACCAGTTTGATCGTCTGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTACTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACATCTCCGGCATTCGCTTGCAGTGATa                                  <  2:762530/150‑1 (MQ=255)
                                                                                                           cgcgCACCAGATTGATCGTCTGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTACTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACATCTCCGGCATTCGCTTGCAGTGATa                                  >  2:111618/1‑150 (MQ=255)
                                                                                                             cgcACCAGATTGATCGTCTGGCCATTGCCCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTACTTTTCCCCGATCCCGCCGCTCCGATGATGCTTATTACATCTCCGGCATTCGCTTGCAGTGATccc                                >  2:523542/1‑150 (MQ=255)
                                                                                                             cgcACCAGATTGATCGTCTGGCCATTGACCCCGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTACTGCAGCGCAGAAAGGTACTTTTCCCCGAGCCCGACGATCCGATGATGCTTATTCCCTCTCCGGCATTCGCTTGCAGTGATAcc                                >  1:594958/1‑150 (MQ=255)
                                                                                                             cgcACCAGATTGATCGTCTGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTACTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACATCTCCGGCATTCGCTTGCAGTGATAcc                                <  2:591737/150‑1 (MQ=255)
                                                                                                                 ccAGATTGATCGTCAGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTACTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACATCTCCGGCATTCGCTTGCAGTGATACCCCtt                            >  1:646079/1‑150 (MQ=255)
                                                                                                                  cAGATTGATCGTCTGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTACTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACATCTCCGGCATTCGCTTGCAGTGATACCCCttt                           <  1:1113659/150‑1 (MQ=255)
                                                                                                                   aGATTGATCGTCTGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGTAAGTTAATGCATCGCAGAAAGGTACTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACATCTCCGGCATTCGCTTGCAGTGATACCCCTTTc                          <  2:29559/150‑1 (MQ=255)
                                                                                                                     aTTGATCGTCTGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTACTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACATCTCCGGCATTCGCTTGCAGTGATACCCCTTTCAg                        >  1:988420/1‑150 (MQ=255)
                                                                                                                       tGATCGTCTGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTACTTTTCCCCGATCCCGACGATCCGATGATTCTAATTACATCTCCGGCATTCGCTTGCAGTGATACCCCTTTCAGCa                      <  1:628736/150‑1 (MQ=255)
                                                                                                                            gTCTGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTACTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACATCTCCGGCATTCGCTTGCAGTGATACCCCTTTCAGCACTTCa                 <  2:1215831/150‑1 (MQ=255)
                                                                                                                               tGGCCATTGACCACGATCGACACTTCAATCGGTTTTTCGAGGAAGTTAATGCAGATCAGAAAGGAACTTTTCCCCGATCCCGACGATCCGATGAAGCTTATTACATCTCCGACATTCGCTTGCAGTGATACCCATTACACCACTTCTTGt              >  2:457907/1‑150 (MQ=255)
                                                                                                                                 gCCATTGACCACGATCGACCCTTCACTCGGTTTAGCGAGGTAGTTAATGCAGCGCAGAAAGTTACTTTTCCCCGATCACGACGATCCGATGATGCTTATTACATCTCCGGCATTCGCTTGCAGTGATACCCCTTTCAGCACTTCATGTTc            <  2:1009824/150‑1 (MQ=255)
                                                                                                                                  ccATTGACCACGATCGCCCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTACTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACATCTCCGGCATTCGCTTGCAGTGAGACCCCTTTCAGCACTTCATGTTCg           >  2:114521/1‑150 (MQ=255)
                                                                                                                                     ttGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTACTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACATCTCCGGCATTCGCTTGCAGTGATACCCATTTCAGCACTACATGTTCGCCg        >  2:1030731/1‑150 (MQ=255)
                                                                                                                                       gACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTACTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACATCTCCGGCATTCGCTTGCAGTGATACCCCTTTCAGCACTTCATGTTCGCCGTa      >  2:150357/1‑150 (MQ=255)
                                                                                                                                        aCCACGCTCGACCCTTCACTCGGGTTTTCGAGGAAGGTAATGCAGCGCAGAAAGGTACTTTTCCCCGCTCCCGCCGATCCGATGATGCTGATTACATCTCCGGCATTCGCTTGCAGTGATACCCCTTTCCGCACTGCCTGTTCGCCGTAg     >  2:354602/1‑150 (MQ=255)
                                                                                                                                        aCCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTACTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACATCTCCGGCATTCGCTTGCAGTGATACCCCTTTCAGAACTTCATGTTCGCCGTAg     >  2:376243/1‑150 (MQ=255)
                                                                                                                                         ccACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTACTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACATCTCCGGCATTCGCTTGCAGTGATACCCCTTTCAGCACTTCATGTTCGCCGTAGc    >  2:246436/1‑150 (MQ=255)
                                                                                                                                           aCGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTACTTTTCCCCGAACCCGACGATCCGATGATCCTTATTACATCTCCGGCATTCGCTTGCAGTGATACCCCTTTCAGCACTTCATGTTCGCCGTATCGt  >  1:624513/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                
CACCGTCATATGGCTCCAGAGATTGAAGTGCTGGAATACCATCGTCAGGCGTGTGCGCAGTAAGCGCAGTTGATTTTTATCGGCGACTTTGAGTTGACCGTCTTTGTCGCGCACCAGATTGATCGTCTGGCCATTGACCACGATCGACCCTTCACTCGGTTTTTCGAGGAAGTTAATGCAGCGCAGAAAGGTACTTTTCCCCGATCCCGACGATCCGATGATGCTTATTACATCTCCGGCATTCGCTTGCAGTGATACCCCTTTCAGCACTTCATGTTCGCCGTAGCGT  >  NC_000913/2424183‑2424471

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