Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,260,724 G→A 5.8% D565D (GAC→GAT tdcE ← pyruvate formate‑lyase 4/2‑ketobutyrate formate‑lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,260,7240GA 5.8% 218.9 / 2.4 69D565D (GAC→GATtdcEpyruvate formate‑lyase 4/2‑ketobutyrate formate‑lyase
Reads supporting (aligned to +/- strand):  ref base G (33/32);  new base A (2/2);  total (35/34)
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.72e-01

TTAGAAGTGATAGTCAGAATCGACTGGGTAGGGACGGCGTTGCGATAGGTTGGCAGCGCTTTAATTTTCTTCATAAAGCGTTCAACCAGGTCGCAGGCAATGCTGTCTACGCGCTCGTCGTTGTTGCCGTACTGCGGATATTCACCGTCGATTTCAAAGTCCACCGCCAGGCCGTTTTCGTCACGGATTGGTTTCACGCGGGCATATTTGATGGCAGACAGGGAGTCCGTCGCCACCGACAGGCCCGCGATGCCGCATGCCATAGTGCGATAGACATCACGATCGTGCAGCGC  >  NC_000913/3260578‑3260870
                                                                                                                                                  |                                                                                                                                                  
ttAGAAGTGATAGTCAGAATCGACTGGGTAGGGACGGCGTTGCGATAGGTTGGCAGCGCTTTAATTTTCTTCATAAAGCGTTCAACCAGGTCGCAGGCAATGCTGTCTACGCGCTCGTCGTTGTTGCCGTACTGCGGATATTCACCGTc                                                                                                                                                  <  3:95810/149‑1 (MQ=255)
   gAAGTGATAGTCAGAATCGACTGGGTAGGGACGGCGTTGCGATAGGTTGGCAGCGCTTTAATTTTCTTCATAAAGCGTTCAACCAGGTCGCAGGCAATGCTGTCTACGCGCTCGTCGTTGTTGCCGTACTGCGGATATTCACCGTCGa                                                                                                                                                <  8:251176/148‑1 (MQ=255)
      gTGATAGTCAGAATCGACTGGGTAGGGACGGCGTTGCGATAGGTTGGCAGCGCTTTAATTTTCTTCATAAAGCGTTCAACCAGGTCGCAGGCAATGCTGTCTACGCGCTCGTCGTTGTTGCCGTACTGCGGATATTCACCGTCGAttt                                                                                                                                             >  5:129400/1‑148 (MQ=255)
       tGATAGTCATAATCGACTGTGTAGTGTCGGCGTTGCGATAGGTTGGCAGCGCTTTAATTTGCTTCATATATCGTTCAACCAGGTCGCAGGCAATGCTGTCTACGCGCTGGGCGTTGTTGCCGTACTGCGGATATTCACCGTCGATTTCa                                                                                                                                           <  7:122735/149‑1 (MQ=255)
          tAGTCAGAATCGACTGGGTAGGGACGGCGTTGCGATAGGTTGGCAGCGCTTTAATTTTCTTCATAAAGCGTTCAACCAGGTCGCAGGCAATGCTGTCTACGCGCTCGTCGTTGTTGCCGTACTGCGGATATTCACCGTCGATTTCaaa                                                                                                                                         <  8:85874/148‑1 (MQ=255)
            gTCAGAATCGACTGGGTAGGGACGGCGTTGCGATAGGTTGGCAGCGCTTTAATTTTCTTCATAAAGCGTTCAACCAGGTCGCAGGCAATGCTGTCTACGCGCTCGTCGTTGTTGCCGTACTGCGGATATTCACCGTCGATTTCAAAGTc                                                                                                                                      >  3:198607/1‑149 (MQ=255)
            gTCAGAATCGACTGGGTAGGGACGGCGTTGCGATAGGTTGGCAGCGCTTTAATTTTCTTCATAAAGCGTTCAACCAGGTCGCAGGCAATGCTGTCTACGCGCTCGTCGTTGTTGCCGTACTGCGGATATTCACCGTCGATTTCAAAGTc                                                                                                                                      <  5:23937/149‑1 (MQ=255)
              cAGAATCGACTGGGTAGGGACGGCGTTGCGATAGGTTGGCAGCGCTTTAATTTTCTTCATAAAGCGTTCAACCAGGTCGCAGGCAATGCTGTCTACGCGCTCGTCGTTGTTGCCGTACTGCGGATATTCACCGTCGATTTCAAAGTCCa                                                                                                                                    >  3:51783/1‑149 (MQ=255)
               aGAATCGACTGGGTAGGGACGGCGTTGCGATAGGTTGGCAGCGCTTTAATTTTCTTCATAAAGCGTTCAACCAGGTCGCAGGCAATGCTGTCTACGCGCTCGTCGTTGTTGCCGTACTGCGGATATTCACCGTCGATATCAAAGTCCAc                                                                                                                                   >  5:70172/1‑149 (MQ=255)
                     gACTGGGTAGGGACGGCGTTGCGATAGGTTGGCAGCGCTTTAATTTTCTTCATAAAGCGTTCAACCAGGTCGCAGGCAATGCTGTCTACGCGCTCGTCGTTGTTGCCGTACTGCGGATATTCACCGTCGATTTCAAAGTCCACCGCCAg                                                                                                                             <  2:15688/149‑1 (MQ=255)
                       cTGGGTAGGGACGGCGTTGCGATAGGTTGGCAGCGCTTTAATTTTCTTCATAAAGCGTTCAACCAGGTCGCAGGCAATGCTGTCTACGCGCTCGTCGTTGTTGCCGTACTGCGGATATTCACCGTCGATTTCAAAGTCCACCGCCAGGc                                                                                                                           >  3:91799/1‑149 (MQ=255)
                        tGGGTAGGGACGGCGTTGCGATAGGTTGGCAGCGCTTTAATTTTCTTCATAAAGCGTTCAACCAGGTCGCAGGCAATGCTGTCTACGCGCTCGTCGTTGTTGCCGTACTGCGGATATTCACCGTCGATTTCAAAGTCCACCGCCAGGcc                                                                                                                          >  7:230/1‑149 (MQ=255)
                           gTAGGGACGGCGTTGCGATAGGTTGGCAGCGCTTTAATTTTCTTCATAAAGCGTTCAACCAGGTCGCAGGCAATGCTGTCTACGCGCTCGTCGTTGTTGCCGTACTGCGGATATTCACCGTCGATTTCAAAGTCCACCGCCAGGCCGtt                                                                                                                       <  4:187917/149‑1 (MQ=255)
                              gggACGGCGTTGCGATAGGTTGGCAGCGCTTTAATTTTCTTCATAAAGCGTTCAACCAGGTCGCAGGCAATGCTGTCTACGCGCTCGTCGTTGTTGCCGTACTGCGGATATTCACCGTCGATTTCAAAGTCCACCGCCAGGCCGtttt                                                                                                                     <  8:230/148‑1 (MQ=255)
                              gggACGGCGTTGCGATAGGTTGGCAGCGCTTTAATTTTCTTCATAAAGCGTTCAACCAGGTCGCAGGCAATGCTGTCTACGCGCTCGTCGTTGTTGCCGTACTGCGGATATTCACCGTCGATTTCAAAGTCCACCGCCAGGCCGTTTTc                                                                                                                    <  4:51783/149‑1 (MQ=255)
                               ggACGGCGTTGCGATAGGTTGGCAGCGCTTTAATTTTCTTCATAAAGCGTTCAACCAGGTCGCAGGCAATGCTGTCTACGCGCTCGTCGTTGTTGCCGTACTGCGGATATTCACCGTCGATTTCAAAGTCCACCGCCAGGCCGTTTTCg                                                                                                                   >  4:154231/1‑149 (MQ=255)
                                   ggCGTTGCGATAGGTTGGCAGCGCTTTAATTTTCTTCATAAAGCGTTCAACCAGGTCGCAGGCAATGCTGTCTACGCGCTCGTCGTTGTTGCCGTACTGCGGATATTCACCGTCGATTTCAAAGTCCACCGCCAGGCCGTTTTCGTCAc                                                                                                               <  2:183734/149‑1 (MQ=255)
                                      gTTGCGATAGGTTGGCAGCGCTTTAATTTTCTTCATAAAGCGTTCAACCAGGTCGCAGGCAATGCTGTCTACGCGCTCGTCGTTGTTGCCGTACTGCGGATATTCACCGTCGATTTCAAAGTCCACCGCCAGGCCGTTTTCGTCACGGa                                                                                                            >  2:116419/1‑149 (MQ=255)
                                      gTTGCGATAGGTTGGCAGCGCTTTAATTTTCTTCATAAAGCGTTCAACCAGGTCGCAGGCAATGCTGTCTACGCGCTCGTCGTTGTTGCCGTACTGCGGATATTCACCGTCGATTTCAAAGTCCACCGCCAGGCCGTTTTCGTCACGGa                                                                                                            <  5:217075/149‑1 (MQ=255)
                                         gCGATTGGTTGGCAGCGCTTTAATTTTCTTCATTAAGCGTTGAACCAGGTCGCAGGCAATGCTGTCTACGCGCTCGGCGTTGTTGCCGTACTGCGGATATTCACCGTCGATTTCAAAGTCCACCGCCAGGACGTTTTCGTCACGGAtt                                                                                                          <  6:78624/148‑1 (MQ=255)
                                          cGATAGGTTGGCAGCGCTTTAATTTTCTTCATAAAGCGTTCAACCAGGTCGCAGGCAATGCTGTCTACGCGCTCGTCGTTGTTGCCGTACTGCGGATATTCACCGTCGATTTCAAAGTCCACCGCCAGGCCGTTTTCGTCACGGATTgg                                                                                                        >  8:159395/1‑149 (MQ=255)
                                           gATAGGTTGGCAGCGCTTTAATTTTCTTCATAAAGCGTTCAACCAGGTCGCAGGCAATGCTGTCTACGCGCTCGTCGTTGTTGCCGTACTGCGGATATTCACCGTCGATTTCAAAGTCCACCGCCAGGCCGTTTTCGTCACGGATTGGt                                                                                                       <  8:208111/149‑1 (MQ=255)
                                                gTTGGCAGCGCTTTAATTTTCTTCATAAAGCGTTCAACCAGGTCGCAGGCACTGCTGTCTCCGCGCTCGTCGTTGTTGCCGTACTGCGGATATTCACCGTCGATTTCAAAGTCCACCGCCAGGCCGTTTTCGTCCCGGATTGGTTTCAc                                                                                                  >  8:192889/1‑149 (MQ=255)
                                                    ccgtctcTTTAATTTTCTTCATAAAGCGTTCAACCAGGTCGCAGGCAATGCTGTCTACGCGCTCGTCGTTGTTGCCGTACTGCGGATATTCACCGTCGATTTCAAAGTCCACCGCCAGGCCGTTTTCGTCACGGATTGGTTTCAcgcg                                                                                               <  5:44698/142‑1 (MQ=255)
                                                         gcTTTAATTTTCTTCATAAAGCGTTCAACCAGGTCGCAGGCAATGCTGTCTACGCGCTCGTCGTTGTTGCCGTACTGCGGATATTCACCGTCGATTTCAAAGTCCACCGCCAGGCCGTTTTCGTCACGGATTGGTTTCACGCGGGCata                                                                                         >  1:86821/1‑149 (MQ=255)
                                                             tAATTTTCTTCATAAAGCGTTCAACCAGGTCGCAGGCAATGCTGTCTACGCGCTCGTCGTTGTTGCCGTACTGCGGATATTCACCGTCGATTTCAAAGTCCACCGCCAGGCCGTTTTCGTCACGGATTGGTTTCACGCGGGCATAttt                                                                                      >  6:57611/1‑148 (MQ=255)
                                                               aTTTTCTTCATAAAGCGTTCAACCAGGTCGCAGGCAATGCTGTCTACGCGCTCGTCGTTGTTGCCGTACTGCGGATATTCACCGTCGATTTCAAAGTCCACCGCCAGGCCGTTTTCGTCACGGATTGGTTTCACGCGGGCATATTTGa                                                                                    <  7:159395/148‑1 (MQ=255)
                                                                     ttcATAAAGCGTTCAACCAGGTCGCAGGCAATGCTGTCTACGCGCTCGTCGTTGTTGCCGTACTGCGGATATTCACCGTCGATTTCAAAGTCCACCGCCAGGCCGTTTTCGTCACGGATTGGTTTCACGCGGGCATATTTGATGGCAGa                                                                             >  7:139455/1‑149 (MQ=255)
                                                                       cATAAAGCGTTCAACCAGGTCGCAGGCAATGCTGTCTACGCGCTCGTAGTTGTTGCCGTACTGCGGATATTCACCGTCGATTTCAAAGACCACCGCCAGGCCGTTTTCGTCACGGATTGGTTTCACGCGGGCATATTTGATGGCAGa                                                                             >  4:156805/1‑147 (MQ=255)
                                                                        aTAAAGCGTTCAACCAGGTCGCAGGCAATGCTGTCTACGCGCTCGTCGTTGTTGCCGTACTGCGGATATTCACCGTCGATTTCAAAGTCCACCGCCAGGCCGTTTTCGTCACGGATTGGTTTCACGCGGGCATATTTGATGGCAGACAg                                                                          <  1:116419/149‑1 (MQ=255)
                                                                        aTAAAGCGTTCAACCAGGTCGCAGGCAATGCTGTCTACGCGCTCGTCGTTGTTGCCGTACTGCGGATATTCACCGTCGATTTCAAAGTCCACCGCCAGGCCGTTTTCGTCACGGATTGGTTTCACGCGGGCATATTTGATGGCAGACAg                                                                          <  6:111265/149‑1 (MQ=255)
                                                                         tAAAGCGTTCAACCAGGTCGCAGGCAATGCTGTCTACGCGCTCGTCGTTGTTGCCGTACTGCGGATATTCACCGTCGATTTCAAAGTCCACCGCCAGGCCGTTTTCGTCACGGATTGGTTTCACGCGGGCATATTTGATGGCAGACAgg                                                                         >  2:83643/1‑149 (MQ=255)
                                                                         tAAAGCGTTCAACCAGGTCGCAGGCAATGCTGTCTACGCGCTCGTCGTTGTTGCCGTACTGCGGATATTCACCGTCGATTTCAAAGTCCACCGCCAGGCCGTTTTCGTCACGGATTGGTTTCACGCGGGCATATTTGATGGCAGACAgg                                                                         >  2:181879/1‑149 (MQ=255)
                                                                         tAAAGCGTTCAACCAGGTCGCAGGCAATGCTGTCTACGCGCTCGTCGTTGTTGCCGTACTGCGGATATTCACCGTCGATTTCAAAGTCCACCGCCAGGCCGTTTTCGTCACGGATTGGTTTCACGCGGGCATATTTGATGGCAGACAgg                                                                         <  6:129400/149‑1 (MQ=255)
                                                                          aaaGCGTTCAACCAGGTCGCAGGCAATGCTGTCTACGCGCTCGTAGTTGTTGCCTTACTGCGGATATTCACAGTCGATTTCAAAGTCCACAGAAAGGCCGTTTTCGTCACGGATTGTTTTCACGCGGGCATATTTTATGGCAGAAAgg                                                                         >  7:108479/1‑148 (MQ=255)
                                                                                ttCAACCAGGTCGCAGGCAATGATGTCTACGCGCTCGTCGTTGTTGCCGTACTGCGGATATTCACCGTCGATTTCAAAGTCCACCGCCAGGCCGTTTTCGTCACGGATTGGTTTCACGCGGGCATATTTGATGGCAGACAGGGAGTCCg                                                                  >  4:107148/1‑149 (MQ=255)
                                                                                   aaCCAGGTCGCAGGCAATGCTGTCTACGCGCTCGTCGTTGTTGCCGTACTGCGGATATTCACCATCGATTTCAAAGTCCACCGCCAGGCCGTTTTCGTCACGGATTGGTTTCACGCGGGCATATTTGATGGCAGACAGGGa                                                                       >  7:183017/1‑141 (MQ=255)
                                                                                   aaCCAGGTCGCAGGCAATGCTGTCTACGCGCTCGTCGTTGTTGCCGTACTGCGGATATTCACCATCGATTTCAAAGTCCACCGCCAGGCCGTTTTCGTCACGGATTGGTTTCACGCGGGCATATTTGATGGCAGACAGGGa                                                                       <  8:183017/141‑1 (MQ=255)
                                                                                   aaCCAGGTCGCAGGCAATGCTGTCTACGCGCTCGTCGTTGTTGCCGTACTGCGGATATTCACCATCGATTTCAAAGTCCACCGCCAGGCCGTTTTCGTCACGGATTGGTTTCACGCGGGCATATTTGATGGCAGACAGGGa                                                                       >  1:163321/1‑141 (MQ=255)
                                                                                   aaCCAGGTCGCAGGCAATGCTGTCTACGCGCTCGTCGTTGTTGCCGTACTGCGGATATTCACCATCGATTTCAAAGTCCACCGCCAGGCCGTTTTCGTCACGGATTGGTTTCACGCGGGCATATTTGATGGCAGACAGGGa                                                                       <  2:163321/141‑1 (MQ=255)
                                                                                      cAGGTCGCAGGCAATGCTGTCTACGCGCTCGTCGTTGTTGCCGTACTGCGGATATTCACCGTCGATTTCAAAGTCCACCGCCAGGCCGTTTTCGTCACGGATTGGTTTCACGCGGGCATATTTGATGGCAGAc                                                                            >  4:211599/1‑133 (MQ=255)
                                                                                      cAGGTCGCAGGCAATGCTGTCTACGCGCTCGTCGTTGTTGCCGTACTGCGGATATTCACCGTCGATTTCAAAGTCCACCGCCAGGCCGTTTTCGTCACGGATTGGTTTCACGCGGGCATATTTGATGGCAGAc                                                                            <  3:211599/133‑1 (MQ=255)
                                                                                         gTCGCAGGCAATGCTGTCTACGCACTCGTCGTTGTTGCCGTACTGCGGATATTCACCGTCGATTTCAAAGTCCACCGCCAGGCCGTTTTCGTCACGGATTGGTTTCACGCGGGCATATTTGATGGCAGACAGGGAGTCCGTCGCCAcc                                                          >  1:245529/1‑148 (MQ=255)
                                                                                                gCAATGCTGTCTACGCGCTCGTCGTTGTTGCCGTACTGCGGATATTCACCGTCGATTTCAAAGTCCACCGCCAGGCCGTTTTCGTCACGGATTGGTTTCAGGCGGGCATATTTGATGGCAGACAGGGAGTCCGTCGCCACCGACAGGcc                                                  <  4:198607/149‑1 (MQ=255)
                                                                                                gCAATGCTGTCTACGCGCTCGTCGTTGTTGCCGTACTGCGGATATTCACCGTCGATTTCAAAGTCCACCGCCAGGCCGTTTTCGTCACGGATTGGTTTCACGCGGGCATATTTGATGGCAGACAGGGa                                                                       >  5:128701/1‑128 (MQ=255)
                                                                                                gCAATGCTGTCTACGCGCTCGTCGTTGTTGCCGTACTGCGGATATTCACCGTCGATTTCAAAGTCCACCGCCAGGCCGTTTTCGTCACGGATTGGTTTCACGCGGGCATATTTGATGGCAGACAGGGa                                                                       <  6:128701/128‑1 (MQ=255)
                                                                                                  aaTGCTGTCTACGCGCTCGTCGTTGTTGCCGTACTGCGGATATTCACCGTCGATTTCAAAGTCCACCGCCAGGCCGTTTTCGTCACGGATTGGTTTCACGCGGGCATATTTGATGGCAGACAGGGAGTCCGTCGCCACCGACAGGCCcg                                                >  7:21130/1‑149 (MQ=255)
                                                                                                         tCTACGCGCTCGTCGTTGTTGCCGTACTGCGGATATTCACCGTCGATTTCAAAGTCCACCGCCAGGCCGTTTTCGTCACGGATTGGTTTCACGCGGGCATATTTGATGGCAGACAGGGAGTCCGTCGCCACCGACAGGCCCGCGATg                                           <  2:86821/147‑1 (MQ=255)
                                                                                                           tACGCGCTCGTCGTTGTTGCCGTACTGCGGATATTCACCGTCGATTTCAAAGTCCACCGCCAGGCCGTTTTCGTCACGGATTGGTTTCACGCGGGCATATTTGATGGTAGACAGGGa                                                                       >  5:198235/1‑117 (MQ=255)
                                                                                                           tACGCGCTCGTCGTTGTTGCCGTACTGCGGATATTCACCGTCGATTTCAAAGTCCACCGCCAGGCCGTTTTCGTCACGGATTGGTTTCACGCGGGCATATTTGATGGTAGACAGGGa                                                                       <  6:198235/117‑1 (MQ=255)
                                                                                                                gcTCGTCGTTGTTGCCGTACTGCGGATATTCACCGTCGATTTCAAAGTCCACCGCCAGGCCGTTTTCGTCACGGATTGGTTTCACGCGGGCATATTTGATGGCAGACAgg                                                                         >  3:132642/1‑110 (MQ=255)
                                                                                                                gcTCGTCGTTGTTGCCGTACTGCGGATATTCACCGTCGATTTCAAAGTCCACCGCCAGGCCGTTTTCGTCACGGATTGGTTTCACGCGGGCATATTTGATGGCAGACAgg                                                                         <  4:132642/110‑1 (MQ=255)
                                                                                                                    gtcgtTGTTGCCGTACTGCGGATATTCACCGTCGATTTCAAAGTCCACCGCCAGGCCGTTTTCGTCACGGATTGGTTTCACGCGGGCATATTTGATGGCAGCCAGGGAGTCCGTCGCCACCGACAGGCCCGCGATGCCGCATGCCATAg                              >  1:97280/1‑149 (MQ=255)
                                                                                                                        ttgttgCCGTACTGCGGATATTCACCGTCGATTTCAAAGTCCACCGCCAGGCCGTTTTCGTCACGGATTGGTTTCACGCGGGCATATTTGATGGCAGACAGGGAGTCCGTCGCCACCGACAGGCCCGCGATGCCGCATGCCATAGTGc                           >  3:131848/1‑148 (MQ=255)
                                                                                                                          gttgCCGTTCTGCTGATATTCACCGTCGATTTCAAATTCCACCGCGAGGCCGTTTTCGTCTCGGATTGGTTTCACGCGGGCATATTTGTTGGCAGACAGGGAGTCCGTCGCCACCGAGAGGCCCGCGATGcc                                         <  6:1199/132‑1 (MQ=255)
                                                                                                                          gttgCCGTACTGCGGATATTCACCGTCGATTTCAAAGTCCACCGCCAGGCCGTTTTCGTCACGGATTGGTTTCACGCGGGCATATTTGATGGCAGACAGGGAGTCCGTCGCCACCGACAGGCCCGCGATGCCGCATGCCATAGTGCGa                         <  4:131848/148‑1 (MQ=255)
                                                                                                                          cttgCCGTACTGCGGATATTCACCGTCGATTTCAAAGTCCACAGCCAGGCCGTTTTCGTCACGGATTGGTTTCACGCGGGAAAATTTGATGGCAGACATGGAGTCCGTCGCCACCGACAGGCCCGCGAAGcc                                         >  5:1199/2‑132 (MQ=255)
                                                                                                                           ttgCCGTACTGCGGATATTCACCGTCGATTTCAAAGTCCACCGCCAGGCCGTTTTCGTCACGGATTGGTTTCACGCGGGCATATTTGATGGCAGACAGGGAGTCCGTCGCCACCGACAGGCCCGCGATGCCGCATGCCATAGTGCGa                         >  1:231455/1‑147 (MQ=255)
                                                                                                                           ttgCCGTACTGCGGATATTCACCGTCGATTTCAAAGTCCACCGCCAGGCCGTTTTCGTCACGGATTGGTTTCACGCGGGCATATTTGATGGCAGACAGGGAGTCCGTCGCCACCGACAGGCCCGCGATGCCGCATGCCATAGTGCGa                         <  2:231455/147‑1 (MQ=255)
                                                                                                                           ttgCCGTACTGCGGATATTCACCGTCGATTTCAAAGTCCACCGCCAGGCCGTTTTCGTCACGGATTGGTTTCACGCGGGCATATTTGATGGCAGACAGGGAGTCCGTCGCCACCGACAGGCCCGCGATGCCGCATGCCATAGTGCGATa                       <  1:181879/149‑1 (MQ=255)
                                                                                                                             gCCGTACTGCGGATATTCACCGTCGATTTCAAAGTCCACCGCCAGGCCGTTTTCGTCACGGATTGGTTTCACGCGGGCATATTTGATGGCAGACAGGGAGTCCGTCGCCACCGACAGGCCCGCGATGCCGCATGCCATAGTGCGATAg                      <  5:57611/148‑1 (MQ=255)
                                                                                                                             gCCGTACTGCGGATATTCACCGTCGATTTCAAAGTCCACCGCCAGGCCGTTTTCGTCACGGATTGGTTTCACGCGGGCATATTTGATGGCAGACAGGGAGTCCGTCGCCACCGACAGGCCCGCGATGCCGCATGCCATAGTGCGATAGa                     >  3:35744/1‑149 (MQ=255)
                                                                                                                                      cGGATATTCACCGTCGATTTCAAAGTCCACCGCCAGGCCGTTTTCGTCACGGATTGGTTTCACGCGGGCATATTTGATGGCAGACAGGGAGTCCGTCGCCACCGACAGGCCCGCGATGCCGCATGCCATAGTGCGATAGACATCACGAt            >  2:106323/1‑149 (MQ=255)
                                                                                                                                      cGGATATTCACCGGCGATTTCAAAGTCCACCGCCAGGCCGTTTTCGTCACGGATTGGTTTCACGCGGGCATATTTGATGGCAGACAGGGAGTCCGTCGCCACCGACAGGCCCGCGATGCCGCATGCCATAGTGCGATAGACATCACGAt            <  6:70172/149‑1 (MQ=255)
                                                                                                                                       ggATATTCACCGTCGATTTCAAAGTCCACCGCCAGGCCGTTTTCGTCACGGATTGGTTTCACGCGGGCATATTTGATGGCAGACAGGGAGTCCGTCGCCACCGACAGGCCCGCGATGCCGCATGCCATAGTGCGATAGACATCACGATc           <  2:97280/149‑1 (MQ=255)
                                                                                                                                            ttCACCGTCGATTTCAAAGTCCACCGCCAGGCCGTTTTCGTCACGGATTGGTTTCACGCGGGCATATTTGATGGCAGACAGGGAGTCCGTCGCCACCGACAGGCCCGCGATGCCGCATGCCATAGTGCGATAGACATCACGATCGTg        >  6:201074/1‑147 (MQ=255)
                                                                                                                                            ttCACCGTCGATTTCAAAGTCCACCGCCAGGCCGTTTTCGTCACGGATTGGTTTCACGCGGGCATATTTGATGGCAGACAGGGAGTCCGTCGCCACCGACAGGCCCGCGATGCCGCATGCCATAGTGCGATAGACATCACGATCGTGCa      <  1:83643/149‑1 (MQ=255)
                                                                                                                                               aCCGTCGATTTCAAAGTCCACCGCCAGGCCGTTTTCGTCACGGATTGGTTTCACGCGGGCATATTTGATGGCAGACAGGGAGTCCGTCGCCACCGACAGGCCCGCGATGCCGCATGCCATAGTGCGATAGACATCACGATCGTGCAgcg   >  2:99199/1‑149 (MQ=255)
                                                                                                                                                ccGTCGATTTCAAAGTCCACCGCCAGGCCGTTTTCGTCACGGATTGGTTTCACGCGGGCATATTTGATGGCAGACAGGGAGTCCGTCGCCACCGACAGGCCCGCGATGCCGCATGCCATAGTGCGATAGACATCACGATCGTGCAgcgc  >  2:23074/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TTAGAAGTGATAGTCAGAATCGACTGGGTAGGGACGGCGTTGCGATAGGTTGGCAGCGCTTTAATTTTCTTCATAAAGCGTTCAACCAGGTCGCAGGCAATGCTGTCTACGCGCTCGTCGTTGTTGCCGTACTGCGGATATTCACCGTCGATTTCAAAGTCCACCGCCAGGCCGTTTTCGTCACGGATTGGTTTCACGCGGGCATATTTGATGGCAGACAGGGAGTCCGTCGCCACCGACAGGCCCGCGATGCCGCATGCCATAGTGCGATAGACATCACGATCGTGCAGCGC  >  NC_000913/3260578‑3260870

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