Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 699,134 A→G 5.6% V15A (GTT→GCT)  asnB ← asparagine synthetase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913699,1340AG 5.6% 216.1 / 2.9 71V15A (GTT→GCT) asnBasparagine synthetase B
Reads supporting (aligned to +/- strand):  ref base A (31/36);  new base G (2/2);  total (33/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

TAGAGAGGTTGCGCCCCCGCGTTAACGTCAACAATTGACAGACGTTCGTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCTTGTTGTGTTTGCTTGTTTAAAAAAATGCCGCAAAGCAGCACTGTGCGCAAGCGATTTGGCGGTGAAAAAATAAAA  >  NC_000913/698987‑699283
                                                                                                                                                   |                                                                                                                                                     
tAGAGAGGTTGCGCCCCCGCGTTAACGTCAACAATTGACAGACGTTCGTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACt                                                                                                                                                     <  1:67859/150‑1 (MQ=255)
   agagGTTGCGCCCCCGCGTTAACGTCAACAATTGACAGACGTTCGTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGTTGAGCTGCg                                                                                                                                                  <  1:1568445/150‑1 (MQ=255)
   agagGTTGCGCCCCCGCGTTAACGTCAACAATTGACAGACGTTCGTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCg                                                                                                                                                  <  2:67609/150‑1 (MQ=255)
   agagGTTGCGCCCCCGCGTTAACGTCAACAATTGACAGACGTTCGTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCg                                                                                                                                                  <  1:1926189/150‑1 (MQ=255)
     agGTTGCGCCCCCGCGTTAACGTCAACAATTGACAGACGTTCGTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTc                                                                                                                                                >  1:265467/1‑150 (MQ=255)
      ggTTGCGCCCCCGCGTTAACGTCAACAATTGACAGACGTTCGTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCt                                                                                                                                               >  1:1519517/1‑150 (MQ=255)
      ggTTGCGCCCCCGCGTTAACGTCAACAATTGACAGACGTTCGTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCt                                                                                                                                               <  1:235719/150‑1 (MQ=255)
          gcgcCCCCGCGTTAACGTCAACAATTGACAGACGTTCGTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGttt                                                                                                                                           <  2:265467/150‑1 (MQ=255)
           cgcCCCCGCGTTAACGTCAACAATTGACAGACGTTCGTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGTTGAGCTGCGtaagaaag                                                                                                                                          <  1:79180/150‑8 (MQ=255)
           cgcCCCCGCGTTAACGTCAACAATTGACAGACGTTCGTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGtttt                                                                                                                                          >  2:1109161/1‑150 (MQ=255)
            gcCCCCGCGTTAACGTCAACAATTGACAGACGTTCGTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTg                                                                                                                                         <  1:31485/150‑1 (MQ=255)
              ccccGCGTTAACGTCAACAATTGACAGACGTTCGTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGat                                                                                                                                       >  2:1116477/1‑150 (MQ=255)
               cccGCGTTAACGTCAACAATTGACAGACGTTCGTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGata                                                                                                                                      >  2:1749182/1‑150 (MQ=255)
                    gTTAACGTCAACAATTGACAGACGTTCGTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCg                                                                                                                                   >  2:686035/1‑148 (MQ=255)
                    gTTAACGTCAACAATTGACAGACGTTCGTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCg                                                                                                                                   <  1:686035/148‑1 (MQ=255)
                     ttAACGTCAACAATTGACAGACGTTCGTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAAt                                                                                                                                <  2:1055463/150‑1 (MQ=255)
                          gTCAACAATTGACAGACGTTCGTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGcc                                                                                                                           >  1:355215/1‑150 (MQ=255)
                                    gACAGACGTTCGTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGaa                                                                                                                 <  1:642927/150‑1 (MQ=255)
                                     aCAGACGTTCGTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAac                                                                                                                <  2:1082872/150‑1 (MQ=255)
                                        gACGTTCGTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATAAGGCGTGACAGCTCGAGGGCTTTCTTACGAAGCTCAACTGCGTCTGTATTGATATCGATTACGGCAAAAATTGagcgca                                                                                                             >  2:1583898/1‑145 (MQ=255)
                                         aCGTTCGTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACAt                                                                                                            <  2:1466580/150‑1 (MQ=255)
                                          cGTTCGTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATa                                                                                                           >  2:1033043/1‑150 (MQ=255)
                                           gTTCGTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATaa                                                                                                          >  2:267774/1‑150 (MQ=255)
                                           gTTCGTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATaa                                                                                                          >  2:242384/1‑150 (MQ=255)
                                            ttCGTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAAc                                                                                                         >  2:1655800/1‑150 (MQ=255)
                                             tCGTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAAcc                                                                                                        >  2:286235/1‑150 (MQ=255)
                                              cGTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCt                                                                                                       <  2:1785616/150‑1 (MQ=255)
                                               gTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCtt                                                                                                      <  1:1280874/150‑1 (MQ=255)
                                                tGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTtc                                                                                                     >  1:980543/1‑150 (MQ=255)
                                                  ggCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTtctc                                                                                                   <  1:1706548/150‑1 (MQ=255)
                                                      agaATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGtt                                                                                               >  1:1425740/1‑150 (MQ=255)
                                                       taaTGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTa                                                                                              <  2:1715946/149‑1 (MQ=255)
                                                       gaATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTa                                                                                              <  1:678760/150‑1 (MQ=255)
                                                              gTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGTCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGt                                                                                       <  1:1304757/150‑1 (MQ=255)
                                                              gTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGt                                                                                       <  1:341887/150‑1 (MQ=255)
                                                                   cGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATttg                                                                                  <  2:679284/150‑1 (MQ=255)
                                                                   cGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATttg                                                                                  <  2:1581416/150‑1 (MQ=255)
                                                                     cTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATttgtt                                                                                <  2:1519517/150‑1 (MQ=255)
                                                                      tGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATttgttg                                                                               <  2:1656475/150‑1 (MQ=255)
                                                                          aTAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCttg                                                                           <  2:1425740/150‑1 (MQ=255)
                                                                           tAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCttgt                                                                          <  1:1116477/150‑1 (MQ=255)
                                                                               tACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCTTGTtgtg                                                                      <  2:1546859/150‑1 (MQ=255)
                                                                                         gTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCTTGTTGTGTTTGCTTGtt                                                            >  1:1891738/1‑150 (MQ=255)
                                                                                           ccGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCTTGTTGTGTTTGCTTGTTTa                                                          <  1:242384/150‑1 (MQ=255)
                                                                                                 cACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCTTGTTGTGTTTGCTTGTTTaaaaaaa                                                    <  2:1387679/150‑1 (MQ=255)
                                                                                                  aCGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCTTGTTGTGTTTGCTTGTTTAAAAAAAt                                                   >  2:387144/1‑150 (MQ=255)
                                                                                                     aTGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCTTGTTGTGTTTGCTTGTTTAAAAAAATGcc                                                >  1:629921/1‑150 (MQ=255)
                                                                                                        aCGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCTTGTTGTGTTTGCTTGTTTAAAAAAATGCCGCa                                             <  2:629921/150‑1 (MQ=255)
                                                                                                         cGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCTAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCTTGTTGTGTTTGCTTGTTTAAAAAAATGCCGCaa                                            >  1:1968006/1‑150 (MQ=255)
                                                                                                           cATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCTTGTTGTGTTTGCTTGTTTAAAAAAATGCCGCAAag                                          >  1:326971/1‑150 (MQ=255)
                                                                                                           cATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTCCTCCGTTAACCTGGTATTTGTTGCTTGTTGTGTTTGCTTGTTTAAAAAAATGCCGCAAag                                          >  2:333968/1‑150 (MQ=255)
                                                                                                            aTCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCTTGTTGTGTTTGCTTGTTTAAAAAAATGCCGCAAagc                                         <  1:1109161/150‑1 (MQ=255)
                                                                                                                      aCAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCTTGTTGTGTTTGCTTGTTTAAAAAAATGCCGCAAAGCAGCACTGTgc                               >  1:1265419/1‑150 (MQ=255)
                                                                                                                       cAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCTTGTTGTGTTTGCTTGTTTAAAAAAATGCCGCAAAGCAGCACTGTgcg                              >  2:384790/1‑150 (MQ=255)
                                                                                                                       cAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCTTGTTGTGTTTGCTTGTTTAAAAAAATGCCGCAAAGCAGCACTGTgcg                              <  2:980543/150‑1 (MQ=255)
                                                                                                                          cTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCTTGTTGTGTTTGCTTGTTTAAAAAAATGCCGCAAAGCAGCACTGTGCGCaa                           >  1:1078612/1‑150 (MQ=255)
                                                                                                                          cTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCTTGTTGTGTTTGCTTGTTTAAAAAAATGCCGCAAAGCAGCACTGTGCGCaa                           >  2:399199/1‑150 (MQ=255)
                                                                                                                          cTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCTTGTTGTGTTTGCTTGTTTAAAAAAATGCCGCAAAGCAGCACTGTGCGCaa                           >  2:44884/1‑150 (MQ=255)
                                                                                                                               gggCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCTTGTTGTGTTTGCTTGTTTAAAAAAATGCCGCAAAGCAGCACTGTGCGCAAGCGAt                      >  1:1099902/1‑150 (MQ=255)
                                                                                                                                  cTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCTTGTTGTGTTTGCTTGTTTAAAAAAATGCCGCAAAGCAGCACTGTGCGCAAGCGATTTg                   >  1:1438027/1‑150 (MQ=255)
                                                                                                                                   tttCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCTTGTTGTGTTTGCTTGTTTAAAAAAATGCCGCAAAGCAGCACTGTGCGCAAGCGATTTgg                  <  1:387144/150‑1 (MQ=255)
                                                                                                                                    ttCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCTTGTTGTGTTTGCTTGTTTAAAAAAATGCCGCAAAGCAGCACTGTGCGCAAGCGATTTGGc                 <  2:1120886/150‑1 (MQ=255)
                                                                                                                                    aacagaCGCAGTTGAGCTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCTTGTTGTGTTTGCTTGTTTAAAAAAATGCCGCAAAGCAGCACTGTGCGCAAGCGATTTGGc                 >  2:1041135/6‑150 (MQ=255)
                                                                                                                                     tCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCTTGTTGTGTTTGCTTGTTTAAAAAAATGCCGCAAAGCAGCACTGTGCGCAAGCGATTTGGCg                >  2:132932/1‑150 (MQ=255)
                                                                                                                                      cagaCGCAGTTGAGCTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCTTGTTGTGTTTGCTTGTTTAAAAAAATGCCGCAAAGCAGCACTGTGCGCAAGCGATTTGGCgg               >  2:1323641/4‑150 (MQ=255)
                                                                                                                                        tACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCTAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCTTGTTGTGTTTGCTTGTTTAAAAAAATGCCGCAAAGCAGCACTGTGCGCAAGCGATTTGGCGGTg             <  2:1968006/150‑1 (MQ=255)
                                                                                                                                          cGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCTTGTTGTGTTTGCTTGTTTAAAAAAATGCCGCAAAGCAGCACTGTGCGCAAGCGATTTGGCGGTGaa           <  2:1891738/150‑1 (MQ=255)
                                                                                                                                           gCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCTTGTTGTGTTTGCTTGTTTAAAAAAATGCCGCAAAGCAGCACTGTGCGCAAGCGATTTGGCGGTGaaa          >  2:1201140/1‑150 (MQ=255)
                                                                                                                                                tCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCTTGTTGTGTTTGCTTGTTTAAAAAAATGCCGCAAAGCAGCACTGTGCGCAAGCGATTTGGCGGTGAAAAAATa     <  1:1655800/150‑1 (MQ=255)
                                                                                                                                                 cAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCTTGTTGTGTTTGCTTGTTTAAAAAAATGCCGCAAAGCAGCACTGTGCGCAAGCGATTTGGCGGTGAAAAAATaa    <  1:1415900/150‑1 (MQ=255)
                                                                                                                                                 cAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCTTGTTGTGTTTGCTTGTTTAAAAAAATGCCGCAAAGCAGCACTGTGCGCAAGCGATTTGGCGGTGAAAAAATaa    >  1:305585/1‑150 (MQ=255)
                                                                                                                                                   aCTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCTTGTTGTGTTTGCTTGTTTAAAAAAATGCCGCAAAGCAGCACTGTGCGCAAGCGATTTGGCGGTGAAAAAATaaaa  >  2:571659/1‑150 (MQ=255)
                                                                                                                                                   aCTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCTTGTTGTGTTTGCTTGTTTAAAAAAATGCCGCAAAGCAGCACTGTGCGCAAGCGATTTGGCGGTGAAAAAATaaaa  <  1:1715398/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TAGAGAGGTTGCGCCCCCGCGTTAACGTCAACAATTGACAGACGTTCGTGGGCGAGAATGGCGTTATCGCTGGCATAAATACCGGACCAGTCCGGGCCACGATGACGCATCAGGCGTGACAGCTCGAGGGCTTTCTTACGCAGCTCAACTGCGTCTGTTTTGATATCGAATACGCCAAAAATTGAACACATAACCTTCTCCGTTAACCTGGTATTTGTTGCTTGTTGTGTTTGCTTGTTTAAAAAAATGCCGCAAAGCAGCACTGTGCGCAAGCGATTTGGCGGTGAAAAAATAAAA  >  NC_000913/698987‑699283

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