Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,147,586 A→C 7.5% intergenic (+46/‑35) rpmF → / → plsX 50S ribosomal subunit protein L32/putative phosphate acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,147,5860AC 7.5% 296.5 / 9.7 99intergenic (+46/‑35)rpmF/plsX50S ribosomal subunit protein L32/putative phosphate acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (51/40);  new base C (6/2);  total (57/42)
Fisher's exact test for biased strand distribution p-value = 4.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCACCAGCCTGTCTGTAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTT  >  NC_000913/1147440‑1147733
                                                                                                                                                  |                                                                                                                                                   
tCACCAGCCTGTCTGTAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCACCTg                                                                                                                                                  >  1:1654371/1‑150 (MQ=255)
 cACCAGCCTGTCTGTAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTgg                                                                                                                                                 <  1:41931/150‑1 (MQ=255)
 cACCAGCCTGTCTGTAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTgg                                                                                                                                                 <  1:42137/150‑1 (MQ=255)
 cACCAGCCTGTCTGTAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTgg                                                                                                                                                 <  1:42140/150‑1 (MQ=255)
  aCCAGCCTGTCTGTAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTggg                                                                                                                                                >  2:1251948/1‑150 (MQ=255)
  aCCAGCCTGTCTGTAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTggg                                                                                                                                                >  2:1479123/1‑150 (MQ=255)
  aCCAGCCTGTCTGTAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTggg                                                                                                                                                >  2:386816/1‑150 (MQ=255)
    cAGCCTGTCTGTAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGa                                                                                                                                              >  1:1117391/1‑150 (MQ=255)
    cAGCCTGTCTGTAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGa                                                                                                                                              >  1:1132865/1‑150 (MQ=255)
    cAGCCTGTCTGTAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGa                                                                                                                                              >  2:1401272/1‑150 (MQ=255)
     aGCCTGTCTGTAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGTTGCCTGGGGaa                                                                                                                                             <  2:167259/150‑1 (MQ=255)
     aGCCTGTCTGTAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGTTGCCTGGGGaa                                                                                                                                             <  1:939904/150‑1 (MQ=255)
     aGCCTGTCTGTAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGaa                                                                                                                                             >  2:946688/1‑150 (MQ=255)
     aGCCTGTCTGTAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGaa                                                                                                                                             >  2:1499234/1‑150 (MQ=255)
     aGCCTGTCTGTAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGaa                                                                                                                                             >  2:946685/1‑150 (MQ=255)
         tgtctgtAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGAc                                                                                                                                         <  2:1160265/150‑1 (MQ=255)
         tgtctgtAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGAc                                                                                                                                         <  2:919991/150‑1 (MQ=255)
                 gACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCgg                                                                                                                                 <  2:389916/150‑1 (MQ=255)
                  aCAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCggg                                                                                                                                >  2:409871/1‑150 (MQ=255)
                    aaaaCTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGggcg                                                                                                                              >  1:1636589/1‑150 (MQ=255)
                    aaaaCTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGggcg                                                                                                                              >  2:1561682/1‑150 (MQ=255)
                     aaaCTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAATCCGggcgg                                                                                                                             >  2:1152758/1‑150 (MQ=255)
                     aaaCTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGggcgg                                                                                                                             >  2:1152754/1‑150 (MQ=255)
                      aaCTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGggcggc                                                                                                                            <  2:311114/150‑1 (MQ=255)
                         ttCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGcgac                                                                                                                         >  1:1345264/1‑150 (MQ=255)
                            tGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGAt                                                                                                                      >  1:939538/1‑150 (MQ=255)
                              gagAAAAACAACTGAGACACCACAACACTCCCCACGGTTACTACCGCCGCCGCAAGGTCATCCCTAAGAAATCACGAATCTGCGGGATGAAGCTTAGTGAGGATTTTCCCAAGTCACCTGGGGAAAGACAAAACCGGGCGGCGACGATAc                                                                                                                    >  2:1542089/3‑150 (MQ=255)
                              gTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATAc                                                                                                                    >  2:1542097/1‑150 (MQ=255)
                               tGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAATTAATCACGCATCTGCGTGATGTAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATAcc                                                                                                                   <  1:909002/150‑1 (MQ=255)
                               tGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATAcc                                                                                                                   >  1:219765/1‑150 (MQ=255)
                                 aaaaaCACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCtt                                                                                                                 >  2:1377935/1‑150 (MQ=255)
                                       acCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACg                                                                                                           >  1:1099851/1‑150 (MQ=255)
                                        ccTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGt                                                                                                          <  2:1651727/150‑1 (MQ=255)
                                         cTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTc                                                                                                         >  1:666815/1‑150 (MQ=255)
                                           gCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTa                                                                                                       >  1:729810/1‑150 (MQ=255)
                                           gCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTa                                                                                                       >  1:1585008/1‑150 (MQ=255)
                                                   acaTCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGc                                                                                               >  1:523313/1‑150 (MQ=255)
                                                     aTCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGt                                                                                             >  2:189953/1‑150 (MQ=255)
                                                     aTCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGt                                                                                             >  2:703237/1‑150 (MQ=255)
                                                       cACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTACGTACTCTCGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGCCCAAACCGGGCGGCGACGATACCTTGTCACGTCTCACCCTGGCGTTa                                                                                           <  2:1685654/150‑1 (MQ=255)
                                                       cACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTa                                                                                           <  1:1154223/150‑1 (MQ=255)
                                                         cTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGa                                                                                         <  1:1667675/150‑1 (MQ=255)
                                                          tGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGAt                                                                                        <  2:1335922/150‑1 (MQ=255)
                                                          tGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGAt                                                                                        <  2:1317728/150‑1 (MQ=255)
                                                          tGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGAt                                                                                        <  2:1275509/150‑1 (MQ=255)
                                                          tGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGAt                                                                                        <  2:83159/150‑1 (MQ=255)
                                                          tGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGAt                                                                                        <  1:199935/150‑1 (MQ=255)
                                                          tGCCGACGGGTACTCCCGCGGCCGCACGGTCATCGCTAAGTCAGCACGCAACTGCGTGCTGAAGCTTAGTGAGGCTTGTCCCCAGGCAACTGGGGAAAGACCAACCCGGGCGGCGACGATACCTGGACACGTCTAACCCGGGCGTTAGAt                                                                                        <  2:948656/150‑1 (MQ=255)
                                                                  gTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATggg                                                                                >  2:1444156/1‑150 (MQ=255)
                                                                  gTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGATTAGTGAGGATTTTCCCCAGGCAACTGGTGAAAGACCAAAACGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATggg                                                                                >  2:196972/1‑150 (MQ=255)
                                                                   ttACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGa                                                                               >  1:1083755/1‑150 (MQ=255)
                                                                       tacCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAgggg                                                                           >  1:1526528/1‑150 (MQ=255)
                                                                         ccGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGAt                                                                         >  1:1213249/1‑150 (MQ=255)
                                                                         ccGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGAt                                                                         >  1:1568861/1‑150 (MQ=255)
                                                                               ccGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGc                                                                   <  2:8101/150‑1 (MQ=255)
                                                                               ccGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGc                                                                   <  2:21738/150‑1 (MQ=255)
                                                                               ccGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGc                                                                   <  2:792579/150‑1 (MQ=255)
                                                                                cGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTTACCCTGGCGTTTGATGTCATGGGAGGGGATTTTGGcc                                                                  <  2:368870/150‑1 (MQ=255)
                                                                                   aaGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCtt                                                               >  1:759732/1‑150 (MQ=255)
                                                                                    aGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTc                                                              <  2:1634888/150‑1 (MQ=255)
                                                                                     ggTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTcc                                                             <  1:795110/150‑1 (MQ=255)
                                                                                      gTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCg                                                            >  1:24558/1‑150 (MQ=255)
                                                                                        cATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTAGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTg                                                          >  2:506062/1‑150 (MQ=255)
                                                                                              tAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTg                                                    <  1:485293/150‑1 (MQ=255)
                                                                                                    aTCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTgca                                              <  2:746980/150‑1 (MQ=255)
                                                                                                    aTCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTgca                                              <  2:938347/150‑1 (MQ=255)
                                                                                                    aTCACGCAGATGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTgca                                              >  2:582332/1‑150 (MQ=255)
                                                                                                       aCGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTgcagca                                           >  1:985748/1‑150 (MQ=255)
                                                                                                       aCGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTgcagca                                           >  2:596622/1‑150 (MQ=255)
                                                                                                           aTCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGc                                       <  1:563153/150‑1 (MQ=255)
                                                                                                                cGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGTGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCa                                  <  1:409871/150‑1 (MQ=255)
                                                                                                                cGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCa                                  <  2:1099851/150‑1 (MQ=255)
                                                                                                                 gTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCAc                                 <  2:1117391/150‑1 (MQ=255)
                                                                                                                 gTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCAc                                 >  2:935456/1‑150 (MQ=255)
                                                                                                                     tgaAGCTTAGTGAGGATGTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAATGCTGCAGGCACTgtc                             <  2:1010199/150‑3 (MQ=255)
                                                                                                                      gaAGCTTAGTGAGGTTTTTCCCCAGGCAACGGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGTCACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCGCTTCCGTGACTGTGCCTGCAGCATTGCAGGCACTGAAt                            <  1:492488/150‑1 (MQ=255)
                                                                                                                       aaGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAAtt                           >  2:1352538/1‑150 (MQ=255)
                                                                                                                       aaGCATAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAAtt                           >  2:1593290/1‑150 (MQ=255)
                                                                                                                       aaGCATAGTGAGGATTTTCCACAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAAtt                           >  2:1593296/1‑150 (MQ=255)
                                                                                                                           ttAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTaa                       <  1:1338248/150‑1 (MQ=255)
                                                                                                                               tGAGGATTTTCCCCAGGCACCTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCGTTGCAGGCACTGAATTCTAATTCg                   >  1:1044252/1‑150 (MQ=255)
                                                                                                                               tGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCg                   >  1:257196/1‑150 (MQ=255)
                                                                                                                               tGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCg                   >  2:352015/1‑150 (MQ=255)
                                                                                                                               tGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCg                   >  2:1415132/1‑150 (MQ=255)
                                                                                                                               tGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCg                   >  1:767026/1‑150 (MQ=255)
                                                                                                                                gAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGc                  >  2:194490/1‑150 (MQ=255)
                                                                                                                                gAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGc                  >  1:1527496/1‑150 (MQ=255)
                                                                                                                                gAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGc                  >  2:170100/1‑150 (MQ=255)
                                                                                                                                  ggATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCaa                <  1:15245/150‑1 (MQ=255)
                                                                                                                                  ggATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCaa                <  1:643630/150‑1 (MQ=255)
                                                                                                                                  ggATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCaa                <  1:861650/150‑1 (MQ=255)
                                                                                                                                    tctttCCCCAGTTGCCTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAact              >  1:1013699/3‑150 (MQ=255)
                                                                                                                                    tctttCCCCAGTTGCCTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCAGGCGTTAGATGTCATGGGAGGGGAGTTTGGCCCTTCCGGGACAGTGCCTGCAGCAGTGCAGGCACTGAATTCTAATTAGCAact              >  1:107739/3‑150 (MQ=255)
                                                                                                                                    aTTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAact              <  2:666815/150‑1 (MQ=255)
                                                                                                                                      tttCCCCAGTTGCCTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAactca            >  2:686109/1‑150 (MQ=255)
                                                                                                                                      tttCCCCAGGCAACTGGGGGAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAactca            >  2:1349281/1‑150 (MQ=255)
                                                                                                                                        tCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAactcact          >  2:733568/1‑150 (MQ=255)
                                                                                                                                                cAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCtttt  <  1:1716690/150‑1 (MQ=255)
                                                                                                                                                cAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCtttt  <  1:1716732/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TCACCAGCCTGTCTGTAGACAAAACTTCTGGTGAAAAACACCTGCGTCACCACATCACTGCCGACGGTTACTACCGCGGCCGCAAGGTCATCGCTAAGTAATCACGCATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTT  >  NC_000913/1147440‑1147733

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