Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,367,700 G→T 6.2% L210M (CTG→ATG)  aspA ← aspartate ammonia‑lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,367,7000GT 6.2% 321.5 / 6.8 96L210M (CTG→ATG) aspAaspartate ammonia‑lyase
Reads supporting (aligned to +/- strand):  ref base G (47/43);  new base T (4/2);  total (51/45)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.70e-01

AGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTAATCGCATCTACCAGCTTAATC  >  NC_000913/4367553‑4367848
                                                                                                                                                   |                                                                                                                                                    
aGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGa                                                                                                                                                    >  1:629812/1‑150 (MQ=255)
aGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGa                                                                                                                                                    >  2:1264706/1‑150 (MQ=255)
   cAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGTTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGc                                                                                                                                                 <  1:266723/150‑1 (MQ=255)
   cAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGc                                                                                                                                                 <  2:378645/150‑1 (MQ=255)
   cAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGc                                                                                                                                                 <  2:1319719/150‑1 (MQ=255)
       aaCTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGaa                                                                                                                                             <  1:154340/150‑1 (MQ=255)
        aCTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGaaa                                                                                                                                            >  1:1053763/1‑150 (MQ=255)
          ttCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCATCCTGCTGAAAga                                                                                                                                          <  2:391480/150‑2 (MQ=255)
          ttCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTCTTTCGGAGTGTTCAGACCAGTACCGTTTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAgc                                                                                                                                          <  1:623359/150‑1 (MQ=255)
           tCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAgcg                                                                                                                                         <  1:1380323/150‑1 (MQ=255)
            cAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCATCCTGCTGAAAgaag                                                                                                                                        <  2:410868/150‑4 (MQ=255)
              gCCAGTTTTTTCACTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCgg                                                                                                                                      <  2:713822/150‑1 (MQ=255)
                 aGTTTTTTCACTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAAt                                                                                                                                   >  1:844866/1‑150 (MQ=255)
                   ttttttCACTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTc                                                                                                                                 <  2:1281054/150‑1 (MQ=255)
                       ttCACTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGa                                                                                                                             >  1:177965/1‑150 (MQ=255)
                              ccAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAggg                                                                                                                      <  1:1413021/150‑1 (MQ=255)
                                aGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTc                                                                                                                    >  1:1592221/1‑150 (MQ=255)
                                   ggAGAGTACTCTTTCGGCGTGTTCTGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCATTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATc                                                                                                                 <  1:461145/150‑1 (MQ=255)
                                        gTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTAc                                                                                                            >  1:77049/1‑150 (MQ=255)
                                          aCTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTg                                                                                                          >  1:1454131/1‑150 (MQ=255)
                                            tctTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCg                                                                                                        <  2:1593380/150‑1 (MQ=255)
                                             ctTTCGGCCTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACCTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGt                                                                                                       >  2:1208660/1‑150 (MQ=255)
                                             ctTTCGGCATGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCACCGGTACCTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGt                                                                                                       >  2:1214332/1‑150 (MQ=255)
                                               ttCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTcc                                                                                                     <  2:511149/150‑1 (MQ=255)
                                                gcGGCGTTGTCAGACCAGTACCGATTGCTGCTGCACCAAGGTTAACCTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCt                                                                                                    <  2:661261/149‑1 (MQ=255)
                                                  ggCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGc                                                                                                  <  2:123087/150‑1 (MQ=255)
                                                     gtgtTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGc                                                                                               <  2:844866/150‑1 (MQ=255)
                                                       gtTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCATCTg                                                                                             >  1:891305/1‑150 (MQ=255)
                                                       gtTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTg                                                                                             >  2:133492/1‑150 (MQ=255)
                                                       gtTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTg                                                                                             >  1:1140240/1‑150 (MQ=255)
                                                         tCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTGAGCAGGATGCTGAAAGCGCGGAATTACTGACCGAGGGTCATCGGTACTTCGTCCTGAAGCTggt                                                                                           <  2:1191695/150‑2 (MQ=255)
                                                         tCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTggg                                                                                           >  2:386997/1‑150 (MQ=255)
                                                            gACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTAc                                                                                        >  1:606171/1‑150 (MQ=255)
                                                              ccAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGa                                                                                      >  2:940772/1‑150 (MQ=255)
                                                               cAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGAc                                                                                     <  1:1555605/150‑1 (MQ=255)
                                                                aGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGAcc                                                                                    >  2:781423/1‑150 (MQ=255)
                                                                   aCCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCAt                                                                                 >  2:372552/1‑150 (MQ=255)
                                                                    ccGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCCGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCAtt                                                                                <  1:707483/150‑1 (MQ=255)
                                                                    ccGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCAATTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCAtt                                                                                >  2:1448637/1‑150 (MQ=255)
                                                                      gATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAACGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCAgttt                                                                              >  1:1499863/1‑150 (MQ=255)
                                                                         tGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAg                                                                           <  1:1615689/150‑1 (MQ=255)
                                                                         tGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAg                                                                           >  1:888055/1‑150 (MQ=255)
                                                                          gCTGTTGCACCAAGGTTAACTTCCAGCAGGAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAgg                                                                          <  1:466750/150‑1 (MQ=255)
                                                                          gCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACCGCGTCCTGCAGCTGGGTACGACCCATTTTCAgg                                                                          >  2:571053/1‑150 (MQ=255)
                                                                             gTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATg                                                                       >  1:302708/1‑150 (MQ=255)
                                                                             gTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATg                                                                       >  1:676082/1‑150 (MQ=255)
                                                                             gTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGAATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATg                                                                       >  2:566138/1‑150 (MQ=255)
                                                                              ttGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGt                                                                      <  1:314095/150‑1 (MQ=255)
                                                                              ttGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGt                                                                      <  1:243244/150‑1 (MQ=255)
                                                                                gCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTcc                                                                    >  1:525406/1‑150 (MQ=255)
                                                                                 cACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTGTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCt                                                                   >  1:783186/1‑150 (MQ=255)
                                                                                   ccAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTgg                                                                 >  1:394919/1‑150 (MQ=255)
                                                                                        gTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTc                                                            <  1:680497/150‑1 (MQ=255)
                                                                                        gTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTc                                                            >  1:984766/1‑150 (MQ=255)
                                                                                         ttAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCg                                                           <  1:1264706/150‑1 (MQ=255)
                                                                                           aaCTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGAc                                                         <  2:629812/150‑1 (MQ=255)
                                                                                             cTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAg                                                       >  2:1127602/1‑150 (MQ=255)
                                                                                             cTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAg                                                       >  2:753256/1‑150 (MQ=255)
                                                                                             cTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTACTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAg                                                       >  1:698221/1‑150 (MQ=255)
                                                                                             cTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAACCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAg                                                       >  1:144038/1‑150 (MQ=255)
                                                                                              ttCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGc                                                      <  2:386999/150‑1 (MQ=255)
                                                                                                   gcagcagTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTAc                                                 <  2:1453658/150‑1 (MQ=255)
                                                                                                    cagcagTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACg                                                >  1:583272/1‑150 (MQ=255)
                                                                                                     agcagTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGt                                               >  2:1001040/1‑150 (MQ=255)
                                                                                                           tCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAg                                         <  1:1188851/150‑1 (MQ=255)
                                                                                                              gcagTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGc                                                      >  2:428561/4‑134 (MQ=255)
                                                                                                              gcagTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGc                                                      <  1:428561/131‑1 (MQ=255)
                                                                                                              gCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCt                                      <  2:212382/150‑1 (MQ=255)
                                                                                                               cGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACCGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCtt                                     <  1:571053/150‑1 (MQ=255)
                                                                                                                 gTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGGCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCa                                   <  1:678735/150‑1 (MQ=255)
                                                                                                                   aCGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACg                                 <  2:606171/150‑1 (MQ=255)
                                                                                                                      ttGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAg                              <  2:984766/150‑1 (MQ=255)
                                                                                                                         gATGTTTTTCACTTCTTCTTTCAGCAGAATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTg                           <  1:566138/150‑1 (MQ=255)
                                                                                                                           tGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGt                         >  2:977377/1‑150 (MQ=255)
                                                                                                                             tttttCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTa                       >  1:40049/1‑150 (MQ=255)
                                                                                                                             tttttCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTa                       >  1:927557/1‑150 (MQ=255)
                                                                                                                                ttCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCCCGCAGTTGGTTAATc                    >  2:512423/1‑150 (MQ=255)
                                                                                                                                 tCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTAATCg                   <  1:71792/150‑1 (MQ=255)
                                                                                                                                   aCTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTAATCGCa                 <  2:1532107/150‑1 (MQ=255)
                                                                                                                                   aCTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTAATCGCa                 >  1:1176590/1‑150 (MQ=255)
                                                                                                                                      cgcgcttTCAGCATCCTGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGCCAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTAATCGCATCt              >  1:719336/5‑150 (MQ=255)
                                                                                                                                       cttcttTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTAATCGCATCtg             >  1:1042281/1‑149 (MQ=255)
                                                                                                                                       cttcttTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTAATCGCATCTa             >  2:551219/1‑150 (MQ=255)
                                                                                                                                        cgcttTCAGCATCCTGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTAATCGCATCTAc            >  1:20290/3‑150 (MQ=255)
                                                                                                                                         tcttTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTAATCGCATCTAcc           <  1:1001040/150‑1 (MQ=255)
                                                                                                                                          cttTCAGCATCCTGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTAATCGCATCTACCa          >  1:1201843/1‑150 (MQ=255)
                                                                                                                                          cttTCAGCATCCTGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTAATCGCATCTACCa          >  1:1365052/1‑150 (MQ=255)
                                                                                                                                          cttTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACg                                 <  1:663759/127‑1 (MQ=255)
                                                                                                                                          cttTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACg                                 >  2:663759/1‑127 (MQ=255)
                                                                                                                                          cttTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTAATCGCATCTACCa          >  2:310694/1‑150 (MQ=255)
                                                                                                                                          cttTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTAATCGCATCTACCa          >  2:491127/1‑150 (MQ=255)
                                                                                                                                           tttCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTAATCGCATCTACCAg         <  1:1448637/150‑1 (MQ=255)
                                                                                                                                            ttCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTAATCGCATCTACCAGc        >  1:1321155/1‑150 (MQ=255)
                                                                                                                                            ttCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTAATCGCATCTACCAGc        <  2:583272/150‑1 (MQ=255)
                                                                                                                                            ttCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTAATCGCATCTACCAGc        <  1:1214332/150‑1 (MQ=255)
                                                                                                                                            ttCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTAATCGCATCTACCAGc        <  1:1208660/150‑1 (MQ=255)
                                                                                                                                              cagcagGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTAATCGCATCTACCAGCtt      >  1:650318/1‑150 (MQ=255)
                                                                                                                                                gcagGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTAATCGCATCTACCAGCTTaa    >  1:226525/1‑150 (MQ=255)
                                                                                                                                                  agGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTAATCGCATCTACCAGCTTAATc  >  2:1078548/1‑150 (MQ=255)
                                                                                                                                                  agGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTAATCGCATCTACCAGCTTAATc  >  2:340105/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTAATCGCATCTACCAGCTTAATC  >  NC_000913/4367553‑4367848

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: ‑