Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,829,556 T→C 5.4% E60G (GAA→GGA)  astD ← succinylglutamic semialdehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,829,5560TC 5.4% 231.5 / 2.2 74E60G (GAA→GGA) astDsuccinylglutamic semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (35/35);  new base C (2/2);  total (37/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCC  >  NC_000913/1829412‑1829701
                                                                                                                                                |                                                                                                                                                 
gATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAg                                                                                                                                                >  4:75985/1‑148 (MQ=255)
    tCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAAt                                                                                                                                           >  3:85989/1‑149 (MQ=255)
     cgcgATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAAt                                                                                                                                           >  1:348551/1‑148 (MQ=255)
        gATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGaaa                                                                                                                                       <  7:157126/149‑1 (MQ=255)
               ttGATCATCGCCGTCACTTCGGTTGACGCTTAACCGCGCGGCTTAACAGTGTCTCTGGCAATAAGCGCGGTTAATTCGTCTTTATTGCTTTCCAGCAGTGCGGCAAAGAGTTCGACAAAGGCATGACGTTCAGCAAATGAAAGAcgctc                                                                                                                                >  6:336041/1‑147 (MQ=255)
                 gATCATCGCCGTCACTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGcc                                                                                                                               >  5:153696/1‑148 (MQ=255)
                  atcatcGCCGTCACTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTCCAGCAAATGAAAGCCGCGcc                                                                                                                               >  3:172797/1‑147 (MQ=255)
                  atcatcGCCGTCACTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTCCAGCAAATGAAAGCCGCGCCCa                                                                                                                             >  3:167390/1‑149 (MQ=255)
                   tcatcGCCGTCACTTCGGTTGCCGCTTCCCAGCGCGGCTTCCCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAg                                                                                                                            >  5:54034/1‑149 (MQ=255)
                      tcGCCGTCACTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAgcgc                                                                                                                         >  1:69075/1‑149 (MQ=255)
                         ccGTCACTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCgg                                                                                                                       >  1:268088/1‑148 (MQ=255)
                            tCACTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAAc                                                                                                                   >  8:330941/1‑149 (MQ=255)
                               cTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCTGGGCGCggctttcat                                                                                                                <  3:75985/149‑8 (MQ=255)
                               cTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAAcgcc                                                                                                                <  7:178798/149‑1 (MQ=255)
                                ttCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAAcgccg                                                                                                               >  4:52149/1‑149 (MQ=255)
                                ttCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAAcgccg                                                                                                               >  7:326092/1‑149 (MQ=255)
                                         cGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTg                                                                                                      >  7:68212/1‑149 (MQ=255)
                                          gCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGc                                                                                                     >  7:113786/1‑149 (MQ=255)
                                              cccAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGccc                                                                                                   <  2:69075/147‑1 (MQ=255)
                                              cccAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGa                                                                                                 >  6:292328/1‑149 (MQ=255)
                                                 aGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACaa                                                                                              <  4:249321/149‑1 (MQ=255)
                                                  gcgcGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAAcgccgc                                                                                                              <  8:260830/132‑1 (MQ=255)
                                                  gcgcGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAAcgccgc                                                                                                              >  7:260830/1‑132 (MQ=255)
                                                   cgcgGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGc                                                                                            <  6:153696/149‑1 (MQ=255)
                                                   cgcgGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGc                                                                                            <  5:292328/149‑1 (MQ=255)
                                                            cccGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGAc                                                                                   <  6:265654/149‑1 (MQ=255)
                                                             ccGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGAcc                                                                                  >  7:215482/1‑149 (MQ=255)
                                                               gTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGAc                                                                                   <  7:37043/146‑1 (MQ=255)
                                                               gTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGAc                                                                                   >  8:37043/1‑146 (MQ=255)
                                                               gTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCAGACAAGCCTGCTCGACCt                                                                                 >  7:260408/1‑148 (MQ=255)
                                                                        gCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGcggcatc                                                                       <  8:80628/149‑1 (MQ=255)
                                                                          aataatCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGcggcatcgg                                                                     >  6:301873/1‑149 (MQ=255)
                                                                           ataatCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACg                                                                                                           >  7:183756/1‑110 (MQ=255)
                                                                           ataatCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACg                                                                                                           <  8:183756/110‑1 (MQ=255)
                                                                                 gctgTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCAtt                                                              <  3:145011/146‑1 (MQ=255)
                                                                                 gcgGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCAtt                                                              <  3:57501/149‑1 (MQ=255)
                                                                                 gcgGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTCCTCGACCTGAGCGGCATCGGCATCAtt                                                              >  4:136363/1‑149 (MQ=255)
                                                                                 gcgGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCCCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGcatcat                                                               <  3:108418/148‑1 (MQ=255)
                                                                                  cgGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCAttg                                                             >  3:2422/1‑149 (MQ=255)
                                                                                     ttAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGtt                                                                                                                                                   >  7:126534/1‑60 (MQ=255)
                                                                                     ttAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGtt                                                                                                                                                   <  8:126534/60‑1 (MQ=255)
                                                                                        aTTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCAttgccttg                                                        <  7:273571/148‑1 (MQ=255)
                                                                                        aTTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACACCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCAttgccttg                                                        <  2:137905/148‑1 (MQ=255)
                                                                                         ttCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCAtt                                                              >  8:17364/1‑141 (MQ=255)
                                                                                         ttCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCAtt                                                              <  7:17364/141‑1 (MQ=255)
                                                                                             gCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGCCAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATaa                                                  >  2:160846/1‑149 (MQ=255)
                                                                                              cTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAac                                                 <  2:348551/149‑1 (MQ=255)
                                                                                                ttATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAacac                                               >  3:43860/1‑149 (MQ=255)
                                                                                                  aTTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTCCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACAcc                                              <  4:167390/148‑1 (MQ=255)
                                                                                                  aTTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTCCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACAcc                                              <  4:172797/148‑1 (MQ=255)
                                                                                                   ttGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTc                                            <  3:52149/149‑1 (MQ=255)
                                                                                                         tCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGa                                      >  6:189443/1‑149 (MQ=255)
                                                                                                           cagcagTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCAttgcctt                                                         >  4:143054/1‑128 (MQ=255)
                                                                                                           cagcagTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCAttgcctt                                                         <  3:143054/128‑1 (MQ=255)
                                                                                                           cagcagTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCAttg                                                             <  3:282677/124‑1 (MQ=255)
                                                                                                           cagcagTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCAttg                                                             >  4:282677/1‑124 (MQ=255)
                                                                                                             gcagTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATAcc                                  >  8:276601/1‑149 (MQ=255)
                                                                                                                 tGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGAt                              <  8:68212/149‑1 (MQ=255)
                                                                                                                        aaGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCtt                       <  7:330941/149‑1 (MQ=255)
                                                                                                                         aGCGTTCGACACCGGGATGACGTTCAGCAAATGAAAGCCGCGCCCCGCGCGGAAACGCCGCGCGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCCTCATTGCCTTGCCATAACACCTCGCCCGAGACCCGATTCCGCTTc                      >  2:128547/1‑149 (MQ=255)
                                                                                                                            gTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGcc                                         <  1:322487/127‑1 (MQ=255)
                                                                                                                            gTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGcc                                         >  2:322487/1‑127 (MQ=255)
                                                                                                                            gTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTcaca                   <  8:215482/149‑1 (MQ=255)
                                                                                                                                  cAACGGCATGACGTTCAGCAAATGAACGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGAt             <  8:113786/149‑1 (MQ=255)
                                                                                                                                   aaCGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACAcgc                >  7:30470/1‑145 (MQ=255)
                                                                                                                                   aaCGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACAcgc                <  8:30470/145‑1 (MQ=255)
                                                                                                                                    aCGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATgc           <  5:301873/149‑1 (MQ=255)
                                                                                                                                        cATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGccc       >  3:185825/1‑149 (MQ=255)
                                                                                                                                         aTGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCt      <  4:43860/149‑1 (MQ=255)
                                                                                                                                           gACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTg     <  6:289288/148‑1 (MQ=255)
                                                                                                                                           gACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTg     >  5:289288/1‑148 (MQ=255)
                                                                                                                                            aCGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAATCGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCAGCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACCCGCGATGCGCCCTgg    <  6:73214/148‑1 (MQ=255)
                                                                                                                                            aCGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGc   <  8:326092/149‑1 (MQ=255)
                                                                                                                                             cGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGcc  <  6:148278/149‑1 (MQ=255)
                                                                                                                                              gTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATaa                                                  <  6:356351/100‑1 (MQ=255)
                                                                                                                                              gTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATaa                                                  >  5:356351/1‑100 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
GATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCTTCCCAGCGCGGCTTACCCGTTTCTCTGGCAATAATCGCGGTTAATTCGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTCACACGCGATGCGCCCTGGCC  >  NC_000913/1829412‑1829701

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