Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,058,222 G→A 5.4% T63T (ACC→ACT)  serA ← D‑3‑phosphoglycerate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,058,2220GA 5.4% 240.5 / 2.2 76T63T (ACC→ACT) serAD‑3‑phosphoglycerate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (29/42);  new base A (2/2);  total (31/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.83e-01

ACAGAGCGCGTATTTGAGAACGGTGCGTTAAATACCGGGATCCCGCGCTTTGCCGCCGCATCCAGATCAACCTGGTTTGTTCCGATACAGAAACAGCCAATAGCGACCAGTTTTTCTGCGGCGTTGATCACGTCTTCAGTCAGATGGGTACGGGATCGCAGGCCGATGAAGTGGGCATCGCGGATGGATTCTTTTAATTGTTCATCATCCAGCGCGCCTTTGTGAAATTCGATGTTGGTGTAACCAGCTGCACGAAGGCTTTCCAGCGCCTTTTGGTGCACGCCTTCTACCAG  >  NC_000913/3058076‑3058368
                                                                                                                                                  |                                                                                                                                                  
aCAGAGCGCGTATTTGAGAACGGTGCGTTAAATACCGGGATCCCGCGCTTTGCCGCCGCATCCAGATCAACCTGGTTTGTTCCGATACAGAAACAGCCAATAGCGACCAGTTTTTCTGCGGCGTTGATCACGTCTTCAGTCAGATGGGt                                                                                                                                                  <  1:203244/149‑1 (MQ=255)
 cAGAGCGCGTATTTGAGAACGGTGCGTTAAATACCGGGATCCCGCGCTTTGCCGCCGCATCCAGATCAACCTGGTTTGTTCCGATACAGAAACAGCCAATAGCGACCAGTTTTTCTGCGGCGTTGATCACGTCTTCAGTCAGATGGGTa                                                                                                                                                 <  6:208651/149‑1 (MQ=255)
   gagCGCTTATTTGAGAACGGTGCGTTAAATACCGGGATCCCGCGCTTTGCCGCCGCATCCAGATCAACGTGGTTTGTTCCGATACAGATACAGCCAATAGCGACCTGTTTTGCTGCGGCGTTGATCACTTCTTAAGTCAGATGGGTTCg                                                                                                                                               <  6:200095/149‑1 (MQ=255)
      cgcgTATTTGAGAACGGTGCGTTAAATACCGGGATCCCGCGCTTTGCCGCCGCATCCAGATCAACCTGGTTTGTTCCGATACAGAAACAGCCAATAGCGACCAGTTTTTCTGCGGCGTTGATCACGTCTTCAGTCAGATGGGTACGGGa                                                                                                                                            <  1:287511/149‑1 (MQ=255)
        cgTATTTGAGAACGGTGCGTTAAATACCGGGATCCCGCGCTTTGCCGCCGCATCCAGATCAACCTGGTTTGTTCCGATACAGAAACAGCCAATAGCGACCAGTTTTTCTGCGGCGTTGATCACGTCTTCAGTCAGATGGGTACGGGATc                                                                                                                                          >  4:151/1‑149 (MQ=255)
            tttGAGAACGGTGCGTTAAATACCGGGATCCCGCGCTTTGCCGCCGCATCCAGATCAACCTGGTTTGTTCCGATACAGAAACAGCCAATAGCGACCAGTTTTTCTGCGGCGTTGATCACGTCTTCAGTCAGATGGGTACGGGATCGCa                                                                                                                                       <  1:175639/148‑1 (MQ=255)
            tttGAGAACGGTGCGTTAAATACCGGGATCCCGCGCTTTGCCGCCGCATCCAGATCAACCTGGTTTGTTCCGATACAGAAACAGCCAATAGCGACCAGTTTTTCTGCGGCGTTGATCACGTCTTCAGTCAGATGGGTACGGGATCGCa                                                                                                                                       <  1:55667/148‑1 (MQ=255)
            tttGAGAACGGTGCGTTAAATACCGGGATCCCGCGCTTTGCCGCCGCATCCAGATCAACCTGGTTTGTTCCGATACAGAAACAGCCAATAGCGACCAGTTTTTCTGCGGCGTTGATCACGTCTTCAGTCAGATGGGTACGGGATCGCa                                                                                                                                       <  2:374004/148‑1 (MQ=255)
              tGAGAACGGTGCGTTAAATACCGGGATCCCGCGCTTTGCCGCCGCAGCCAGATCAACCCGGTTTGTTCCGATACAGAAACAGCCAATAGCGACCAGTTTTTCTGCGGCGTTGATCACGTCTTCAGTCAGATGGGTACGGGATCGCAgg                                                                                                                                     <  7:11410/148‑1 (MQ=255)
               gagaACGGTGCGTTAAATACCGGGATCCCGCGCTTTGCCGCCGCATCCAGATCAACCTGGTTTGTTCCGATACAGAAACAGCCAATAGCGACCAGTTTTTCTGCGGCGTTGATCACGTCTTCAGTCAGATGGGTACGGGATCGCAGGc                                                                                                                                    <  8:124913/148‑1 (MQ=255)
                     ggTGCGTTAAATACCGGGATCCCGCGCTTTGCCGCCGCATCCAGATCAACCTGGTTTGTTCCGATACAGAAACAGCCAATAGCGACCAGTTTTTCTGCGGCGTTGATCACGTCTTCAGTCAGATGGGTACGGGATCGCAGGCCGATGaa                                                                                                                             <  7:83912/149‑1 (MQ=255)
                         cGTTAAATACCGGGATCCCGCGCTTTGCCGCCGCATCCAGATCAACCTGGTTTGTTCCGATACAGAAACAGCCAATAGCGACCAGTTTTTCTGCGGCGTTGATCACGTCTTCAGTCAGATGGGTACGGGATCGCAGGCCGATGAAGTgg                                                                                                                         >  8:226078/1‑149 (MQ=255)
                         cGTTAAATACCGGGATCCCGCGCTTTGCCGCCGCATCCAGAGCAACCTGGTTTGTTCCGATACAGAAACAGCCAATAGCGACCAGTTTTTCTGCGGCGTTGATCACGTCTTCAGTCAGATGGGTACGGGATCGCAGGCCGATGAAGTgg                                                                                                                         <  2:38135/149‑1 (MQ=255)
                           ttAAATACCGGGATCCCGCGCTTTGCCGCCGCATCCAGATCAACCTGGTTTGTTCCGATACAGAAACAGCCAATAGCGACCAGTTTTTCTGCGGCGTTGATCACGTCTTCAGTCAGATGGGTACGGGATCGCAGGCCGATGAAGTggg                                                                                                                        >  4:286257/1‑148 (MQ=255)
                           ttAAATACCGGGATCCCGCGCTTTGCCGCCGCATCCAGATCAACCTGGTTTGTTCCGATACAGAAACAGCCAATAGCGACCAGTTTTTCTGCGGCGTTGATCACGTCTTCAGTCAGATGGGTACGGGATCGCAGGCCGATGAAGTggg                                                                                                                        >  6:246585/1‑148 (MQ=255)
                           ttAAATACCGGGATCCCGCGCTTTGCCGCCGCATCCAGATCAACCTGGTTTGTTCCGATACAGAAACAGCCAATAGCGACCAGTTTTTCTGCGGCGTTGATCACGTCTTCAGTCAGATGGGTACGGGATCGCAGGCCGATGAAGTggg                                                                                                                        >  8:386677/1‑148 (MQ=255)
                               aTACCGGGATCCCGCGCTTTGCCGCCGCATCCAGATCAACCTGGTTTGTTCCGATACAGAAACAGCCAATAGCGACCAGTTTTTCTGCGGCGTTGATCACGTCTTCAGTCAGATGGGTACGGGATCGCAGGCCGATGAAGTGGGCATcg                                                                                                                   <  5:301132/149‑1 (MQ=255)
                                tACCGGGATCCCGCGCTTTGCCGCCGCATCCAGATCAACCTGGTTTGTTCCGATACAGAAACAGCCAATAGCGACCAGTTTTTCTGCGGCGTTGATCACGTCTTCAGTCAGATGGGTACGGGATCGCAGGCCGATGAAGTGGGCATcgc                                                                                                                  <  7:394985/149‑1 (MQ=255)
                                  ccGGGATCCCGCGCTTTGCCGCCGCATCCAGATCAACCTGGTTTGTTCCGATACAGAAACAGCCAATAGCGACCAGTTTTTCTGCGGCGTTGATCACGTCTTCAGTCAGATGGGTACGGGATCGCAGGCCGATGAAGTGGGCATcgcg                                                                                                                 >  1:344953/1‑148 (MQ=255)
                                    gggATCCCGCGCTTTGCCGCCGCATCCAGATCAACCTGGTTTGTTCCGATACAGAAACAGCCAATAGCGACCAGTTTTTCTGCGGCGTTGATCACGTCTTCAGTCAGATGGGTACGGGATCGCAGGCCGATGAAGTGGGCATCGCggat                                                                                                              <  2:161015/149‑1 (MQ=255)
                                    gggATCCCGCGCTTTGCCGCCGCATCCAGATCAACCTGGTTTGTTCCGATACAGAAACAGCCAATAGCGACCAGTTTTTCTGCGGCGTTGATCACGTCTTCAGTCAGATGGGTACGGGATCGCAGGCCGATGAAGTGGGCATCGCggat                                                                                                              <  3:368902/149‑1 (MQ=255)
                                    gggATCCCGCGCTTTGCCGCCGCATCCAGATCAACCTGGTTTGTTCCGATACAGAAACAGCCAATAGCGACCAGTTTTTCTGCGGCGTTGATCACGTCTTCAGACAGATGGGTACGGGATCGCAGGCCGATGAAGTGGGCATCGCggat                                                                                                              <  2:66495/149‑1 (MQ=255)
                                     ggATCCCGCGCTTTGCCGCCGCATCCAGATCAACCTGGTTTGTTCCGATACAGAAACAGCCAATAGCGACCAGTTTTTCTGCGGCGTTGATCACGTCTTCAGTCAGATGGGTACGGGATCGCAGGCCGa                                                                                                                                 >  6:275248/1‑129 (MQ=255)
                                     ggATCCCGCGCTTTGCCGCCGCATCCAGATCAACCTGGTTTGTTCCGATACAGAAACAGCCAATAGCGACCAGTTTTTCTGCGGCGTTGATCACGTCTTCAGTCAGATGGGTACGGGATCGCAGGCCGa                                                                                                                                 <  5:275248/129‑1 (MQ=255)
                                      gATCCCGCGCTTTGCCGCCGCATCCAGATCAACCTGGTTTGTTCCGATACAGAAACAGCCAATAGCGACCAGTTTTTCTGCGGCGTTGATCACGTCTTCAGTCAGATGGGTACGGGATCGCAGGCCGATGAAGTGGGCATCGCggatgg                                                                                                            <  3:198423/149‑1 (MQ=255)
                                           cgcgCTTTGCCGCCGCATCCAGATCAACCTGGTTTGTTCCGATACAGAAACAGCCAATAGCGACCAGTTTTTCTGCGGCGTTGATCACGTCTTCAGTCAGATGGGTACGGGATCGCAGGCCGATGAAGTGGGCATCGCGGATGGATTCt                                                                                                       >  4:303532/1‑149 (MQ=255)
                                           cgcgCTTTGCCGCCGCATCCAGATCAACCTGGTTTGTTCCGATACAGAAACAGCCAATAGCGACCAGTTTTTCTGCGGCGTTGATCACGTCTTCAGTCAGATGGGTACGGGATCGCAGGCCGATGAAGTGGGCATCGCGGATGGATTCt                                                                                                       <  2:344953/149‑1 (MQ=255)
                                           cgcgCTTTGCCGCCGCATCCAGATCAACCTGGTTTGTTCCGATACAGAAACAGCCAATAGCGACCAGTTTTTCTGCGGCGTTGATCACGTCTTCAGTCAGATGGGTACGGGATCGCAGGCCGATGAAGTGGGCATCGCGGATGGATTCt                                                                                                       >  6:237690/1‑149 (MQ=255)
                                            gcgcTTTGCCGCCGCATCCAGATCAACCTGGTTTGTTCCGATACAGAAACAGCCAATAGCGACCAGTTTTTCTGCGGCGTTGATCACGTCTTCAGTCAGATGGGTACGGGATCGCAGGCCGATGAAGTGGGCATCGCGGATGGATTCtt                                                                                                      <  3:303532/149‑1 (MQ=255)
                                              gcTTTGCCGCCGCATCCAGATCAACCTGGATTTTTACGATACCGAAACAGCCAATAGCGACCAGTTTTTCTGCGGCGTTGATCACGTCGTCAGTCAGATGGGTACGGGATCGCAGTCCGATGAAGTGGGCATCGCGGATGGATTCtttt                                                                                                    <  2:96256/149‑1 (MQ=255)
                                                tttGCCGCCGCATCCAGATCAACCTGGTTTGTTCCGATACAGAAACAGCCAATAGCGACCAGTTTTTCTGCGGCGTTGATCACGTCTTCAGTCAGATGGGTACGGGATCGCAGGCCGATGAAGTGGGCATCGCGGATGGATTCTTTTaa                                                                                                  >  8:382492/1‑149 (MQ=255)
                                                   gccgccGCATCCAGATCAACCTGGTTTGTTCCGATACAGAAACAGCCAATAGCGACCAGTTTTTCTGCGGCGTTGATCACGTCTTCAGTCAGATGGGTACGGGATCGCAGGCCGATGCAGTGGGCATCGCGGATGGATTCTTTTAATTg                                                                                               >  1:137487/1‑149 (MQ=255)
                                                    ccgccgCATCCAGATCAACCTGGTTTGTTCCGATACAGAAACAGCCAATAGCGACCAGTTTTTCTGCGGCGTTGATCACGTCTTCAGTCAGATGGGTACGGGATCGCAGGCCGATGAAGTGGGCATCGCGGATGGATTCTTTTAATTGt                                                                                              <  3:151/149‑1 (MQ=255)
                                                     cgccgcATCCAGATCAACCTGGTTTGTTCCGATACAGAAACAGCCAATAGCGACCAGTTTTTCTGCGGCGTTGATCACGTCTTCAGTCAGATGGGTACGGGATCGCAGGCCGATGAAGTGGGCATCGCGGATGGATTCTTTTAATTGtt                                                                                             <  4:173147/149‑1 (MQ=255)
                                                     cgccgcATCCAGATCAACCTGGTTTGTTCCGATACAGAAACAGCCAATAGCGACCAGTTTTTCTGCGGCGTTGATCACGTCTTCAGTCAGATGGGTACGGGATCGCAGGCCGATGAAGTGGGCATCGCGGATGGATTCTTTTAATTGtt                                                                                             <  5:237690/149‑1 (MQ=255)
                                                      gccgcATCCAGATCAACCTGGTTTGTTCCGATACAGAAACAGCCAATAGCGACCAGTTTTTCTGCGGCGTTGATCACGTCTTCAGTCAGATGGGTACGGGATCGCAGGCCGATGAAGTGGGCATCGCGGATGGATTCTTTTAATTGTtc                                                                                            <  8:57610/149‑1 (MQ=255)
                                                       ccgcATCCAGATCAACCTGGTTTGTTCCGATACAGAAACAGCCAATAGCGACCAGTTTTTCTGCGGCGTTGATCACGTCTTCAGTCAGATGAGTACGGGATCGCAGGCCGATGAAGTg                                                                                                                          <  4:237662/118‑1 (MQ=255)
                                                       ccgcATCCAGATCAACCTGGTTTGTTCCGATACAGAAACAGCCAATAGCGACCAGTTTTTCTGCGGCGTTGATCACGTCTTCAGTCAGATGAGTACGGGATCGCAGGCCGATGAAGTg                                                                                                                          >  3:21701/1‑118 (MQ=255)
                                                       ccgcATCCAGATCAACCTGGTTTGTTCCGATACAGAAACAGCCAATAGCGACCAGTTTTTCTGCGGCGTTGATCACGTCTTCAGTCAGATGAGTACGGGATCGCAGGCCGATGAAGTg                                                                                                                          <  4:21701/118‑1 (MQ=255)
                                                       ccgcATCCAGATCAACCTGGTTTGTTCCGATACAGAAACAGCCAATAGCGACCAGTTTTTCTGCGGCGTTGATCACGTCTTCAGTCAGATGAGTACGGGATCGCAGGCCGATGAAGTg                                                                                                                          >  3:237662/1‑118 (MQ=255)
                                                          cATCCAGATCAACCTGGTTTGTTCCGATACAGAAACAGCCAATAGCGACCAGTTTTTCTGCGGCGTTGATCACGTCTTCAGTCAGATGGGTACGGGATCGCAGGCCGATGAAGTGGGCATCGCGGATGGATTCTTTTAATTGTtcatca                                                                                        >  4:126274/1‑149 (MQ=255)
                                                            tCCAGATCAACCTGGTTTGTTCCGATACAGAAACAGCCAATAGCGACCAGTTTTTCTGCGGCGTTGATCACGTCTTCAGTCAGATGGGTACGGGATCGCAGGCCGATGAAGTGGGCATCGCggat                                                                                                              >  7:254108/1‑125 (MQ=255)
                                                            tCCAGATCAACCTGGTTTGTTCCGATACAGAAACAGCCAATAGCGACCAGTTTTTCTGCGGCGTTGATCACGTCTTCAGTCAGATGGGTACGGGATCGCAGGCCGATGAAGTGGGCATCGCggat                                                                                                              <  8:254108/125‑1 (MQ=255)
                                                                 aTCAACCTGGTTTGTTCCGATACAGAAACAGCCAATAGCGACCAGTTTTTCTGCGGCGTTGATCACGTCTTCAGTCAGATGGGTACGGGATCGCAGGCCGATGAAGTGGGCATCGCGGATGGATTCTTTTAATTGTTCATCATCCAgcg                                                                                 >  6:204663/1‑149 (MQ=255)
                                                                      ccTGGTTTGTTCCGATACAGAAACAGCCAATAGCGACCAGTTTTTCTGCGGCGTTGATCACGTCTTCAGTCAGATGGGTACGGGATCGCAGGCCGATGAAGTGGGCATCGCGGATGGATTCTTTTAATTGTTCATCATCCAGCGCGCCt                                                                            >  6:173667/1‑149 (MQ=255)
                                                                       cTGGTTTGTTCCGATACAGAAACAGCCAATAGCGACCAGTTTTTCTGCGGCGTTGATCACGTCTTCAGTCAGATGGGTACGGGATCGCAGGCCGATGAAGTGGGCATCGCGGATGGATTCtt                                                                                                      >  2:218808/1‑122 (MQ=255)
                                                                       cTGGTTTGTTCCGATACAGAAACAGCCAATAGCGACCAGTTTTTCTGCGGCGTTGATCACGTCTTCAGTCAGATGGGTACGGGATCGCAGGCCGATGAAGTGGGCATCGCGGATGGATTCtt                                                                                                      <  1:218808/122‑1 (MQ=255)
                                                                          gTTTGTTCCGATACAGAAACAGCCAATAGCGACCAGTTTTTCTGCGGCGTTGATCACGTCTTCAGTCAGATGGGTACGGGATCGCAGGCCGATGAAGTGGGCATCGCGGATGGATTCTTTTAATTGTTCATCATCCAGCGCGCCTTtg                                                                         <  5:246585/148‑1 (MQ=255)
                                                                               ttCCGATACAGAAACAGCCAATAGCGACCAGTTTTTCTGCGGCGTGGATCACGTCTTCACTCAGATGGGTACGGGATCGCAGGCCGATGAAGTGGGCATCGCGGATGGATTCTTTTAATTGTGCATCATCCAGCGCGCCTTTGTGaaag                                                                   <  5:267713/149‑2 (MQ=255)
                                                                                 ccGATACAGAAACAGCCAATAGCGACCAGTTTTTCTGCGGCGTTGATCACGTCTTCAGTCAGATGGGTACGGGATCGCAGGCCGATGAAGTGGGCATCGCGGATGGATTCTTTTAATTGTTCATCATCCAGCGCGCCTTTGTGAAATTc                                                                 >  8:361787/1‑149 (MQ=255)
                                                                                      aCAGAAACAGCCAATAGCGACCAGTTTTTCTGCGGCGTTGATCACGTCTTCAGTCAGATGGGTACGGGATCGCAGGCCGATGAAGTGGGCATCGCGGATGGATTCTTTTAATTGTTCATCATCCAGCGCGCCTTTGTGAAATTCGATGt                                                            >  7:238196/1‑149 (MQ=255)
                                                                                           aaCATCCAATATCGACCAGTTTTTCTGCGGGGTTGATGACGTCTTCAGGCAGATGGGTAAGGTATCGCATGACGATGAAGTGTGCATCGCGGATGGATTCTTTTAATTGTTCATCATCCAGCGCGCCTTTGTGAAATTCGCTGTTGgt                                                        <  2:137487/148‑1 (MQ=255)
                                                                                           aaCAGCCAATAGCGACCAGTTTTTCTGCGGCGTTGATCACGTCTTCAGTCAGATGGGTACGGGATCGCAGGCCGATGAAGTGGGCATCGCGGATGGATTCTTTTAATTGTTCATCATCCAGCGCGCCTTTGTGAAATTCGa                                                               >  8:156623/1‑141 (MQ=255)
                                                                                           aaCAGCCAATAGCGACCAGTTTTTCTGCGGCGTTGATCACGTCTTCAGTCAGATGGGTACGGGATCGCAGGCCGATGAAGTGGGCATCGCGGATGGATTCTTTTAATTGTTCATCATCCAGCGCGCCTTTGTGAAATTCGa                                                               <  7:156623/141‑1 (MQ=255)
                                                                                              aGCCAATAGCGACCAGTTTTTCTGCGGCGTTGATCACGTCTTCAGTCAGATGGGTACGGGATCGCAGGCCGATGAAGTGGGCATCGCGGATGGATTCTTTTAATTGTTCATCATCCAGCGCGCCTTTGTGAAATTCGATGTTGGTGTaa                                                    >  4:166117/1‑149 (MQ=255)
                                                                                                 cAATAGCGACCAGTTTTTCTGCGGCGTTGATCACGTCTTCAGTCAGATGGGTACGGGATCGCAGGCCGATGAAGTGGGCATCGCGGATGGATTCTTTTAATTGTTCATCATCCAGCGCGCCTTtgtg                                                                       >  3:128807/1‑127 (MQ=255)
                                                                                                 cAATAGCGACCAGTTTTTCTGCGGCGTTGATCACGTCTTCAGTCAGATGGGTACGGGATCGCAGGCCGATGAAGTGGGCATCGCGGATGGATTCTTTTAATTGTTCATCATCCAGCGCGCCTTtgtg                                                                       <  4:128807/127‑1 (MQ=255)
                                                                                                    tAGCGACCAGTTTTTCTGCGGCGTTGATCACGTCTTCAGTCAGATGGGTACGGGATCGCAGGCCGATGAAGTGGGCATCGCGGATGGATTCTTTTAATTGTTCATCATCCAGCGCGCCTTTGTGAAATTCGATGTTGGTGTAACCAGCt                                              <  5:173667/149‑1 (MQ=255)
                                                                                                      gCGACCAGTTTTTCTGCGGCGTTGATCACGTCTTCAGTCAGATGGGTACGGGATCGCAGGCCGATGAAGTGGGCATCGCGGATGGATTCTTTTAATTGTTCATCATCCAGCGCGCCTTTGTGAAATTCGATGTTGGTGTAACCAGCTGc                                            <  3:286257/149‑1 (MQ=255)
                                                                                                              tttttCTGCGGCGTTGATCACGTCTTCAGTCAGATGGGTACGGGATCGCAGGCCGATGAAGTGGGCATCGCGGATGGATTCTTTTAATTGTTCATCATCCAGCGCGCCTTTGTGAAATTCGATGTTGGTGTAACCAGCTGCACGAAGGc                                    >  5:334576/1‑149 (MQ=255)
                                                                                                              tttttCTGCGGCGTTGATCACGTCTTCAGTCAGATGGGTACGGGATCGCAGGCCGATGAAGTGGGCATCGCGGATGGATTCTTTTAATTGTTCATCATCCAGCGCGCCTTTGTGAAATTCGATGTTGGTGTAACCAGCTGCACGAAGGc                                    >  1:43922/1‑149 (MQ=255)
                                                                                                                   cTGCGGCGTTGATCACGTCTTCAGTCAGATGGGTACGGGATCGCAGGCCGATGAAGTGGGCATCGCGGATGGATTCTTTTAATTGTTCATCATCCAGCGCGCCTTTGTGAAATTCGATGTTGGTGTAACCAGCTGCACGAAGGCTTTc                                <  7:314834/148‑1 (MQ=255)
                                                                                                                   cTGCGGCGTTGATCACGTCTTCAGTCAGATGGGTACGGGATCGCAGGCCGATGAAGTGGGCATCGCGGATGGAGTCTTTTAATTGTTCATCATCCAGCGCGCCTTTGTGAAATTCGATGTTGGTGTAACCAGCTGCACGAAGGCTTTc                                <  5:54506/148‑1 (MQ=255)
                                                                                                                    tGCGGCGTTGATCACGTCTTCAGTCAGATGGGTACGGGATCGCAGGCCGATGAAGTGGGCATCGCGGATGGATTCTTTTAATTGTTCATCATCCAGCGCGCCTTTGTGAAATTCGATGTTGGTGTAACCAGCTGCACGAAGGCTTTc                                >  2:194532/1‑147 (MQ=255)
                                                                                                                         cggagATCACGTCTTCAGTCAGATGCGTACGGGATCGCAGGCCGATGAAGTGGGCATCGCGGATGGATTCTTTTAATTGTTCCTCATCCAGCGCGCCTTTGTGAAATTCGATGTTGGTGTAACCAGCTGCACGAAGGCTTTCCAGCGcc                         <  2:361643/145‑1 (MQ=255)
                                                                                                                               tCACGTCTTCAGTCAGATGGGTACGGGATCGCAGGCCGATGAAGTGGGCATCGCGGATGGATTCTTTTAATTGTTCATCATCCAGCGCGCCTTTGTGAAATTCGATGTTGGTGTAACCAGCTGCACGaa                                       >  1:95282/1‑129 (MQ=255)
                                                                                                                               tCACGTCTTCAGTCAGATGGGTACGGGATCGCAGGCCGATGAAGTGGGCATCGCGGATGGATTCTTTTAATTGTTCATCATCCAGCGCGCCTTTGTGAAATTCGATGTTGGTGTAACCAGCTGCACGaa                                       <  2:95282/129‑1 (MQ=255)
                                                                                                                                   gTCTTCAGTCAGATGGGTACGGGATCGCAGGCCGATGAAGTGGGCATCGCGGATGGATTCTTTTAATTGTTCATCATCCAGCGCGCCTTTGTGAAATTCGATGTTGGTGTAACCAGCTGCACGAAGGCTTTCCAGCGCCTTTTGGTGCa               <  7:226078/149‑1 (MQ=255)
                                                                                                                                        cagtcagATGGGTACGGGATCTCAGGCCGATGAAGTGGGCATCGAGGATGGATTATTTTAATTGTTAATCAACCAGCGCGCCTTTGTGAAATTCGATTTTGGTGTAACCAGCTGCACGAAGGCTTTACATCGACTTTTGGTGCAAGcat          >  4:291279/1‑147 (MQ=255)
                                                                                                                                         agtcagATGGGTACGGGATCGCAGGCCGATGAAGTGGGCATCGCGGATGGATTCTTTTAATTGTTCATCATCCAGCGCGCCTTTGTGAAATTCGATGTTGGTGTAACCAGCTGCACGAAGGCTTTCCAGCGCCTTTTGGTGCACGCCtt         <  5:204663/149‑1 (MQ=255)
                                                                                                                                         agtcagATGGGTACGGGATCGCAGGCCGATGAAGTGGGCATCGCGGATGGATTCTTTTAATTGTTCATCATCCAGCGCGCCTTTGTGAAATTCGATGTTGGTGTAACCAGCTGCACGAAGGCTTTCCAGCGCCTTTTGGTGCACGCCtt         >  1:1735/1‑149 (MQ=255)
                                                                                                                                           tcagATGGGTACGGGATCGCAGGCCGATGAAGTGGGCATCGCGGATGGATTCTTTTAATTGTTCATCATCCAGCGCGCCTTTGTGAAATTCGATGTTGGTGTAACCAGCTGCACGAAGGCTTTCCAGCGCCTTTTGGTGCACGCCTTc        >  4:357042/1‑148 (MQ=255)
                                                                                                                                           tcagATGGGTACGGGATCGCAGGCCGATGAAGTGGGCATCGCGGATGGATTCTTTTAATTGTTCATCATCCAGCGCGCCTTTGTGAAATTCGATGTTGGTGTAACCAGCTGCACGAAGGCTTTCCAGCGCCTTTTGGTGCACGCCTTCt       <  7:386677/149‑1 (MQ=255)
                                                                                                                                             agATGGGTACGGGATCGCAGGCCGATGAAGTGGGCATCGCGGATGGATTCTTTTAATTGTTCATCATCCAGCGCGCCTTTGTGAAATTCGATGTTGGTGTAACCAGCTGCACGAAGGCTTTCCAGCGCCTTTTGGTGCACGCCTTCTAc     <  2:43922/149‑1 (MQ=255)
                                                                                                                                               aTGGGTACGGGATCGCAGGCCGATGAAGTGGGCATCGCGGATGGATTCTTTTAATTGTTCATCATCCAGCGCGCCTTTGTGAAATTCGATGTTGGTGTAACCAGCTGCACGAAGGCTTTCCAGCGCCTTTTGGTGCACGCCTTCTAc     >  5:293975/1‑147 (MQ=255)
                                                                                                                                                tGGGTACGGGATCGCAGGCCGATGAAGTGGGCATCGCGGATGGATTCTTTTAATTGTTCATCATCCAGCGCGCCTTTGTGAAATTCGATGTTGGTGTAACCAGCTGCACGAAGGCTTTCCAGCGCCTTTTGGTGCACGCCTTCTACcag  <  2:401739/149‑1 (MQ=255)
                                                                                                                                                 gggTACGGGATCGCAGGCCGATGAAGTGGGCATCGCGGATGGATTCTTTTAATTGTTCATCATCCAGCGCGCCTTTGTGAAATTCGATGTTGGTGTAACCAGCTGCACGAAGGCTTTCCAGCGCCTTTTGGTGCACGCCTTCTACcag  >  3:4724/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
ACAGAGCGCGTATTTGAGAACGGTGCGTTAAATACCGGGATCCCGCGCTTTGCCGCCGCATCCAGATCAACCTGGTTTGTTCCGATACAGAAACAGCCAATAGCGACCAGTTTTTCTGCGGCGTTGATCACGTCTTCAGTCAGATGGGTACGGGATCGCAGGCCGATGAAGTGGGCATCGCGGATGGATTCTTTTAATTGTTCATCATCCAGCGCGCCTTTGTGAAATTCGATGTTGGTGTAACCAGCTGCACGAAGGCTTTCCAGCGCCTTTTGGTGCACGCCTTCTACCAG  >  NC_000913/3058076‑3058368

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