Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,969,428 T→C A519A (GCA→GCG tap ← methyl‑accepting protein IV

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,969,4280TC99.2% 441.0 / NA 129A519A (GCA→GCGtapmethyl‑accepting protein IV
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (54/74);  minor base G (0/1);  total (54/75)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 129 total aligned reads displayed.

CCGAAAATGCGCGTCGGAAAGCGCCAGGCGCTCGGTCATCTGTAACAATAAAGACGTTTGCCCACAGGGCAGAGATGAAGTCATAGCGCCTTCTCAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGTGCTACTTCATGTTCTTCAAGGGTAAATACCGCCACACGCGACGAAAGATGGTCGGCCTGATTCGCCAGTTGTTCCGTTGCCACCGCCGCTTCTTCTACCAGCGAGGCGTTCTGCTGAGTCACCTGATCCATCTGGCTGACAGCCT  >  NC_000913/1969281‑1969573
                                                                                                                                                   |                                                                                                                                                 
ccgAAAATGCGCGTCGGAAAGCGCCAGGCGCTCGGTCATCTGTAACAATAAAGACGTTTGCCCACAGGGCAGAGATGAAGTCATAGCGCCTTCTCAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCt                                                                                                                                                >  1:1330904/1‑151 (MQ=255)
  gAAAATGCGCGTCGGAAAGCGCCAGGCGCTCGGTCATCTGTAACAATAAAGACGTTTGCCCACAGGGCAGAGATGAAGTCATAGCGCCTTCTCAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTAc                                                                                                                                              >  1:872201/1‑151 (MQ=255)
  gAAAATGCGCGTCGGAAAGCGCCAGGCGCTCGGTCATCTGTAACAATAAAGACGTTTGCCCACAGGGCAGAGATGAAGACATAGCGCCTTCTCAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTAc                                                                                                                                              >  2:578403/1‑151 (MQ=255)
   aaaaTGCGCGTCGGAAAGCGCCAGGCGCTCGGTCATCTGTAACAATAAAGACGTTTGCCCACAGGGCAGAGATGAAGTCATAGCGCCTTCTCAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACt                                                                                                                                             >  1:1964682/1‑151 (MQ=255)
   aaaaTGCGCGTCGGAAAGCGCCAGGCGCTCGGTCATCTGTAACAATAAAGACGTTTGCCCACAGGGCAGAGATGAAGTCATAGCGCCTTCTCAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACt                                                                                                                                             >  1:447977/1‑151 (MQ=255)
    aaaTGCGCGTCGGAAAGCGCCAGGCGCTCGGTCATCTGTAACAATAAAGACGTTTGCCCACAGGGCAGAGATGAAGTCATAGCGCCTTCTCAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACtt                                                                                                                                            >  2:1665494/1‑151 (MQ=255)
     aaTGCGCGTCGGAAAGCGCCAGGCGCTCGGTCATCTGTAACAATAAAGACGTTTGCCCACAGGGCAGAGATGAAGTCATAGCGCCTTCTCAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACTTc                                                                                                                                           <  2:2068960/151‑1 (MQ=255)
     aaTGCGCGTCGGAAAGCGCCAGGCGCTCGGTCATCTGTAACAATAAAGACGTTTGCCCACAGGGCAGAGATGAAGTCATAGCGCCTTCTCAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACTTc                                                                                                                                           >  2:2112249/1‑151 (MQ=255)
     aaTGCGCGTCGGAAAGCGCCAGGCGCTCGGTCATCTGTAACAATAAAGACGTTTGCCCACAGGGCAGAGATGAAGTCATAGCGCCTTCTCAATCACTTCAGGATACCACTGGCGCAATTCGTAACTGCACCGACTCATGTCGCGCTACTTc                                                                                                                                           <  2:1171699/151‑1 (MQ=255)
      aTGCGCGTCGGAAAGCGCCAGGCGCTCGGTCATCTGTAACAATAAAGACGTTTGCCCACAGGGCAGAGATGAAGTCATAGCGCCTTCTCAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACTTc                                                                                                                                           >  2:2144193/1‑150 (MQ=255)
      aTGCGCGTCGGAAAGCGCCAGGCGCTCGGTCATCTGTAACAATAAAGACGTTTGCCCACAGGGCAGAGATGAAGTCATAGCGCCTTCTCAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACTTCa                                                                                                                                          <  2:10287/151‑1 (MQ=255)
        gcgcGTCGGAAAGCGCCAGGCGCTCGGTCATCTGTAACAATAAAGACGTTTGCCCACAGGGCAGAGATGAAGTCATAGCGCCTTCTCAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACTTCATg                                                                                                                                        >  1:1314827/1‑151 (MQ=255)
         cgcgTCGGAAAGCGCCAGGCGCTCGGTCATCTGTAACAATAAAGACGTTTGCCCACAGGGCAGAGATGAAGTCATAGCGCCTTCTCAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACTTCATGt                                                                                                                                       >  1:784669/1‑151 (MQ=255)
          gcgTCGGAAAGCGCCAGGCGCTCGGTCATCTGTAACAATAAAGACGTTTGCCCACAGGGCAGAGATGAAGTCATAGCGCCTTCTCAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACTTCATGtt                                                                                                                                      <  2:1829001/151‑1 (MQ=255)
           cgTCGGAAAGCGCCAGGCGCTCGGTCATCTGTAACAATAAAGACGTTTGCCCACAGGGCAGAGATGAAGTCATAGCGCCTTCTCAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACTTCATGttc                                                                                                                                     <  1:697970/151‑1 (MQ=255)
              cGGAAAGCGCCAGGCGCTCGGTCATCTGTAACAATAAAGACGTTTGCCCACAGGGCAGAGATGAAGTCATAGCGCCTTCTCAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACTTCATGttctt                                                                                                                                   <  2:1582719/150‑1 (MQ=255)
               ggAAAGCGCCAGGCGCTCGGTCATCTGTAACAATAAAGACGTTTGCCCACAGGGCAGAGATGAAGTCATAGCGCCTTCTCAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACTTCATGttcttc                                                                                                                                  >  2:2042189/1‑150 (MQ=255)
               ggAAAGCGCCAGGCGCTCGGTCATCTGTAACAATAAAGACGTTTGCCCACAGGGCAGAGATGAAGTCATAGCGCCTTCTCAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACTTCATGTTCTTCa                                                                                                                                 <  1:2103901/151‑1 (MQ=255)
                gAAAGCGCCAGGCGCTCGGTCATCTGTAACAATAAAGACGTTTGCCCACAGGGCAGAGATGAAGTCATAGCGCCTTCTCAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACTTCATGTTCTTCaa                                                                                                                                <  2:913496/151‑1 (MQ=255)
                gAAAGCGCCAGGCGCTCGGTCATCTGTAACAATAAAGACGTTTGCCCACAGGGCAGAGATGAAGTCATAGCGCCTTCTCAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACTTCATGTTCTTCaa                                                                                                                                <  2:764074/151‑1 (MQ=255)
                 aaaGCGCCAGGCGCTCGGTCATCTGTAACAATAAAGACGTTTGCCCACAGGGCAGAGATGAAGTCATAGCGCCTTCTCAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACTTCATGTTCTTCAAg                                                                                                                               <  1:2158662/151‑1 (MQ=255)
                  aaGCGCCAGGCGCTCGGTCATCTGTAACAATAAAGACGTTTGCCCACAGGGCAGAGATGAAGTCATAGCGCCTTCTCAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACTTCATGTTCTTCAAgg                                                                                                                              <  1:1453883/151‑1 (MQ=255)
                  aaGCGCCAGGCGCTCGGTCATCTGTAACAATAAAGACGTTTGCCCACAGGGCAGAGATGAAGTCATAGCGCCTTCTCAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACTTCATGTTCTTCAAgg                                                                                                                              >  1:1267444/1‑151 (MQ=255)
                   aGCGCCAGGCGCTCGGTCATCTGTAACAATAAAGACGTTTGCCCACAGGGCAGAGATGAAGTCATAGCGCCTTCTCAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACTTCATGTTCTTCAAggg                                                                                                                             <  1:203550/151‑1 (MQ=255)
                   aGCGCCAGGCGCTCGGTCATCTGTAACAATAAAGACGTTTGCCCACAGGGCAGAGATGAAGTCATAGCGCCTTCTCAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACTTCATGTTCTTCAAggg                                                                                                                             >  2:102111/1‑151 (MQ=255)
                       ccAGGCGCTCGGTCATCTGTAACAATAAAGACGTTTGCCCACAGGGCAGAGATGAAGTCATAGCGCCTTCTCAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACTTCATGTTCTTCAAGGGTaaa                                                                                                                         <  1:1540092/151‑1 (MQ=255)
                         aGGCGCTCGGTCATCTGTAACAATAAAGACGTTTGCCCACAGGGCAGAGATGAAGTCATAGCGCCTTCTCAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACTTCATGTTCTTCAAGGGTAAATa                                                                                                                       <  1:864060/151‑1 (MQ=255)
                          ggCGCTCGGTCATCTGTAACAATAAAGACGTTTGCCCACAGGGCAGAGATGAAGTCATAGCGCCTTCTCAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACTTCATGTTCTTCAAGGGTAAATAc                                                                                                                      <  2:475237/151‑1 (MQ=255)
                           gcgcTCGGTCATCTGTAACAATAAAGACGTTTGCCCACAGGGCAGAGATGAAGTCATAGCGCCTTCTCAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACTTCATGTTCTTCAAGGGTAAATAcc                                                                                                                     <  2:1168851/151‑1 (MQ=255)
                             gcTCGGTCATCTGTAACAATAAAGACGTTTGCCCACAGGGCAGAGATGAAGTCATAGCGCCTTCTCAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACTTCATGTTCTTCAAGGGTAAATACCGc                                                                                                                   >  2:1801635/1‑151 (MQ=255)
                                 ggTCATCTGTAAAAAAAAAGACTTTTGCCCACGGGGCAGAGATAAATTCATAGCCCCTTCTCAATCACTTAAGTATCCCCCTGCCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACTTCATTTTCTTCAAGGGTAAATACCGCcaca                                                                                                               <  2:2087613/151‑1 (MQ=255)
                                    cATCTGTAACAATAAAGACGTTTGCCCACAGGGCAGAGATGAAGTCATAGCGCCTTCTCAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACTTCATGTTCTTCAAGGGTAAATACCGCCACAcgc                                                                                                            >  2:1034976/1‑151 (MQ=255)
                                      tCTGTAACAATAAAGACGTTTGCCCACAGGGCAGAGATGAAGTCATAGCGCCTTCTCAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACTTCATGTTCTTCAAGGGTAAATACCGCCACAcgcg                                                                                                           >  1:1298584/1‑150 (MQ=255)
                                         gTAACAATAAAGACGTTTGCCCACAGGGCAGAGATGAAGTCATAGCGCCTTCTCAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACTTCATGTTCTTCAAGGGTAAATACCGCCACACGcgac                                                                                                         <  2:1162985/149‑1 (MQ=255)
                                         gTAACAATAAAGACGTTTGCCCACAGGGCAGAGATGAAGTCATAGCGCCTTCTCAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACTTCATGTTCTTCAAGGGTAAATACCGCCACACGcga                                                                                                          <  2:872545/148‑1 (MQ=255)
                                         gTAACAATAAAGACGTTTGCCCACAGGGCAGAGATGAAGTCATAGCGCCTTCTCAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACTTCATGTTCTTCAAGGGTAAATACCGCCACACGcga                                                                                                          <  1:1682680/148‑1 (MQ=255)
                                           aaCAATAAAGACGTTTGCCCACAGGGCAGAGATGAAGTCATAGCGCCTTCTCACTCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACTTCATGTTCTTCAAGGGTAAATACCGCCACACGCGACGaa                                                                                                      <  2:1232996/150‑1 (MQ=255)
                                           aaCAATAAAGACGTTTGCCCACAGGGCAGAGATGAAGTCATAGCGCCTTCTCAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACTTCATGTTCTTCAAGGGTAAATACCGCCACACGCGACGaaa                                                                                                     <  1:2118295/151‑1 (MQ=255)
                                           aaCAATAAAGACGTTTGCCCACAGGGCAGAGATGAAGTCATAGCGCCTTCTCAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACTTCATGTTCTTCAAGGGTAAATACCGCCACACGCGACGaa                                                                                                      >  1:648622/1‑150 (MQ=255)
                                             cAATAAAGACGTTTGCCCACAGGGCAGAGATGAAGTCATAGCGCCTTCTCAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACTTCATGTTCTTCAAGGGTAAATACCGCCACACGCGACGAAAGa                                                                                                   <  2:161119/151‑1 (MQ=255)
                                              aaTAAAGACGTTTGCCCACAGGGCAGAGATGAAGTCATAGCGCCTTCTCAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACTTCATGTTCTTCAAGGGTAAATACCGCCACACGCGACGAAAGa                                                                                                   <  2:1414869/150‑1 (MQ=255)
                                              aaTAAAGACGTTTGCCCACAGGGCAGAGATGAAGTCATAGCGCCTTCTCAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACTTCATGTTCTTCAAGGGTAAATACCGCCACACGCGACGAAAGAt                                                                                                  <  2:231381/151‑1 (MQ=255)
                                               aTAAAGACGTTTGCCCACAGGGCAGAGATGAAGTCATAGCGCCTTCTCAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACTTCATGTTCTTCAAGGGTAAATACCGCCACACGCGACGAAAGATg                                                                                                 >  1:721480/1‑151 (MQ=255)
                                                  aaGACGTTTGCCCACAGGGCAGAGATGAAGTCATAGCGCCTTCTCAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACTTCATGTTCTTCAAGGGTAAATACCGCCACACGCGACGAAAGATGGTc                                                                                              <  1:701323/151‑1 (MQ=255)
                                                  aaGACGTTTGCCCACAGGGCAGAGATGAAGTCATAGCGCCTTCTCAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACTTCATGTTCTTCAAGGGTAAATACCGCCACACGCGACGAAAGATGGTc                                                                                              <  1:166606/151‑1 (MQ=255)
                                                   aGACGTTTGCCCACAGGGCAGAGATGAAGTCATAGCGCCTTCTCAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACTTCATGTTCTTCAAGGGTAAATACCGCCACACGCGACGAAAGATGGTCg                                                                                             >  1:2006897/1‑151 (MQ=255)
                                                   aGACGTTTGCCCACAGGGCAGAGATGAAGTCATAGCGCCTTCTCAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACTTCATGTTCTTCAAGGGTAAATACCGCCACACGCGACGAAAGATGGTCg                                                                                             <  2:841637/151‑1 (MQ=255)
                                                      cGTTTGCCCACAGGGCAGAGATGAAGTCATAGCGCCTTCTCAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACTTCATGTTCTTCAAGGGTAAATACCGCCACACGCGACGAAAGATGGTCGGcc                                                                                          >  2:262702/1‑151 (MQ=255)
                                                       gTTTGCCCACAGGGCAGAGATGAAGTCATAGCGCCTTCTCAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACTTCATGTTCTTCAAGGGTAAATACCGCCACACGCGACGAAAGATGGTCGGCCt                                                                                         >  2:1993990/1‑151 (MQ=255)
                                                         ttGCCCACAGGGCAGAGATGAAGTCATAGCGCCTTCTCAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACTTCATGTTCTTCAAGGGTAAATACCGCCACACGCGACGAAAGATGGTCGGCCTg                                                                                        <  2:1034193/150‑1 (MQ=255)
                                                            cccACAGGGCAGAGATGAAGTCATAGCGCCTTCTCAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACTTCATGTTCTTCAAGGGTAAATACCGCCACACGCGACGAAAGATGGTCGGCCTGATTc                                                                                    >  2:1020561/1‑151 (MQ=255)
                                                                caGGGCAGAGATGAAGTCATAGCGCCTTCTCAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACTTCATGTTCTTCAAGGGTAAATACCGCCACACGCGACGAAAGATGGTCGGCCTGATTCGCCa                                                                                <  1:1216926/151‑1 (MQ=255)
                                                                  gggCAGAGATGAAGTCATAGCGCCTTCTCAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACTTCATGTTCTTCAAGGGTAAATACCGCCACACGCGACGAAAGATGGTCGGCCTGATTCGCCAgt                                                                              <  2:2134465/151‑1 (MQ=255)
                                                                   ggCAGAGATGAAGTCATAGCGCCTTCTCAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACTTCATGTTCTTCAAGGGTAAATACCGCCACACGCGACGAAAGATGGTCGGCCTGATTCGCCAgtt                                                                             >  1:505862/1‑151 (MQ=255)
                                                                   ggCAGAGATGAAGTCATAGCGCCTTCTCAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACCTCATGTTCTTCAAGGGTAAATACCGCCACACGCGACGAAAGATGGTCGGCCTGATTCGCCAgtt                                                                             <  2:46433/151‑1 (MQ=255)
                                                                    gCAGAGATGAAGTCATAGCGCCTTCTCAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACTTCATGTTCTTCAAGGGTAAATACCGCCACACGCGACGAAAGATGGTCGGCCTGATTCGCCAgttg                                                                            >  2:224352/1‑151 (MQ=255)
                                                                      agagATGAAGTCATAGCGCCTTCTCAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACTTCATGTTCTTCAAGGGTAAATACCGCCACACGCGACGAAAGATGGTCGGCCTGATTCGCCAgttgt                                                                           <  2:1659533/150‑1 (MQ=255)
                                                                         gaTGAAGTCATAGCGCCTTCTCAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACTTCATGTTCTTCAAGGGTAAATACCGCCACACGCGACGAAAGATGGTCGGCCTGATTCGCCAGTTGTTCCg                                                                       >  1:2010515/1‑151 (MQ=255)
                                                                           tGAAGTCATAGCGCCTTCTCAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACTTCATGTTCTTCAAGGGTAAATACCGCCACACGCGACGAAAGATGGTCGGCCTGATTCGCCAGTTGTTcc                                                                        <  1:1117875/148‑1 (MQ=255)
                                                                            gAAGTCATAGCGCCTTCTCAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACTTCATGTTCTTCAAGGGTAAATACCGCCACACGCGACGAAAGATGGTCGGCCTGATTCGCCAGTTGTTCCGTTg                                                                    <  1:910438/151‑1 (MQ=255)
                                                                            gAAGTCATAGCGCCTTCTCAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACTTCATGTTCTTCAAGGGTAAATACCGCCACACGCGACGAAAGATGGTCGGCCTGATTCGCCAGTTGTTCCGTTg                                                                    <  2:1528305/151‑1 (MQ=255)
                                                                             aaGTCATAGCGCCTTCTCAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACTTCATGTTCTTCAAGGGTAAATACCGCCACACGCGACGAAAGATGGTCGGCCTGATTCGCCAGTTGTTCCGTTGc                                                                   <  1:1575077/151‑1 (MQ=255)
                                                                               gTCATAGCGCCTTCTCAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACTTCATGTTCTTCAAGGGTAAATACCGCCACACGCGACGAAAGATGGTCGGCCTGATTCGCCAGTTGTTCCGTTGCCa                                                                 <  1:1684234/151‑1 (MQ=255)
                                                                                  aTAGCGCCTTCTCAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACTTCATGTTCTTCAAGGGTAAATACCGCCACACGCGACGAAAGATGGTCGGCCTGATTCGCCAGTTGTTCCGTTGCCAcc                                                               <  1:1145069/150‑1 (MQ=255)
                                                                                      cgcCTTCTCAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCCCCGACTCATGTCGCGCTACTTCATTTTCTTAAAGGGTAAATACCGCCACACGCGACGAAAGATGGTCGGCCTGATTCGCCAGTTGTTCCGTTGCCAccgc                                                             <  2:1263818/148‑1 (MQ=255)
                                                                                      cgcCTTCTCAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATTTCGCGCTACTTCATGTTCTTCAAGGGTAAATACCGCCACACGCGACGAAAGATGGTCGGCCTGATTCGCCAGTTGTTCCGTTGCCAccgcc                                                            <  1:1905438/149‑1 (MQ=255)
                                                                                       gcCTTCTCAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACTTCATGTTCTTCAAGGGTAAATACCGCCACACGCGACGAAAGATGGTCGGCCTGATTCGCCAGTTGTTCCGTTGCCACCGCCGct                                                         >  2:1157668/1‑151 (MQ=255)
                                                                                           tctcAATCCCTTCAGGATACCACTGGCGCATTTTGTAACTGCCCCGACTCATGTCGCGCTCCTTCATGTTCTTCAAGGGTAAATCCCGCCACACGCGACGAAAGATGGTCGGCCTGATTCGCCAGTTGTTCCGTTGCCACCGCCGcttct                                                      <  1:854046/150‑1 (MQ=255)
                                                                                            ctcAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACTTCATGTTCTTCAAGGGTAAATACCGCCACACGCGACGAAAGATGGTCGGCCTGATTCGCCAGTTGTTCCGTTGCCACCGCCGCttcttc                                                    >  2:1428740/1‑151 (MQ=255)
                                                                                            ctcAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACTTCATGTTCTTCAAGGGTAAATACCGCCACACGCGACGAAAGATGGTCGGCCTGATTCGCCAGTTGTTCCGTTGCCACCGCCGCttcttc                                                    >  2:1151719/1‑151 (MQ=255)
                                                                                            ctcAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACTTCATGTTCTTCAAGGGTAAATACCGCCACACGCGACGAAAGATGGTCGGCCTGATTCGCCAGTTGTTCCGTTGCCACCGCCGCttcttc                                                    >  2:1694173/1‑151 (MQ=255)
                                                                                            ctcAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACTTCATGTTCTTCAAGGGTAAATACCGCCACACGCGACGAAAGATGGTCGGCCTGATTCGCCAGTTGTTCCGTTGCCACCGCCGCttcttc                                                    >  1:56353/1‑151 (MQ=255)
                                                                                             tcAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACTTCATGTTCTTCAAGGGTAAATACCGCCACACGCGACGAAAGATGGTCGGCCTGATTCGCCAGTTGTTCCGTTGCCACCGCCGCTtcttct                                                   >  1:806156/1‑151 (MQ=255)
                                                                                              cAATCCCTTCAGGATCCCCCTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTCCTTCATGTTCTTAAAGGGTAAATACCGCCACACGCGACGAAAGATGGTCGGCCTGATTCGCCAGTTGTTCCGTTGCCACCGCCGcttctt                                                     <  1:1151556/148‑1 (MQ=255)
                                                                                               aaTCCCTTCAGGATACCACGGGCGCAATTTGTAACTGCCCCGACTCATGTCGCGCTACTTCATGTTCTTAAAGGGTAAATCCCGCCACACGCGACGAAAGATGGTCGGCCTGATTCGCCAGTTGTTCCGTTGCCACCGCCGcttc                                                       <  1:133487/145‑1 (MQ=255)
                                                                                               aaTCACTTCGGAATCCCACTGGCGCATTTTGTAACTGCACCGACTCATGTCGGGCTACTTCATGTTCTTCAAGGGTAAATCCCGCCACACGCGACGAAAGATGGTCGGCCTGATTCGCCAGTTGTTCCGTTGCCACCGCCGCttcttc                                                    <  1:369965/148‑1 (MQ=255)
                                                                                               aaTCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACTTCATGTTCTTCAAGGGTAAATACCGCCACACGCGACGAAAGATGGTCGGCCTGATTCGCCAGTTGTTCCGTTGCCACCGCCGCTTCTTCTAc                                                 >  2:1998618/1‑151 (MQ=255)
                                                                                                  cACTTCAGGATACCCCTGGCGCATTTTGTAACTGCCCCGACTCATGTCGCGCTACTTCATTTTCTTCAAGGTTAAATCCCGCCACACGCGACGAAAGATGGTCGGCCTGATTCGCCAGTTGTTCCGTTGCCACCGCCGCTTCTTCTACCa                                               <  2:594855/150‑1 (MQ=255)
                                                                                                       cAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACTTCATGTTCTTCAAGGGTAAATACCGCCACACGCGACGAAAGATGGTCGGCCTGATTCGCCAGTTGTTCCGTTGCCACCGCCGCTTCTTCTACCAGCGAgg                                         >  1:1882739/1‑151 (MQ=255)
                                                                                                             aCCCCTGGCCAAATTTGTACCTGCACCGACCCATGTCGCGCTTCTCCATTTTCTTCAAGGGTAAAACCCGCCACCCGCGACGAAAGATGGTCGGCCTGATTCGCCAGTTGTTCCGTTGCCACCGCCGCTTCTTCTACCAGCGAGGCGtt                                     <  2:130818/149‑1 (MQ=255)
                                                                                                             aCCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACTTCATGTTCTTCAGGGGTAAATACCGCCACACGCGACGAAAGATGGTCGGCCTGATTCGCCAGTTGTTCCGTTGCCACCGCCGCTTCTTCTACCAGCGAGGCGTTct                                   <  1:409413/151‑1 (MQ=255)
                                                                                                              ccccTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACTTCATGTTCTTCAAGGGTAAATACCGCCACACGCGACGAAAGATGGTCGGCCTGATTCGCCAGTTGTTCCGTTGCCACCGCCGCTTCTTCTACCAGCGAGGCGTTct                                   <  1:1545784/147‑1 (MQ=255)
                                                                                                              ccACTGGCGCAATTTGTAACTGCCCCGACTCATGTCGCGCTACTTCATGTTCTTCAAGGGTAAATCCCGCCACACGCGACGAAAGATGGTCGGCCTGATTCGCCAGTTGTTCCGTTGCCAccgcc                                                            <  2:870365/125‑1 (MQ=255)
                                                                                                              ccACTGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACTTCATGTTCTTCAAGGGTAAATACCGCCACACGCGACGAAAGATGGTCGGCCTGATTCGCCAGTTGTTCCGTTGCCAccgc                                                             >  1:870377/1‑124 (MQ=255)
                                                                                                                  tGGCGCAATTTGTAACTGCACCGACTCATGTCGCGCTACTTCATGTTCTTCAAGGGTAAATACCGCCACACGCGACGAAAGATGGTCGGCCTGATTCGCCAGTTGTTCCGTTGCCACCGCCGCTTCTTCTACCAGCGAGGCGTTctgctg                               >  1:832525/1‑150 (MQ=255)
                                                                                                                     cgaaaTTTGTAACTGCCCCGACTCATGTCGCGCTACTTCATGTTCTTCAAGGGTAAATACCGCCACACGCGACGAAAGATGGTCGGCCTGATTCGCCAGTTGTTCCGTTGCCACCGCCGCTTCTTCTACCAGCGAGGCGTTCTGCTGAGt                            <  1:1852543/147‑1 (MQ=255)
                                                                                                                       cAATTTGTAACTGCACCGACTCATGTCGCGCTACTTCATGTTCTTCAAGGGTAAATACCGCCACACGCGACGAAAGATGGTCGGCCTGATTCGCCAGTTGTTCCGTTGCCACCGCCGCTTCTTCTACCAGCGAGGCGTTCTGCTGAGTCAc                         >  1:147548/1‑151 (MQ=255)
                                                                                                                         aTTTGTAACTGCACCGACTCATGTCGCGCTACTTCATGTTCTTCAAGGGTAAATACCGCCACACGCGACGAAAGATGGTCGGCCTGATTCGCCAGTTGTTCCGTTGCCACCGCCGCTTCTTCTACCAGCGAGGCGTTCTGCTGAGTCACCt                       >  1:303159/1‑151 (MQ=255)
                                                                                                                          tttGTAACTGCACCGACTCATGTCGCGCTACTTCATGTTCTTCAAGGGTAAATACCGCCACACGCGACGAAAGATGGTCGGCCTGATTCGCCAGTTGTTCCGTTGCCACCGCCGCTTCTTCTACCAGCGAGGCGTTCTGCTGAGTCACCTg                      >  2:818738/1‑151 (MQ=255)
                                                                                                                            tGTAACTGCACCGACTCATGTCGCGCTACTTCATGTTCTTCAAGGGTAAATACCGCCACACGCGACGAAAGATGGTCGGCCTGATTCGCCAGTTGTTCCGTTGCCACCGCCGCTTCTTCTACCAGCGAGGCGTTCTGCTGAGTCACCTGAt                    >  1:1525742/1‑151 (MQ=255)
                                                                                                                            tGTAACTGCACCGACTCATGTCGCGCTACTTCATGTTCTTCAAGGGTAAATACCGCCACACGCGACGAAAGATGGTCGGCCTGATTCGCCAGTTGTTCCGTTGCCACCGCCGCTTCTTCTACCAGCGAGGCGTTCTGCTGAGTCACCTGAt                    >  2:58986/1‑151 (MQ=255)
                                                                                                                             gTAACTGCACCGACTCATGTCGCGCTACTTCATGTTCTTCAAGGGTAAATACCGCCACACGCGACGAAAGATGGTCGGCCTGATTCGCCAGTTGTTCCGTTGCCACCGCCGCTTCTTCTACCAGCGAGGCGTTCTGCTGAGTCACCTGATc                   >  1:388195/1‑151 (MQ=255)
                                                                                                                              taccTGCACCGCCTCATGTCGCGCTACTTCATTTTCTTCAAGGGTAAATACCGCCACACGCGACGAAAGATGGTCGGCCTGATTCGCCAGTTGTTCCGTTGCCACCGCCGCTTCTTCTACCAGCGAGGCGTTCTGCTGAGTCACCTGAt                    <  1:1147455/146‑1 (MQ=255)
                                                                                                                               aaCTGCACCGACTCATGTCGCGCTACTTCATGTTCTTCAAGGGTAAATACCGCCACACGCGACGAAAGATGGTCGGCCTGATTCGCCAGTTGTTCCGTTGCCACCGCCGCTTCTTCTACCAGCGAGGCGTTCTGCTGAGTCACCTGATcc                  >  2:65985/1‑150 (MQ=255)
                                                                                                                               aaCAGCACCGACTCATGTCGCGCTACTTCATGTTCTTCAAGGGTAAATACCGCCACACGCGACGAAAGATGGTCGGCCTGATTCGCCAGTTGTTCCGTTGCCACCGCCGCTTCTTCTACCAGCGAGGCGTTCTGCTGAGTCACCTGATcc                  >  1:1222674/1‑150 (MQ=255)
                                                                                                                                            cATGTCGCGCTACTTCATGTTCTTCAAGGGTAAATACCGCCACACGCGACGAAAGATGGTCGGCCTGATTCGCCAGTTGTTCCGTTGCCACCGCCGCTTCTTCTACCAGCGAGGCGTTCTGCTGAGTCACCTGATCCATCTGGCTGACAGc    >  2:1930647/1‑151 (MQ=255)
                                                                                                                                             aTGTCGCGCTACTTCATGTTCTTCAAGGGTAAATACCGCCACACGCGACGAAAGATGGTCGGCCTGATTCGCCAGTTGTTCCGTTGCCACCGCCGCTTCTTCTACCAGCGAGGCGTTCTGCTGAGTCACCTGATCCATCTGGCTGACAGcc   >  2:1029554/1‑151 (MQ=255)
                                                                                                                                             aTGCCGCGCTACTTCATGTTTTTAAAGGGAAAATCCCGCCACACGCGACGAAAGATGGTCGGCCTGATTCGCCAGTTGTTCCGTTGCCACCGCCGCTTCTTCTACCAGCGAGGCGTTCTGCTGAGTCACCTGATCCATCTGGCTGa        <  2:472223/146‑1 (MQ=255)
                                                                                                                                              tGTCGCGCTACTTCATGTTCTTCAAGGGTAAATCCCGCCACACGCGACGAAAGATGGTCGGCCTGATTCGCCAGTTGTTCCGTTGCCACCGCCGCTTCTTCTACCAGCGAGGCGTTCTGCTGAGTCACCTGATCCATCTGGCTGACAGcc   <  2:240384/150‑1 (MQ=255)
                                                                                                                                               gTCGCGCTACTTCATTTTCTTCAAGGGTAAATACCGCCACACGCGACGAAAGATGGTCGGCCTGATTCGCCAGTTGTTCCGTTGCCACCGCCGCTTCTTCTACCAGCGAGGCGTTCTGCTGAGTCACCTGATCCATCTGGCTGACAGCCt  <  1:1371803/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CCGAAAATGCGCGTCGGAAAGCGCCAGGCGCTCGGTCATCTGTAACAATAAAGACGTTTGCCCACAGGGCAGAGATGAAGTCATAGCGCCTTCTCAATCACTTCAGGATACCACTGGCGCAATTTGTAACTGCACCGACTCATGTCGTGCTACTTCATGTTCTTCAAGGGTAAATACCGCCACACGCGACGAAAGATGGTCGGCCTGATTCGCCAGTTGTTCCGTTGCCACCGCCGCTTCTTCTACCAGCGAGGCGTTCTGCTGAGTCACCTGATCCATCTGGCTGACAGCCT  >  NC_000913/1969281‑1969573

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