Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,885,300 T→G 5.4% R75R (CGT→CGG yidC → membrane protein insertase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,885,3000TG 5.4% 318.7 / 4.9 115R75R (CGT→CGGyidCmembrane protein insertase
Reads supporting (aligned to +/- strand):  ref base T (46/62);  new base G (4/3);  total (50/65)
Fisher's exact test for biased strand distribution p-value = 4.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 128 total aligned reads displayed.

ATAAAAACCCGCAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGATAACCCGGCT  >  NC_000913/3885152‑3885447
                                                                                                                                                    |                                                                                                                                                   
aTAAAAACCCGCAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACCCCCgtg                                                                                                                                                    >  2:670978/1‑150 (MQ=255)
 tAAAAACCCGCAACCTCAGGCCCAACAGACCACGCAGACAACGACCCCCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCgtgg                                                                                                                                                   >  1:1372383/1‑150 (MQ=255)
  aaaaaCCCGCAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCgtggt                                                                                                                                                  <  2:1149462/150‑1 (MQ=255)
  aaaaaCCCGCAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCgtggt                                                                                                                                                  <  1:1366141/150‑1 (MQ=255)
   aaaaCCCGCAACCTCAGGCCCAACAGACCACGCAGCCAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGATACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCgtggtg                                                                                                                                                 >  1:727635/1‑150 (MQ=255)
    acaCCCGCAACCTCAGGCCCAACAGACCCCGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCgtggtgg                                                                                                                                                <  2:1347742/148‑1 (MQ=255)
    aaaCCCGCAACCTCAGGCCCAACAGACCACGCAGACAACGCCCACCGCCGCGGGTAGCGCCGCCGACCAGGGCGTTCCGGCCAGGGGCCAGGGGAAACTGATCGCGGGTAAGAACGCCGTGCTTGATCTGCCCATCAACACCCgtggtgg                                                                                                                                                >  1:157050/1‑150 (MQ=255)
    aaaCCCGCAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCgtggtgg                                                                                                                                                <  2:1249976/150‑1 (MQ=255)
    aaaCCCGCAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCgtggtgg                                                                                                                                                <  1:879672/150‑1 (MQ=255)
     aaCCCGCAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGTGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCgtggtggt                                                                                                                                               <  2:1075554/150‑1 (MQ=255)
     aaCCCGCAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCgtggtggt                                                                                                                                               <  1:195586/150‑1 (MQ=255)
     aaCCCGCAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCgtggtggt                                                                                                                                               >  1:792091/1‑150 (MQ=255)
      aCCCGCAACCTCAGGCCCAACACACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGCGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCgtggtggtg                                                                                                                                              <  2:181489/150‑1 (MQ=255)
       cccGCAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGa                                                                                                                                             <  1:1776384/150‑1 (MQ=255)
         cGCAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGCTCTGACCATCAACACCCGTGGTGGTGATg                                                                                                                                           >  1:1673768/1‑150 (MQ=255)
           cAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGtt                                                                                                                                         >  2:502472/1‑150 (MQ=255)
                  aGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGCTCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGcaag                                                                                                                                  <  1:1281527‑M2/150‑5 (MQ=255)
                     cccAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGAGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGctc                                                                                                                               >  1:1612976/1‑150 (MQ=255)
                       cAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTctg                                                                                                                             <  2:11028/150‑1 (MQ=255)
                            gACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGctgcc                                                                                                                        <  1:1667988/150‑1 (MQ=255)
                              ccACGCAGCCAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGctgcctg                                                                                                                      <  2:1461730/150‑1 (MQ=255)
                                     gACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGCCGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCGTAccc                                                                                                               >  1:543753/1‑150 (MQ=255)
                                      aCAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTCCCCg                                                                                                              >  1:211241/1‑150 (MQ=255)
                                       cAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGCCCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGa                                                                                                             >  1:764295/1‑150 (MQ=255)
                                       cAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGa                                                                                                             >  1:467485/1‑150 (MQ=255)
                                       cAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGa                                                                                                             >  2:826211/1‑150 (MQ=255)
                                       cAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGAGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGa                                                                                                             <  1:673765/150‑1 (MQ=255)
                                          cGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAag                                                                                                          >  1:162354/1‑150 (MQ=255)
                                          cGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAag                                                                                                          <  2:1264297/150‑1 (MQ=255)
                                            accaccGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGCCCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAagag                                                                                                        >  2:550545/1‑150 (MQ=255)
                                            accaccGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAagag                                                                                                        <  1:1200163/150‑1 (MQ=255)
                                             ccaccGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGc                                                                                                       >  2:405526/1‑150 (MQ=255)
                                                  gCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCCAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAAc                                                                                                  >  1:1840217/1‑150 (MQ=255)
                                                  gCAGCGGGTAGCGCCGCCGACCAGGGCGAACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTAAAGACCGACGTGCTTGATCTTACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAAc                                                                                                  <  2:378399/150‑1 (MQ=255)
                                                     gCGGGTATCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGGTACGACCGACGGGCGGGAGCTGCCCATCAACACCCGGGGGGGGGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAActct                                                                                               <  1:973580/150‑1 (MQ=255)
                                                      cGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTa                                                                                              <  2:1832560/150‑1 (MQ=255)
                                                       gggTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTAc                                                                                             >  2:592784/1‑150 (MQ=255)
                                                        ggTAGCGCCGCCGACCAGGGCGTACCGGACAGTGGCCAGGGGAAAGTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTAcc                                                                                            <  2:64841/150‑1 (MQ=255)
                                                          tAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGTGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTTATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTAccct                                                                                          <  2:736569/150‑2 (MQ=255)
                                                            gcgcCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGc                                                                                        <  1:623041/150‑1 (MQ=255)
                                                                cgccgACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGtt                                                                                    >  2:437284/1‑150 (MQ=255)
                                                                  ccgACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTcc                                                                                  >  2:1534256/1‑150 (MQ=255)
                                                                      ccAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCt                                                                              >  1:12417/1‑150 (MQ=255)
                                                                      ccAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCt                                                                              <  1:385247/150‑1 (MQ=255)
                                                                      ccAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCt                                                                              >  1:299603/1‑150 (MQ=255)
                                                                         gggCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGtt                                                                           <  2:1523272/150‑1 (MQ=255)
                                                                         gggCGTACCGGCCAGTGGCCAGGGGAAACTGATCACGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGtt                                                                           <  2:564776/150‑1 (MQ=255)
                                                                          ggCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTg                                                                          <  1:1802597/150‑1 (MQ=255)
                                                                            cGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGa                                                                        <  2:1339620/150‑1 (MQ=255)
                                                                              tACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGaaa                                                                      >  2:815907/1‑150 (MQ=255)
                                                                              tACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGaaa                                                                      <  2:464199/150‑1 (MQ=255)
                                                                                  ggCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCCGCCGTTCCAGCTGTTGGAAACTTc                                                                  >  2:1404441/1‑150 (MQ=255)
                                                                                  ggCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTc                                                                  <  2:878937/150‑1 (MQ=255)
                                                                                    ccAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCAc                                                                >  1:1844882/1‑150 (MQ=255)
                                                                                    ccAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCAc                                                                <  2:1327708/150‑1 (MQ=255)
                                                                                     cAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCAcc                                                               <  2:1815220/150‑1 (MQ=255)
                                                                                     cAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCAcc                                                               <  1:65550/150‑1 (MQ=255)
                                                                                      aGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCg                                                              <  1:1831650/150‑1 (MQ=255)
                                                                                      aGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCg                                                              <  1:59095/150‑1 (MQ=255)
                                                                                          gCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGt                                                          >  2:1539213/1‑150 (MQ=255)
                                                                                           ccAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGtt                                                         >  1:1371294/1‑150 (MQ=255)
                                                                                                 gAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGtttattta                                                   <  1:1801287/150‑1 (MQ=255)
                                                                                                  aaaCTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTttatttat                                                  >  1:698110/1‑150 (MQ=255)
                                                                                                    aCTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCa                                                <  1:960429/150‑1 (MQ=255)
                                                                                                     cTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAg                                               <  1:1480937/150‑1 (MQ=255)
                                                                                                      tGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAgg                                              >  1:753327/1‑150 (MQ=255)
                                                                                                       gATCTCGGTTAAGACCGACGAGCGTGAGCTGACCATCAACACCCGTGGGGGGGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGc                                             <  1:1724381/150‑1 (MQ=255)
                                                                                                           tcGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACag                                         >  1:1319513/1‑150 (MQ=255)
                                                                                                            cGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACaga                                        <  2:752579/150‑1 (MQ=255)
                                                                                                              gTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGc                                      <  2:362507/150‑1 (MQ=255)
                                                                                                                tAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCgg                                    <  2:848619/150‑1 (MQ=255)
                                                                                                                   taCCGACGTGCTTTATCTGACCATCAACACCCGTGTTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAATAGCTGAACTCTACCCAGCCGATCCAGCTGTTGGAAACTTCACCTCAGTTTATTTATCAGGCACAGAGCGGTCt                                 <  1:220616/149‑1 (MQ=255)
                                                                                                                     ccGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGcaagcacgtcggtcttaaccgagatcagtttcccctggccactggccggtacgccctggtcggcggcgctacccgctgcggtggtcgttgtctgcgtggtctg                               <  1:771334‑M2/150‑104 (MQ=255)
                                                                                                                      cGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGAc                              <  1:189583/150‑1 (MQ=255)
                                                                                                                       gACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGcaagcacgtcggtcttaaccgagatcagtttcccctggccactggccggtacgccctggtcggcggcgctacccgctgcggtgggcgttgtctgcgtggtctgtt                             <  2:1688076‑M2/150‑106 (MQ=255)
                                                                                                                        aCGTGCTTGATCTGACCATCAACACCCGTGGTTGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCg                            <  1:1278118/150‑1 (MQ=255)
                                                                                                                          gTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGt                          <  2:162354/150‑1 (MQ=255)
                                                                                                                           tGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAAATTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTc                         >  2:1553122/1‑150 (MQ=255)
                                                                                                                            gCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCg                        <  2:764295/150‑1 (MQ=255)
                                                                                                                                gATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGAt                    >  2:1734442/1‑150 (MQ=255)
                                                                                                                                  tCTGACCATGTACACCCGTGGTGGTGATGTTGAGCAAGCGCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATgg                  <  2:251214/150‑1 (MQ=255)
                                                                                                                                  gctcaacATCACCACCACGGGTGTTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATgg                  >  1:680158/7‑150 (MQ=255)
                                                                                                                                   ctcaacATCACCACCACGGGTGTTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGc                 >  1:1801964/6‑150 (MQ=255)
                                                                                                                                   ctcaacATCACCACCACGGGTGTTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGc                 >  1:693260/6‑150 (MQ=255)
                                                                                                                                    tGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGcc                >  2:840773/1‑150 (MQ=255)
                                                                                                                                     gACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGccc               <  1:1115739/150‑1 (MQ=255)
                                                                                                                                       ccATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCTCAGAGCGGTCTGACCGGTCGTGATGGCCCgg             >  1:502169/1‑150 (MQ=255)
                                                                                                                                       ccATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGCGCGGTCTGACCGGTCGTGATGGCCCgg             >  2:218921/1‑150 (MQ=255)
                                                                                                                                       ccATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCgg             >  1:1137385/1‑150 (MQ=255)
                                                                                                                                       ccATCAACACCCGTGGTGGTGATGGTCAGatcaagcacgtcggtcttaaccgagatcagtttcccctggccactggccggtacgccctggtcggcggcgctacccgctgcgggggtcggtgtgtgcgtggtctggtgggcctgaggttgg             <  2:739316‑M2/150‑122 (MQ=255)
                                                                                                                                        cATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTAAACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGAGGGCCCggt            >  2:1442454/1‑149 (MQ=255)
                                                                                                                                         aTCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGAt           >  2:211694/1‑150 (MQ=255)
                                                                                                                                         aTCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGAt           >  1:233946/1‑150 (MQ=255)
                                                                                                                                           cAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGATaa         >  2:1813094/1‑150 (MQ=255)
                                                                                                                                            accacCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTCCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTGGTGATGGCCCGGATAAc        <  2:532649/148‑1 (MQ=255)
                                                                                                                                            aaCACCCGTGGTGGTGATGTTGAGCAATCTCTGCTGCCTTCTTACCCGTAAGAGGTTAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGATAAc        <  1:1056652/150‑1 (MQ=255)
                                                                                                                                            aaCACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGATAAc        >  2:948744/1‑150 (MQ=255)
                                                                                                                                             acacCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGATAAcc       >  2:777823/1‑150 (MQ=255)
                                                                                                                                              cacCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGATAAccc      <  2:875008/150‑1 (MQ=255)
                                                                                                                                               acCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGATAACCCg     <  2:600263/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ATAAAAACCCGCAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGATAACCCGGCT  >  NC_000913/3885152‑3885447

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