Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,896,352 A→C 12.2% V86V (GTT→GTG)  adeP ← adenine:H(+) symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,896,3520AC12.2% 199.1 / 21.1 77V86V (GTT→GTG) adePadenine:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base A (40/27);  new base C (6/3);  total (47/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.51e-01

GCGAACGCGGAAAATCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAGGTGG  >  NC_000913/3896211‑3896500
                                                                                                                                             |                                                                                                                                                    
gCGAACGCGGAAAATCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAg                                                                                                                                              >  1:1315390/1‑150 (MQ=255)
 cGAACGCGGATAATCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAgg                                                                                                                                             >  1:1203835/1‑150 (MQ=255)
 cGAACGCGGAAAATCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTTGCACTGGCacc                                                                                                                                             <  2:287168/150‑3 (MQ=255)
 cGAACGCGGAAAATCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGCACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTTGCACTGGCacc                                                                                                                                             <  1:623276/150‑3 (MQ=255)
  gAACGCGGAAAATCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTTGCACTGGCaccc                                                                                                                                            <  1:1489165/150‑4 (MQ=255)
  gAACGCGGAAAATCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGt                                                                                                                                            >  1:487809/1‑150 (MQ=255)
      gcgGAAAATCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGc                                                                                                                                        >  2:744190/1‑150 (MQ=255)
         gAAAATCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCaaa                                                                                                                                     >  1:985845/1‑150 (MQ=255)
          aaaaTCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAAc                                                                                                                                    <  1:313936/150‑1 (MQ=255)
            aaTCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAg                                                                                                                                  <  1:434438/150‑1 (MQ=255)
                  cAGTAACAGCAGACCTGTCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCcat                                                                                                                            >  2:957235/1‑150 (MQ=255)
                   aGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCcatc                                                                                                                           >  1:1087699/1‑150 (MQ=255)
                   aGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCcatc                                                                                                                           <  2:985845/150‑1 (MQ=255)
                     tAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCcatcat                                                                                                                         <  1:20744/150‑1 (MQ=255)
                      aaCAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCata                                                                                                                        <  1:223818/150‑1 (MQ=255)
                        cagcagACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCataat                                                                                                                      >  2:424601/1‑150 (MQ=255)
                        cagcagACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCataat                                                                                                                      >  1:473247/1‑150 (MQ=255)
                         agcagACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCataata                                                                                                                     >  2:1133965/1‑150 (MQ=255)
                          gcagACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATAc                                                                                                                    <  2:191742/150‑1 (MQ=255)
                               ccTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCg                                                                                                               >  1:1012984/1‑150 (MQ=255)
                               ccTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCg                                                                                                               >  2:228191/1‑150 (MQ=255)
                               ccTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCg                                                                                                               >  2:1418205/1‑150 (MQ=255)
                                 tATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGaa                                                                                                             >  2:951188/1‑150 (MQ=255)
                                  aTCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCa                                                                                                                                               <  1:677278/115‑1 (MQ=255)
                                  aTCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCa                                                                                                                                               >  2:677278/1‑115 (MQ=255)
                                  aTCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAAt                                                                                                            <  1:682154/150‑1 (MQ=255)
                                   tCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCATTGCAACTGGCAGGTTTGCAAAGAGTCCCATCATAATACTGCCGAATg                                                                                                           <  2:849866/150‑1 (MQ=255)
                                    cgcgCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCAATCGCCTCTACAACGACAAAAGCGAAGAACGCATACAGGCCCATAGGCTGTGCCAGTGAAACTGGCAGGATAGCAAACAGTCCCAGCATAATACTGCCGAAGGc                                                                                                          >  2:776311/1‑150 (MQ=255)
                                      cgcCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAg                                                                                                        >  1:1469165/1‑150 (MQ=255)
                                      cgcCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAg                                                                                                        <  2:1315390/150‑1 (MQ=255)
                                      cgcCCAAGAAGAATGCGCCCATCCCGACCTGCCACGGCACGCCCATCGCCTGTACAACGACAACAGCGAAGAACGCAATCCGGCCCATAGCGGCTGCCAGTGCAACTGGCAGCTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAg                                                                                                        >  2:697461/1‑150 (MQ=255)
                                        ccccAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCg                                                                                                      <  2:247022/150‑1 (MQ=255)
                                        ccccAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCg                                                                                                      <  2:252546/150‑1 (MQ=255)
                                        ccccAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCg                                                                                                      >  2:1088717/1‑150 (MQ=255)
                                           cAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATc                                                                                                   <  2:1012984/150‑1 (MQ=255)
                                             taagaTTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATTGCGGGTGCCAGTGCAACTGGCAGGGTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAg                                                                                                 <  2:831982/149‑1 (MQ=255)
                                                gaTTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACa                                                                                              <  1:661831/150‑1 (MQ=255)
                                                    gcgcCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTa                                                                                          >  2:537031/1‑150 (MQ=255)
                                                         cATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGAc                                                                                     >  2:252150/1‑150 (MQ=255)
                                                              cGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGa                                                                                >  1:1381052/1‑150 (MQ=255)
                                                               gACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGAc                                                                               <  2:1161156/150‑1 (MQ=255)
                                                                 ccTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACgg                                                                             >  1:303982/1‑150 (MQ=255)
                                                                    gCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGgcgc                                                                          >  1:965121/1‑150 (MQ=255)
                                                                       aCGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTgg                                                                       <  2:555317/150‑1 (MQ=255)
                                                                          gCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTAt                                                                    <  2:663540/150‑1 (MQ=255)
                                                                               cccATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATg                                                               <  1:635492/150‑1 (MQ=255)
                                                                                 cATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGcc                                                             >  1:1366545/1‑150 (MQ=255)
                                                                                   tCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAg                                                           <  2:487809/150‑1 (MQ=255)
                                                                                      ccTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCGACTGGCAGGGTAGCGTACAGTCCCATCATAATACTGCCGGCTGCAGCGAGCCGACAGGTAGTGACGAGGGCGGCGCTGGTGGCCAGGCCCGccc                                                        >  2:1428222/1‑148 (MQ=255)
                                                                                        tGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACg                                                      >  1:1106477/1‑150 (MQ=255)
                                                                                          tACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGcc                                                    <  1:997252/150‑1 (MQ=255)
                                                                                           aCAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCa                                                   >  2:586090/1‑150 (MQ=255)
                                                                                                gACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAAt                                              >  2:163862/1‑150 (MQ=255)
                                                                                                 aCAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAAtt                                             >  2:434642/1‑150 (MQ=255)
                                                                                                   aaaaGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCACGTTAGCAAACAGTCCCATCATAAAACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTg                                           >  1:699777/1‑150 (MQ=255)
                                                                                                     aaGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCg                                         <  1:252150/150‑1 (MQ=255)
                                                                                                       gCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCggg                                       <  2:1469165/150‑1 (MQ=255)
                                                                                                       gCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCggg                                       >  2:27407/1‑150 (MQ=255)
                                                                                                          aagaaCGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTa                                    <  1:1088717/150‑1 (MQ=255)
                                                                                                          aagaaCGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTa                                    >  2:743909/1‑150 (MQ=255)
                                                                                                                 cATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACaaaaa                             >  1:841488/1‑150 (MQ=255)
                                                                                                                     cAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGAt                         >  2:896163/1‑150 (MQ=255)
                                                                                                                      aGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATg                        <  1:1133965/150‑1 (MQ=255)
                                                                                                                       ggCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGt                       <  1:744190/150‑1 (MQ=255)
                                                                                                                          ccATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACTAAAACGATGTaaa                    <  1:361258/150‑1 (MQ=255)
                                                                                                                             tAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCa                 >  1:1435248/1‑150 (MQ=255)
                                                                                                                              ctaacctGCCAGTTGCACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCAt                >  1:1318056/7‑150 (MQ=255)
                                                                                                                               gCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATc               >  2:895739/1‑150 (MQ=255)
                                                                                                                                 gggTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGt             >  1:501035/1‑150 (MQ=255)
                                                                                                                                 acctGCCAGTTGCACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGt             >  2:110690/4‑150 (MQ=255)
                                                                                                                                 acctGCCAGTTGCACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGt             >  1:757207/4‑150 (MQ=255)
                                                                                                                                  cctGCCAGTTGCACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTc            >  2:58515/3‑150 (MQ=255)
                                                                                                                                  cctGCCAGTTGCACTGGCAGGTTAGCAAACAGTCCCATCATAATACTCCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTGAACAAAAACGATGTAAACCATCGTc            >  2:115768/3‑150 (MQ=255)
                                                                                                                                    tGCCAGTTGCACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAg          >  1:1291801/1‑150 (MQ=255)
                                                                                                                                           gCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAggtg   >  2:1002210/1‑150 (MQ=255)
                                                                                                                                           gCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAggtg   >  1:549671/1‑150 (MQ=255)
                                                                                                                                            cAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAggtgg  >  1:514874/1‑150 (MQ=255)
                                                                                                                                             |                                                                                                                                                    
GCGAACGCGGAAAATCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAGAATTTGCGGGTTAACAAAAACGATGTAAACCATCGTCAGGAAGGTGG  >  NC_000913/3896211‑3896500

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