Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,987,174 C→G 11.4% D237H (GAT→CAT)  hemY ← protein HemY

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,987,1740CG11.4% 264.1 / 22.1 88D237H (GAT→CAT) hemYprotein HemY
Reads supporting (aligned to +/- strand):  ref base C (38/40);  new base G (4/6);  total (42/46)
Fisher's exact test for biased strand distribution p-value = 7.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.74e-01

GATAATTTGCTGGGCAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCGCCAGACGTAATACTTC  >  NC_000913/3987025‑3987321
                                                                                                                                                     |                                                                                                                                                   
gatAATTTGCTGGGCAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATc                                                                                                                                                     <  1:234748/150‑1 (MQ=255)
gatAATTTGCTGGGCAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATc                                                                                                                                                     <  1:1477547/150‑1 (MQ=255)
 atAATTTGCTGGGCAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATcc                                                                                                                                                    <  2:1557974/150‑1 (MQ=255)
 atAATTTGCTGGGCAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATcc                                                                                                                                                    >  1:515485/1‑150 (MQ=255)
  tAATTTGCTGGGCAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCa                                                                                                                                                   <  2:546840/150‑1 (MQ=255)
    aTTTGCTGGGCAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATTACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATc                                                                                                                                                 <  2:1503594/150‑1 (MQ=255)
    aTTTGCTGGGCAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATc                                                                                                                                                 <  2:852464/150‑1 (MQ=255)
       tGCTGGGCAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAgg                                                                                                                                              <  1:267221/150‑1 (MQ=255)
       tGCTGGGCAGTATCATGATCGGCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAgg                                                                                                                                              <  1:1632938/150‑1 (MQ=255)
        gCTGGGCAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTAATCCATCAGGc                                                                                                                                             <  2:940493/150‑1 (MQ=255)
           gggCAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATGGATCAGGcgcg                                                                                                                                          <  2:1226271/150‑4 (MQ=255)
           gggCAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATGGATCAGGcgcg                                                                                                                                          <  2:1387893/150‑4 (MQ=255)
           gggCAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATGGATCAGGcgcg                                                                                                                                          <  1:1185569/150‑4 (MQ=255)
            ggCAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATGGATCAGGcgcgt                                                                                                                                         <  1:144868/150‑5 (MQ=255)
             gCAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATGGATCAGGcgcgtg                                                                                                                                        <  2:104334/150‑6 (MQ=255)
               aGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATGGATCAGGcgcgtgcc                                                                                                                                      <  1:970654/150‑8 (MQ=255)
                gTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCAt                                                                                                                                     >  1:709065/1‑150 (MQ=255)
                gTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCAt                                                                                                                                     <  2:483450/150‑1 (MQ=255)
                      tGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATcc                                                                                                                                                    <  2:86287/129‑1 (MQ=255)
                      tGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATcc                                                                                                                                                    >  1:86287/1‑129 (MQ=255)
                       gATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTATCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCtgtt                                                                                                                              >  1:917062/1‑150 (MQ=255)
                            tCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTcc                                                                                                                         <  2:486131/150‑1 (MQ=255)
                             cacaTTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCa                                                                                                                        <  1:1201089/150‑1 (MQ=255)
                               caTTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGc                                                                                                                      >  1:1410734/1‑150 (MQ=255)
                               caTTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGc                                                                                                                      <  2:1129518/150‑1 (MQ=255)
                                 ttCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCAt                                                                                                                    <  2:515485/150‑1 (MQ=255)
                                   cAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTg                                                                                                                  >  1:696573/1‑150 (MQ=255)
                                      tAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCAc                                                                                                               <  2:303558/150‑1 (MQ=255)
                                        aGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGa                                                                                                             >  2:1322072/1‑150 (MQ=255)
                                         gATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGAt                                                                                                            <  1:638964/150‑1 (MQ=255)
                                            gTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGtt                                                                                                         >  2:859850/1‑150 (MQ=255)
                                                     ttGCCACCTGCATCGCTAGGTGTTTTGGCGTTTTGCTTCTTTGGTTTTTCCACGAGTTACGCATAGCTTCGCTACCGTTATCTTCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTGCAGCATTGCACGATGTTCTtcatctc                                                                                                <  2:1468507/150‑3 (MQ=255)
                                                          aCCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACa                                                                                           <  2:1467801/150‑1 (MQ=255)
                                                            cTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATg                                                                                         <  1:672832/150‑1 (MQ=255)
                                                            cTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATg                                                                                         >  2:1370039/1‑150 (MQ=255)
                                                              gCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATggg                                                                                       <  2:917062/150‑1 (MQ=255)
                                                                   gCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTc                                                                                  <  1:1322072/150‑1 (MQ=255)
                                                                     tACCTGATGACGCGTTTTCAGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGc                                                                                >  2:1184448/1‑150 (MQ=255)
                                                                      aCCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATggg                                                                                       <  2:848004/142‑1 (MQ=255)
                                                                      aCCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATggg                                                                                       >  1:848004/1‑142 (MQ=255)
                                                                       ccTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCa                                                                              >  2:453120/1‑150 (MQ=255)
                                                                          gatgaCGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTg                                                                           >  1:1516958/1‑150 (MQ=255)
                                                                            tgaCGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGAt                                                                         <  1:1557700/150‑1 (MQ=255)
                                                                             gaCGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATg                                                                        >  2:728773/1‑150 (MQ=255)
                                                                             gaCGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATg                                                                        >  2:70098/1‑150 (MQ=255)
                                                                               cgcgTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTcc                                                                                                                         >  2:1342173/1‑99 (MQ=255)
                                                                               cgcgTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTcc                                                                                                                         <  1:1342173/99‑1 (MQ=255)
                                                                               cgcgTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATggg                                                                      <  1:859850/150‑1 (MQ=255)
                                                                               cgcgTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATggg                                                                      <  1:202301/150‑1 (MQ=255)
                                                                               cgcgTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATggg                                                                      <  1:1486878/150‑1 (MQ=255)
                                                                                  gTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGata                                                                   <  1:487107/150‑1 (MQ=255)
                                                                                       ccGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATc                                                              >  2:97154/1‑150 (MQ=255)
                                                                                       ccGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATc                                                              >  2:644251/1‑150 (MQ=255)
                                                                                        cGGCTTTGGTTTTTCCACCAGTTTCGCAAACCTTCGCTACCGTTATCTGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCAAGCATTGCACGATGTTCTTCATGACCAACATGGGCTTTCGCCATTGATGGGATAATTTcc                                                             <  2:710588/150‑1 (MQ=255)
                                                                                         ggCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCca                                                            >  2:214502/1‑150 (MQ=255)
                                                                                           cTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCAt                                                                             <  2:91912/131‑1 (MQ=255)
                                                                                           cTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCAt                                                                             >  1:91912/1‑131 (MQ=255)
                                                                                            tttGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCcagca                                                         >  1:683038/1‑150 (MQ=255)
                                                                                               ggTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGTCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCg                                                      <  1:451199/150‑1 (MQ=255)
                                                                                               ggTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATgg                                                                       >  2:839392/1‑133 (MQ=255)
                                                                                               ggTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATgg                                                                       <  1:839392/133‑1 (MQ=255)
                                                                                               ggTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCg                                                      >  1:1276080/1‑150 (MQ=255)
                                                                                               ggTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCg                                                      >  2:1285118/1‑150 (MQ=255)
                                                                                               ggTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCg                                                      >  1:320323/1‑150 (MQ=255)
                                                                                                gTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGa                                                     >  1:12643/1‑150 (MQ=255)
                                                                                                gTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGa                                                     >  1:1577461/1‑150 (MQ=255)
                                                                                                  ttttCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAAc                                                   <  1:258888/150‑1 (MQ=255)
                                                                                                  ttttCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAAc                                                   >  2:1248413/1‑150 (MQ=255)
                                                                                                  ttttCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAAc                                                   <  1:496979/150‑1 (MQ=255)
                                                                                                      ccaccaGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCa                                               <  1:1285118/150‑1 (MQ=255)
                                                                                                            gTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCAcc                                         <  1:786506/150‑1 (MQ=255)
                                                                                                               aCGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCtgt                                      >  2:227856/1‑150 (MQ=255)
                                                                                                                cGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCtgtg                                     >  1:649354/1‑150 (MQ=255)
                                                                                                                  cAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCg                                   >  2:442633/1‑150 (MQ=255)
                                                                                                                       cTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGatat                              <  2:1410734/150‑1 (MQ=255)
                                                                                                                         tCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATAc                            >  1:1225666/1‑150 (MQ=255)
                                                                                                                         tCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATAc                            >  1:988298/1‑150 (MQ=255)
                                                                                                                            cTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGcc                         <  2:1272486/150‑1 (MQ=255)
                                                                                                                                cGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGtt                     <  1:70098/150‑1 (MQ=255)
                                                                                                                                 gTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTc                    <  2:683038/150‑1 (MQ=255)
                                                                                                                                   tATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCg                  <  1:1248413/150‑1 (MQ=255)
                                                                                                                                      cGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCGCCa               >  2:990051/1‑150 (MQ=255)
                                                                                                                                       tggattggcCTGATGGATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCGCCAg              >  1:1604263/8‑150 (MQ=255)
                                                                                                                                       ggCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCGCCAg              >  2:1129002/1‑150 (MQ=255)
                                                                                                                                       ggCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCGCCAg              >  1:987472/1‑150 (MQ=255)
                                                                                                                                        ggattggcCTGATGGATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCGCCAGa             >  2:621082/7‑150 (MQ=255)
                                                                                                                                        gCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCGCCAGa             <  2:137955/150‑1 (MQ=255)
                                                                                                                                         gattggcCTGATGGATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCGCCAGAc            >  2:37447/6‑150 (MQ=255)
                                                                                                                                         cACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCGCCAGAc            >  1:1524215/1‑150 (MQ=255)
                                                                                                                                           ttggcCTGATGGATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCGCCAGACGt          >  1:707761/4‑150 (MQ=255)
                                                                                                                                            gcgcCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCGCCAGACGTa         <  2:297701/150‑1 (MQ=255)
                                                                                                                                            gcgcCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCGCCAGACGTa         >  1:417981/1‑150 (MQ=255)
                                                                                                                                            gcgcCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCAGTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCGCCAGACGTa         >  1:1111580/1‑150 (MQ=255)
                                                                                                                                                   aTCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCGCCAGACGTAATACTTc  >  1:959087/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GATAATTTGCTGGGCAGTATCATGATCGTCACATTCAATAAGATGTTCCGCCATTGCCACCTGCAACGCTACCTGATGACGCGTTTTCCGGCTTTGGTTTTTCCACCAGTTACGCAAACCTTCGCTACCGTTATCGGCACGCGCCTGATCCATCAGGCCAATCCATGCCTGTTGTTCCAGCATTGCACGATGTTCTTCATCACCAACATGGGCTTTCGCCATTGATGGGATAATATCCAGCAGCGAACTCCATGCACCTGTGCGGATATACGCCTGTTCCGCCAGACGTAATACTTC  >  NC_000913/3987025‑3987321

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: ‑