Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,201,688 A→C T64T (ACA→ACC ymfI → e14 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,201,6880AC98.5% 258.5 / ‑2.7 68coding (192/342 nt)ymfIe14 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base C (31/36);  minor base . (0/1);  total (31/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.64e-01

GAACAGACCAAAAAATTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACATTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCTTACTTCCGTCTGCTTCCTTAGTTGGTAACGATTTCGTTCATTACGATCTTACGCCTAAGCCCATCAAGCTCGATTCT  >  NC_000913/1201542‑1201835
                                                                                                                                                  |                                                                                                                                                   
gAACAGACCAAAAAATTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCtttt                                                                                                                                                 >  1:910851/1‑151 (MQ=255)
 aaCAGACCAAAAAATTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTc                                                                                                                                                <  2:257862/151‑1 (MQ=255)
   cAGACCAAAAAATTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCa                                                                                                                                               >  2:800517/1‑150 (MQ=255)
          aaaaaTTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTaa                                                                                                                                       <  2:883583/151‑1 (MQ=255)
            aaaTTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAAc                                                                                                                                      >  1:518814/1‑150 (MQ=255)
              aTTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGt                                                                                                                                   <  1:349278/151‑1 (MQ=255)
                ttGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTg                                                                                                                                  <  1:755846/150‑1 (MQ=255)
                ttGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGc                                                                                                                                 >  2:174959/1‑151 (MQ=255)
                 tGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGcc                                                                                                                                >  1:282424/1‑151 (MQ=255)
                  gaagaaTTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCg                                                                                                                               <  2:511155/151‑1 (MQ=255)
                       aTTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGtt                                                                                                                          <  2:84002/151‑1 (MQ=255)
                         tAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTg                                                                                                                         <  2:912658/150‑1 (MQ=255)
                          aaTGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGAc                                                                                                                       >  1:447494/1‑151 (MQ=255)
                           aTGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACg                                                                                                                      >  2:98884/1‑151 (MQ=255)
                           aTGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACg                                                                                                                      >  1:630304/1‑151 (MQ=255)
                                tAACTTTCTGTACCAGAAAACAGTTGAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCaaa                                                                                                                  <  2:740271/150‑1 (MQ=255)
                                tAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCaaa                                                                                                                  <  2:10571/150‑1 (MQ=255)
                                tAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCaa                                                                                                                   <  2:16326/149‑1 (MQ=255)
                                tAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAat                                                                                                                 >  1:66220/1‑151 (MQ=255)
                                 aaCTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAat                                                                                                                 <  1:754136/150‑1 (MQ=255)
                                     ttCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAatatat                                                                                                             <  1:193276/150‑1 (MQ=255)
                                      tCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATAtt                                                                                                            <  1:504451/150‑1 (MQ=255)
                                         gTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTg                                                                                                          <  1:199373/149‑1 (MQ=255)
                                           aCCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTg                                                                                                      <  2:511551/151‑1 (MQ=255)
                                              aGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTaa                                                                                                   <  2:229970/151‑1 (MQ=255)
                                              aGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTa                                                                                                    >  2:554305/1‑150 (MQ=255)
                                                  aaCAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAAttt                                                                                               >  2:382211/1‑151 (MQ=255)
                                                   aCAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTg                                                                                              <  1:109920/151‑1 (MQ=255)
                                                   aCAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTg                                                                                              <  2:896199/151‑1 (MQ=255)
                                                     aGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGaa                                                                                            <  1:194390/151‑1 (MQ=255)
                                                     aGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGaa                                                                                            <  2:552098/151‑1 (MQ=255)
                                                       ttAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAAc                                                                                           >  2:788832/1‑150 (MQ=255)
                                                        tAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACa                                                                                          >  1:175309/1‑150 (MQ=255)
                                                           agagTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAgc                                                                                      <  1:333587/151‑1 (MQ=255)
                                                            gagTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAgcg                                                                                     <  1:689630/151‑1 (MQ=255)
                                                             agTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAgcgc                                                                                    >  2:267520/1‑151 (MQ=255)
                                                                  gACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATg                                                                                 <  1:465542/149‑1 (MQ=255)
                                                                   aCGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCtt                                                                              <  2:165941/151‑1 (MQ=255)
                                                                        aGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCTTACTTc                                                                         >  2:545428/1‑151 (MQ=255)
                                                                         gctgctGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCTTACTTcc                                                                        >  2:727642/1‑151 (MQ=255)
                                                                         gctgctGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCTTACTTcc                                                                        <  1:879053/151‑1 (MQ=255)
                                                                         gctgctGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCTTACTTcc                                                                        <  1:793332/151‑1 (MQ=255)
                                                                            gctgcAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCTTACTTCCGTCt                                                                    <  1:860181/151‑1 (MQ=255)
                                                                              tgcAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCTTACTTCCGTCTg                                                                   >  2:840832/1‑151 (MQ=255)
                                                                                          ccAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCTTACTTCCGTCTGCTTCCTTAGTTg                                                       <  2:545075/151‑1 (MQ=255)
                                                                                          ccAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCTTACTTCCGTCTGCTTCCTTAGTTg                                                       >  1:155777/1‑151 (MQ=255)
                                                                                                ggATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCTTACTTCCGTCTGCTTCCTTAGTTGGTAACg                                                 >  1:140084/1‑151 (MQ=255)
                                                                                                 gATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCTTACTTCCGTCTGCTTCCTTAGTTGGTAACg                                                 >  1:63180/1‑150 (MQ=255)
                                                                                                  atatatCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCTTACTTCCATCTGCTTCCTTAGTTGGTAACGAt                                               >  1:874979/1‑151 (MQ=255)
                                                                                                         cTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCTTACTTCCGTCTGCTTCCTTAGTTGGTAACGATTTCGTTc                                        >  2:668103/1‑151 (MQ=255)
                                                                                                             aaaGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCTTACTTCCGTCTGCTTCCTTAGTTGGTAACGATTTCGTTCATTa                                    <  1:214612/151‑1 (MQ=255)
                                                                                                             aaaGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCTTACTTCCGTCTGCTTCCTTAGTTGGTAACGATTTCGTTCATTa                                    <  1:346248/151‑1 (MQ=255)
                                                                                                               aGTACAATGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCTTACTTCCGTCTGCTTCCTTAGTTGGTAACGATTTCGTTCATTACg                                  <  2:81613/151‑1 (MQ=255)
                                                                                                                     aTGAACTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCTTACTTCCGTCTGCTTCCTTAGTTGGTAACGATTTCGTTCATTACGATCTTa                            <  2:792242/151‑1 (MQ=255)
                                                                                                                        aaCTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCTTACTTCCGTCTGCTTCCTTAGTTGGTAACGATTTCGTTCATTACGATCTTACGc                         >  1:225824/1‑151 (MQ=255)
                                                                                                                          cTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCTTACTTCCGTCTGCTTCCTTAGTTGGTAACGATTTCGTTCATTACGATCTTACGCCt                       >  2:305923/1‑151 (MQ=255)
                                                                                                                          cTGCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCTTACTTCCGTCTGCTTCCTTAGTTGGTAACGATTTCGTTCATTACGATCTTACGCCt                       >  2:255257/1‑151 (MQ=255)
                                                                                                                            gCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCTTACTTCCGTCTGCTTCCTTAGTTGGTAACGATTTCGTTCATTACGATCTTACGCCTaa                     <  2:252531/151‑1 (MQ=255)
                                                                                                                            gCCAACAAATCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCTTACTTCCGTCTGCTTCCTTAGTTGGTAACGATTTCGTTCATTACGATCTTACGCCTaa                     <  1:499370/151‑1 (MQ=255)
                                                                                                                                  aaaTCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCTTACTTCCGTCTGCTTCCTTAGTTGGTAACGATTTCGTTCATTACGATCTTACGCCTAAGCCCAt               >  1:96376/1‑151 (MQ=255)
                                                                                                                                   aaTCCTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCTTACTTCCGTCTGCTTCCTTAGTTGGTAACGATTTCGTTCATTACGATCTTACGCCTAAGCCCATc              <  2:543616/151‑1 (MQ=255)
                                                                                                                                       cTTGATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCTTACTTCCGTCTGCTTCCTTAGTTGGTAACGATTTCGTTCATTACGATCTTACGCCTAAGCCCATCAAg           >  1:572552/1‑150 (MQ=255)
                                                                                                                                          gATCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCTTACTTCCGTCTGCTTCCTTAGTTGGTAACGATTTCGTTCATTACGATCTTACGCCTAAGCCCATCAAGCTc        <  1:382219/150‑1 (MQ=255)
                                                                                                                                            tCAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCTTACTTCCGTCTGCTTCCTTAGTTGGTAACGATTTCGTTCATTACGATCTTACGCCTAAGCCCATCAAGCTCGAt     >  2:762232/1‑151 (MQ=255)
                                                                                                                                             cAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCTTACTTCCGTCTGCTTCCTTAGTTGGTAACGATTTCGTTCATTACGATCTTACGCCTAAGCCCATCAAGCTCGAtt    >  2:384257/1‑151 (MQ=255)
                                                                                                                                             cAAACCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCTTACTTCCGTCTGCTTCCTTAGTTGGTAACGATTTCGTTCATTACGATCTTACGCCTAAGCCCATCAAGCTCGAt     <  1:806905/150‑1 (MQ=255)
                                                                                                                                              aaaCCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCTTACTTCCGTCTGCTTCCTTAGTTGGTAACGATTTCGTTCATTACGATCTTACGCCTAAGCCCATCAAGCTCGATTc   >  1:841667/1‑151 (MQ=255)
                                                                                                                                               aaCCTTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCTTACTTCCGTCTGCTTCCTTAGTTGGTAACGATTTCGTTCATTACGATCTTACGCCTAAGCCCATCAAGCTCGATTCt  >  1:766524/1‑151 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GAACAGACCAAAAAATTTGAAGAATTAATGCTTAACTTTCTGTACCAGAAAACAGTTAAAGAGTCTGACGATAGCTGCTGCAGACTGATTCCAGAGGGATATATCCTCAAAAGTACAATGAACTGCCAACAAATCCTTGATCAAACATTTTCAATTGCTAACAGTGCCGGTGTTGACGCAAATATATTTGTCTGTAAATTTGAACAAAGCGCATGCTTACTTCCGTCTGCTTCCTTAGTTGGTAACGATTTCGTTCATTACGATCTTACGCCTAAGCCCATCAAGCTCGATTCT  >  NC_000913/1201542‑1201835

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: