Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,483,900 T→A 7.6% V49D (GTC→GAC)  evgA → response regulator in two‑component regulatory system with EvgS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,483,9000TA 7.6% 245.9 / 7.3 91V49D (GTC→GAC) evgAresponse regulator in two‑component regulatory system with EvgS
Reads supporting (aligned to +/- strand):  ref base T (35/48);  new base A (2/5);  total (38/53)
Fisher's exact test for biased strand distribution p-value = 6.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

ATCTATGAACGCAATAATTATTGATGACCATCCTCTTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACGGGAAACATTGTGCTGATGCTGGCGCTAATGGTT  >  NC_000913/2483751‑2484049
                                                                                                                                                     |                                                                                                                                                     
atCTATGAACGCAATAATTATTGATGACCATCCTCTTGCTATCGCAGCAATTCGTACTTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGt                                                                                                                                                       <  1:214145/150‑1 (MQ=255)
atCTATGAACGCAATAATTATTGATGACCATCCTCTTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGt                                                                                                                                                       >  2:134685/1‑150 (MQ=255)
  atATGAACGCAATAATTATTGATGACCATCCTCTTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGtca                                                                                                                                                     <  1:554272/149‑1 (MQ=255)
    aTGAACGCAATAATTATTGATGACCATCCTCTTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGtcatc                                                                                                                                                   <  2:1190761/150‑1 (MQ=255)
     tGAACGCAATAATTATTGATGACCATCCTCTTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGtcatca                                                                                                                                                  <  1:373385/150‑1 (MQ=255)
       aaCGCAATAATTATTGATGACCATCCTCTTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCAtt                                                                                                                                                >  1:1194049/1‑150 (MQ=255)
        aCGCAATAATTATTGATGACCATCCTCTTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTg                                                                                                                                               <  2:207364/150‑1 (MQ=255)
         cGCAATAATTATTGATGCCCATCCTCTTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGa                                                                                                                                              <  1:724515/150‑1 (MQ=255)
         cGCAATAATTATTGATGACCATCCTCTTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGa                                                                                                                                              >  1:1649723/1‑150 (MQ=255)
          gCAATAATTATTGATGACCATCCTCTTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGAt                                                                                                                                             >  2:817083/1‑150 (MQ=255)
             ataatTATTGATGACCATCCTATTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTc                                                                                                                                          <  2:315250/150‑1 (MQ=255)
             ataatTATTGATGACCATACTCTTTCTATCGAAGCAATTCTTAATTTATTGATCAAAAACGATAATGAAATATTAGCAGAGTTGACTGAAGGCTGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCATTATATCGTCATCATTGATTTc                                                                                                                                          <  2:180439/150‑1 (MQ=255)
               aatTATTGATGACCATCCTCTTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGa                                                                                                                                        <  2:1019318/150‑1 (MQ=255)
                 ttattGATGACCATCCTCTTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGata                                                                                                                                      <  1:1331321/150‑1 (MQ=255)
                     tgatgaCCATCCTCTTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGACATCAATGATGACGATATcag                                                                                                                                  <  1:50019/150‑3 (MQ=255)
                      gatgaCCATCCTCTTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATcccc                                                                                                                                 >  2:526613/1‑150 (MQ=255)
                      gatgaCCATCCTCTTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGACATCAATGATGACGATATcagg                                                                                                                                 <  1:139573/150‑4 (MQ=255)
                      gatgaCCATCCTCTTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGACATCAATGATGACGATATcagg                                                                                                                                 <  2:633147/150‑4 (MQ=255)
                       atgaCCATCCTCTTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCg                                                                                                                                <  2:1001925/150‑1 (MQ=255)
                       atgaCCATCCTCTTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCg                                                                                                                                >  1:749160/1‑150 (MQ=255)
                       atgaCCATCCTCTTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGACATCAATGATGACGATATcaggc                                                                                                                                <  2:1516970/150‑5 (MQ=255)
                        tgaCCATCCTCTTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCgg                                                                                                                               <  2:1204849/150‑1 (MQ=255)
                          aCCATCCTCTTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAg                                                                                                                             >  2:1578629/1‑150 (MQ=255)
                           ccATCCTCTTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGt                                                                                                                            <  1:67803/150‑1 (MQ=255)
                            cATCCTCTTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGtt                                                                                                                           >  1:19998/1‑150 (MQ=255)
                               ccTCTTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAAc                                                                                                                        <  1:774307/150‑1 (MQ=255)
                               ccTCTTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAAc                                                                                                                        <  1:922626/150‑1 (MQ=255)
                                ctctTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACg                                                                                                                       <  1:1254028/150‑1 (MQ=255)
                                 tctTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACgg                                                                                                                      <  2:356038/150‑1 (MQ=255)
                                 tctTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACgg                                                                                                                      <  1:517208/150‑1 (MQ=255)
                                   ttGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTa                                                                                                                    >  2:1356853/1‑150 (MQ=255)
                                    tGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACgg                                                                                                                      >  2:1579645/1‑147 (MQ=255)
                                    tGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACgg                                                                                                                      <  1:1579645/147‑1 (MQ=255)
                                      cTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATcc                                                                                                                 >  1:904685/1‑150 (MQ=255)
                                       tATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCa                                                                                                                <  1:643269/150‑1 (MQ=255)
                                          cGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCTGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCATGATATCGTCATCATTGATGTCGATATCACCGGAGTTAACGGTATCAAGgt                                                                                                             >  1:1300367/1‑150 (MQ=255)
                                           gcagcaATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGgtg                                                                                                            >  2:495380/1‑150 (MQ=255)
                                           gcagcaATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGgtg                                                                                                            <  2:212049/150‑1 (MQ=255)
                                             agcaATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGtt                                                                                                          >  1:402049/1‑150 (MQ=255)
                                                aaTTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGa                                                                                                       <  1:1653822/150‑1 (MQ=255)
                                                     gTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGc                                                                                                  >  1:187543/1‑150 (MQ=255)
                                                      tAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCt                                                                                                 >  2:180285/1‑150 (MQ=255)
                                                       aaTTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTg                                                                                                <  2:183214/150‑1 (MQ=255)
                                                        aTTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGa                                                                                               <  1:1626355/150‑1 (MQ=255)
                                                         tttATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAg                                                                                              <  1:834384/150‑1 (MQ=255)
                                                          ttATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAgg                                                                                             >  1:981292/1‑150 (MQ=255)
                                                             ttGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAg                                                                                          >  1:458800/1‑150 (MQ=255)
                                                              tGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAgc                                                                                         >  2:480537/1‑150 (MQ=255)
                                                               gATCAAAAACGATATTGAAATCATAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAgcg                                                                                        >  1:1139350/1‑150 (MQ=255)
                                                                aTCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAgcgc                                                                                       >  2:1384371/1‑150 (MQ=255)
                                                                  cAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCa                                                                                     <  1:1468638/150‑1 (MQ=255)
                                                                   aaaaaCGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCaa                                                                                    >  2:370433/1‑150 (MQ=255)
                                                                   aaaaaCGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCaa                                                                                    <  1:495380/150‑1 (MQ=255)
                                                                    agaaCGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAat                                                                                   <  2:440208/148‑1 (MQ=255)
                                                                         gATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGc                                                                              <  2:746122/150‑1 (MQ=255)
                                                                               gAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAatt                                                                        <  2:1194049/150‑1 (MQ=255)
                                                                               gAAATCTTAGCAGAGTTGACTGAAGGCGGAAGGGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAatt                                                                        <  2:1281101/150‑1 (MQ=255)
                                                                                   tcatAGCAGAGTTGTCTGAATGCGGAAGTGCCGTGCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCTATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAattatta                                                                    <  1:330280/147‑1 (MQ=255)
                                                                                   tCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAattatta                                                                    <  1:187680/150‑1 (MQ=255)
                                                                                   tCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAattatta                                                                    <  1:526613/150‑1 (MQ=255)
                                                                                       aGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATtattattat                                                                >  1:915499/1‑150 (MQ=255)
                                                                                        gCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGGTATCATTGATGTCGATATCCCCGGAGTTTACGGTATCCAGGTGTTTGAAACGCTGAGGAAGCGCCTATATAGCGGAATTATGATTATc                                                               >  1:265602/1‑150 (MQ=255)
                                                                                        ccAGAGTTGACTGAAGGCGGAAGTGCCGATCCGCGGGTGGAAACACTTAAGCCTGATATCGATATCATTGATCTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAATCGCCAATAAAGCTCAATTATTATTATc                                                               >  1:265603/2‑150 (MQ=255)
                                                                                         cAGAGTTGACTGAAGGCTGTAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGACATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCg                                                              <  1:134685/150‑1 (MQ=255)
                                                                                           gagTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGtc                                                            <  2:599899/150‑1 (MQ=255)
                                                                                             gTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGtctc                                                          >  1:1505039/1‑150 (MQ=255)
                                                                                             gTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGtctc                                                          >  1:953703/1‑150 (MQ=255)
                                                                                               tGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCg                                                        >  2:1209623/1‑150 (MQ=255)
                                                                                                gACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGc                                                       <  2:402049/150‑1 (MQ=255)
                                                                                                    gAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTaaa                                                   >  2:1645375/1‑150 (MQ=255)
                                                                                                          tgAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGAc                                             <  1:1037553/149‑1 (MQ=255)
                                                                                                          ggAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGAc                                             <  2:978814/150‑1 (MQ=255)
                                                                                                               tGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCAttt                                        >  2:323829/1‑150 (MQ=255)
                                                                                                                    ttCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACg                                   <  1:1578629/150‑1 (MQ=255)
                                                                                                                     tCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACgg                                  <  1:389040/150‑1 (MQ=255)
                                                                                                                      tatagtgTGGAAACACTTACGCCTGATATCGTTATCATTGATGTCGATATCGCCGGAGTTAAAGGTATCCAGGTGTTAGAACCGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACggg                                 <  2:265602/144‑1 (MQ=255)
                                                                                                                      taGCGGGTGGAAAGACTTAAGCCTGAGATGGTTATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTTGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACggg                                 <  2:265603/149‑1 (MQ=255)
                                                                                                                       aGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACGGGa                                <  2:19998/150‑1 (MQ=255)
                                                                                                                         cGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACGGGaaa                              >  1:1496712/1‑150 (MQ=255)
                                                                                                                          gggTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACGGGAAAc                             >  2:393256/1‑150 (MQ=255)
                                                                                                                               gAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTCTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACGGGAAACATtgt                        >  2:757992/1‑150 (MQ=255)
                                                                                                                                  acacTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGGTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACGGGAAACATTGTGCt                     <  1:1282678/150‑1 (MQ=255)
                                                                                                                                       tAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGTCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACGGGAAACATTGTGCTGATGc                <  2:1285729/150‑1 (MQ=255)
                                                                                                                                       tAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACGGGAAACATTGTGCTGATGc                >  2:370925/1‑150 (MQ=255)
                                                                                                                                        aaGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACGGGAAACATTGTGCTGATGCt               <  1:1209623/150‑1 (MQ=255)
                                                                                                                                          gCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACGGGAAACATTGTGCTGATGCTgg             >  1:1710917/1‑150 (MQ=255)
                                                                                                                                          gCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACGGGAAACATTGTGCTGATGCTgg             <  2:1534726/150‑1 (MQ=255)
                                                                                                                                          gCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACGGGAAACATTGTGCTGATGCTgg             <  2:1433357/150‑1 (MQ=255)
                                                                                                                                            gtGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACGGGAAACATTGTGCTGATGCTGgcg           <  2:499020/149‑1 (MQ=255)
                                                                                                                                            cTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGACGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACGGGAAACATTGTGCTGATGCTGgcg           >  2:1525358/1‑150 (MQ=255)
                                                                                                                                              gATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACGGGAAACATTGTGCTGATGCTGGCGCt         <  2:1386516/150‑1 (MQ=255)
                                                                                                                                              gATATAGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACGGGAAACATTGTGCTGATGCTGGCGCt         <  1:1384371/150‑1 (MQ=255)
                                                                                                                                               atatCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACGGGAAACATTGTGCTGATGCTGGCGCTa        <  2:953703/150‑1 (MQ=255)
                                                                                                                                               atatCGACATCAATGATGACGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACGGGAAACATTGTGCTGATGCTGGCGCTa        >  1:1281633/1‑150 (MQ=255)
                                                                                                                                                   cGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACGGGAAACATTGTGCTGATGCTGGCGCTAATgg    >  2:643563/1‑150 (MQ=255)
                                                                                                                                                   cGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACGGGAAACATTGTGCTGATGCTGGCGCTAATgg    >  1:1133919/1‑150 (MQ=255)
                                                                                                                                                     tcatcaTTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACGGGAAACATTGTGCTGATGCTGGCGCTAATGGtt  >  1:691125/1‑150 (MQ=255)
                                                                                                                                                     tcatcaTTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACGGGAAACATTGTGCTGATGCTGGCGCTAATGGtt  >  1:1103123/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATCTATGAACGCAATAATTATTGATGACCATCCTCTTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACGGGAAACATTGTGCTGATGCTGGCGCTAATGGTT  >  NC_000913/2483751‑2484049

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