Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,627,717 G→T 10.4% S67S (TCG→TCT ydfG → 3‑hydroxy acid dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,627,7170GT10.4% 318.1 / 15.8 99S67S (TCG→TCTydfG3‑hydroxy acid dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (47/41);  new base T (3/8);  total (50/49)
Fisher's exact test for biased strand distribution p-value = 1.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGAAACGATGATTGATACCAACAACAAAGGCCTGGTATATATGACGCGCGC  >  NC_000913/1627570‑1627863
                                                                                                                                                   |                                                                                                                                                  
gcgtTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCt                                                                                                                                                  <  2:1075239/149‑1 (MQ=255)
  gtTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGCGATGCCAGCATCTCTTc                                                                                                                                                <  1:1320745/150‑1 (MQ=255)
  gtTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTc                                                                                                                                                >  1:946515/1‑150 (MQ=255)
   ttttATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCATCTGGACGTTCGCATCCGCGCCGCTATTGAAGCGATGCCAGCATCTCTTca                                                                                                                                               <  2:386209/150‑2 (MQ=255)
   ttttATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTcc                                                                                                                                               <  2:1130748/150‑1 (MQ=255)
   ttttATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTcc                                                                                                                                               >  2:256908/1‑150 (MQ=255)
   ttttATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGGCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGCGATGCCAGCATCTCTTca                                                                                                                                               <  2:801602/150‑2 (MQ=255)
    tttATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTACAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCt                                                                                                                                              >  1:412184/1‑150 (MQ=255)
     ttATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGCGATGCCAGCATCTCTTcaat                                                                                                                                             <  1:290803/150‑4 (MQ=255)
     ttATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGCGATGCCAGCATCTCTTcaat                                                                                                                                             <  1:836326/150‑4 (MQ=255)
     ttATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTg                                                                                                                                             <  2:1563878/150‑1 (MQ=255)
      tATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGc                                                                                                                                            <  2:1173867/150‑1 (MQ=255)
       aTTCAACAAGGGCATAAAGTTATCGCCACTTTCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTTGCATCGCTTCCTGcc                                                                                                                                           <  2:1496149/150‑1 (MQ=255)
        ttCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGCGATGCCAGCATCTCTTcaatagc                                                                                                                                          <  2:1358631/150‑7 (MQ=255)
        ttCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCg                                                                                                                                          <  2:920162/150‑1 (MQ=255)
        ttCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCg                                                                                                                                          >  2:1628222/1‑150 (MQ=255)
         tCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGCGATGCCAGCATCTCTTcaatagcg                                                                                                                                         <  2:1593830/150‑8 (MQ=255)
          caacaaGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAg                                                                                                                                        <  1:313855/150‑1 (MQ=255)
          caacaaGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGGTAAAAGACGACCTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAg                                                                                                                                        >  2:517230/1‑150 (MQ=255)
          caacaaGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTCGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAg                                                                                                                                        >  2:695623/1‑150 (MQ=255)
           aacaaGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAgt                                                                                                                                       >  1:829106/1‑150 (MQ=255)
            acaaGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCATGCCGAgtg                                                                                                                                      <  2:55480/150‑1 (MQ=255)
             caaGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCg                                                                                                                                          <  2:272807/145‑1 (MQ=255)
             caaGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCg                                                                                                                                          >  1:272807/1‑145 (MQ=255)
              aatgccctAAAGTTATCGCCACTGGCCGTCGCCAGGAAGGGTTGCAGGCGTAAACAGACGAACTGGGTGATAATCTTTATATCGCCGAACTGGAGGTTCGCACCCGCGGCGCTATTGAATAGATGCTGTCATCTCTTCCTGGCGAgttgt                                                                                                                                    <  2:820446/143‑1 (MQ=255)
              aaGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAgtggt                                                                                                                                    >  2:511049/1‑150 (MQ=255)
               aGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAgtggtg                                                                                                                                   >  1:224882/1‑150 (MQ=255)
                  gCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAg                                                                                                                                        <  2:1437514/142‑1 (MQ=255)
                  gCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAg                                                                                                                                        >  1:1437514/1‑142 (MQ=255)
                   cATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAat                                                                                                                               >  2:340935/1‑150 (MQ=255)
                    aTAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAata                                                                                                                              <  2:731646/150‑1 (MQ=255)
                       aaGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTg                                                                                                                           <  1:928015/150‑1 (MQ=255)
                        aGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGa                                                                                                                          >  1:1024090/1‑150 (MQ=255)
                             tCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGatatac                                                                                                                     >  2:1530878/1‑148 (MQ=255)
                                ccACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTgg                                                                                                                  <  1:1432053/150‑1 (MQ=255)
                                 cACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGt                                                                                                                 >  1:1060615/1‑150 (MQ=255)
                                  aCTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCTCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTa                                                                                                                <  1:1525679/150‑1 (MQ=255)
                                  aCTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTa                                                                                                                <  2:1346053/150‑1 (MQ=255)
                                  aCTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTa                                                                                                                <  1:1312228/150‑1 (MQ=255)
                                  aCTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTa                                                                                                                <  2:1527457/150‑1 (MQ=255)
                                   cTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTaa                                                                                                               >  2:251389/1‑150 (MQ=255)
                                     ggCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAaat                                                                                                             >  2:954569/1‑150 (MQ=255)
                                       ccGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAaataa                                                                                                           >  1:763190/1‑150 (MQ=255)
                                             ccAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCgg                                                                                                     >  1:1530397/1‑150 (MQ=255)
                                             ccAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCgg                                                                                                     <  1:227331/150‑1 (MQ=255)
                                                ggAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAgtgg                                                                                                                                     <  2:688246/115‑1 (MQ=255)
                                                ggAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAgtgg                                                                                                                                     >  1:688246/1‑115 (MQ=255)
                                                     ggTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGt                                                                                             >  2:150922/1‑150 (MQ=255)
                                                      gTTGCAGGAGTTAAAAGACGAACTGGTAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGtt                                                                                            <  1:1400331/150‑1 (MQ=255)
                                                         gCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTggg                                                                                         <  1:256908/150‑1 (MQ=255)
                                                           aGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCa                                                                                       <  2:486202/150‑1 (MQ=255)
                                                               gTTAAAATAAAAACTGAGAGATAATCAGTATATCGCACAACTGGAAGTTCGCAACATAGCAGATATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGa                                                                                   >  2:568265/1‑150 (MQ=255)
                                                                 tAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGc                                                                                 >  1:1443871/1‑150 (MQ=255)
                                                                          aaCTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATa                                                                        <  1:944074/150‑1 (MQ=255)
                                                                            cTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATaaa                                                                      <  2:144368/150‑1 (MQ=255)
                                                                            cTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATaaa                                                                      <  2:1060615/150‑1 (MQ=255)
                                                                               ggAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGcc                                                                   <  1:1626708/150‑1 (MQ=255)
                                                                                    tAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTa                                                                                                                <  1:1498706/100‑1 (MQ=255)
                                                                                    tAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTa                                                                                                                >  2:1498706/1‑100 (MQ=255)
                                                                                     aaTCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGtt                                                             >  2:461945/1‑150 (MQ=255)
                                                                                     aaTCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGtt                                                             >  1:1135450/1‑150 (MQ=255)
                                                                                       tCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGa                                                           >  2:794853/1‑150 (MQ=255)
                                                                                       tCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGa                                                           >  2:833021/1‑150 (MQ=255)
                                                                                       tCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGa                                                           <  2:36237/150‑1 (MQ=255)
                                                                                         tGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAg                                                         <  2:375690/150‑1 (MQ=255)
                                                                                            atatCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACt                                                      >  1:1595931/1‑150 (MQ=255)
                                                                                              atCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTgg                                                    <  2:267950/150‑1 (MQ=255)
                                                                                              atCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTgg                                                    >  2:311214/1‑150 (MQ=255)
                                                                                                cGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGa                                                  >  2:1414795/1‑150 (MQ=255)
                                                                                                cGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGa                                                  >  2:633554/1‑150 (MQ=255)
                                                                                                 gCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGaa                                                 >  2:66434/1‑150 (MQ=255)
                                                                                                   ccAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGAAAc                                               <  2:829106/150‑1 (MQ=255)
                                                                                                         ggACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGaa                                                          <  1:1160833/133‑1 (MQ=255)
                                                                                                         ggACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGaa                                                          >  2:1160833/1‑133 (MQ=255)
                                                                                                         ggACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGAAACgatgat                                         >  1:733035/1‑150 (MQ=255)
                                                                                                         ggACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGAAACgatgat                                         >  1:358333/1‑150 (MQ=255)
                                                                                                                cGCAACCGCGCCGCTATTGTAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGAAACGATGATTGATAcc                                  <  2:662846/150‑1 (MQ=255)
                                                                                                                cGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGAAACGATGATTGATAcc                                  >  2:461620/1‑150 (MQ=255)
                                                                                                                cGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGAAACGATGATTGATAcc                                  <  1:326946/150‑1 (MQ=255)
                                                                                                                cGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGAAACGATGATTGATAcc                                  <  1:340935/150‑1 (MQ=255)
                                                                                                                     ccGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGTCTGCGCATAAAGCCAGCGTTGAAGACTGGGAAAc                                               <  2:212041/132‑1 (MQ=255)
                                                                                                                     ccGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGTCTGCGCATAAAGCCAGCGTTGAAGACTGGGAAAc                                               >  1:212041/1‑132 (MQ=255)
                                                                                                                          ccgcTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGAAACGATGATTGATACCAACAACAAAg                        >  2:1221362/1‑150 (MQ=255)
                                                                                                                           tcggcaggAAGCGATGCCAGCATCTCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGAAACGATGATTGATACCAACAACAAAgg                       >  2:597090/8‑150 (MQ=255)
                                                                                                                             ggcaggAAGCGATGCCAGCATCTCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGAAACGATGATTGATACCAACAACAAAGGcc                     >  1:1625203/6‑150 (MQ=255)
                                                                                                                              gcaggAAGCGATGCCAGCATCTCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGAAACGATGATTGATACCAACAACAAAGGCCt                    >  2:1407732/5‑150 (MQ=255)
                                                                                                                               aTTGAAGAGATGCTGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGAAACGATGATTGATACCAACAACAAAGGCCTgg                  >  2:1545855/1‑150 (MQ=255)
                                                                                                                                ttGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGAAACGATGATTGATACCAACAACAAAGGCCTgg                  >  1:625961/1‑150 (MQ=255)
                                                                                                                                ttGAAGAGATGCTGGCACCGCTTCCTGCCGAGTGGGGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGAAACGATGATTGATACCAACAACAAAGGCCTgg                  <  2:406222/150‑1 (MQ=255)
                                                                                                                                  gAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATa                                                                        >  2:438151/1‑94 (MQ=255)
                                                                                                                                  gAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATa                                                                        <  1:438151/94‑1 (MQ=255)
                                                                                                                                  gAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGAAACGATGATTGATACCAACAACAAAGGCCTGGta                >  2:908116/1‑150 (MQ=255)
                                                                                                                                  gAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGAAACGATGATTGATACCAACAACAAAGGCCTGGta                >  1:105503/1‑150 (MQ=255)
                                                                                                                                    agagATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGTAACGATGATTGATACCAACAACAAAGGCCTGGtata              >  1:830951/1‑150 (MQ=255)
                                                                                                                                      agaTGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGAAACGATGATTGATACCAACAACAAAGGCCTGGtatata            <  2:412184/150‑1 (MQ=255)
                                                                                                                                        aTGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGAAACGATGATTGATACCAACAACAAAGGCCTGGTATATATg          <  2:224882/150‑1 (MQ=255)
                                                                                                                                            tGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGAAACGATGATTGATACCAACAACAAAGGCCTGGTATATATGAcgc      <  1:908116/150‑1 (MQ=255)
                                                                                                                                              gCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGCCTGGGAAACGATGAGTGCTACCAACAACAACGGCCTGGTATATATGAcgcgc    >  2:902158/1‑150 (MQ=255)
                                                                                                                                                aTCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGAAACGATGATTGATACCAACAACAAAGGCCTGGTATATATGACgcgcgc  >  1:94191/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGAAACGATGATTGATACCAACAACAAAGGCCTGGTATATATGACGCGCGC  >  NC_000913/1627570‑1627863

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