Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,416,941 A→G 5.4% intergenic (+50/‑140) pta → / → yfcC phosphate acetyltransferase/putative inner membrane transporter; C4‑dicarboxylate anaerobic carrier family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,416,9410AG 5.4% 225.9 / 2.5 74intergenic (+50/‑140)pta/yfcCphosphate acetyltransferase/putative inner membrane transporter; C4‑dicarboxylate anaerobic carrier family protein
Reads supporting (aligned to +/- strand):  ref base A (35/34);  new base G (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 = 9.09e-01

TGCGCAAGCCGGTTAACGACCTGTCCCGTGGCGCACTGGTTGACGATATCGTCTACACCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCCCACTTTGCATACTTACCACTTTGTTTTGTGCAAGGGAATATTTGCGCTATGTCCG  >  NC_000913/2416794‑2417087
                                                                                                                                                   |                                                                                                                                                  
tGCGCAAGCCGGTTAACGACCTGTCCCGTGGCGCACTGGTTGACGATATCGTCTACACCATCGCGCTGACTGCGATTCAGTCTGGACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAaataat                                                                                                                                                   <  7:187514/149‑1 (MQ=255)
  cgcAAGCCGGTTAACGACCTGTCCCGTGGCGCACTGGTTGACGATATCGTCTACACCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCa                                                                                                                                                 <  8:121727/149‑1 (MQ=255)
   gcAAGCCGGTTAACGACCTGTCCCGTGGCGCACTGGTTGACGATATCGTCTACACCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCAc                                                                                                                                                >  8:185774/1‑149 (MQ=255)
     aaGCCGGTTAACGACCTGTCCCGTGGCGCACTGGTTGACGATATCGTCTACACCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTa                                                                                                                                              >  6:287749/1‑149 (MQ=255)
         cGGTTAACGACCTGTCCCGTGGCGCACTGGTTGACGATATCGTCTACACCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTc                                                                                                                                          <  6:233989/149‑1 (MQ=255)
              aaCGACCTGTCCCGTGGCGCACTGGTTGACGATATCGTCTACACCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGtt                                                                                                                                     >  8:54361/1‑149 (MQ=255)
                 gACCTGTCCCGTGGCGCACTGGTTGACGATATCGTCTACACCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGttttt                                                                                                                                  >  5:14860/1‑149 (MQ=255)
                 gACCTGTCCCGTGGCGCACTGGTTGACGATATCGTCTACACCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACAATTTCCGGttttt                                                                                                                                  >  8:65938/1‑149 (MQ=255)
                   ccTGTCCCGTGGCGCACTGGTTGACGATATCGTCTACACCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTa                                                                                                                                >  8:259631/1‑149 (MQ=255)
                     tGTCCCGTGGCGCACTGGTTGACGATATCGTCTACACCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTAtt                                                                                                                              >  1:40850/1‑149 (MQ=255)
                          cGTGGCGCACTGGTTGACGATATCGTCTACACCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTctct                                                                                                                          >  4:3149/1‑148 (MQ=255)
                           gTGGCGCACTGGTTGACGATATCGTCTACACCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTa                                                                                                                        <  1:281838/149‑1 (MQ=255)
                            tGGCGCACTGGTTGACGATATCGTCTACACCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTaa                                                                                                                       <  1:289737/149‑1 (MQ=255)
                                      gTTGACGATATCGTCTACACCATCGCGCTGACTGCTATTCAGTCTACACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTaa                                                                                                             <  1:249115/149‑1 (MQ=255)
                                       ttGACGATATCGTCTACACCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAAt                                                                                                            <  6:14860/149‑1 (MQ=255)
                                         gACGATATCGTCTACACCATCGCGCGGACTGCGATGCAGTCTGCACAGCAGCAGTACTCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATc                                                                                                           <  2:280246/148‑1 (MQ=255)
                                          aCGATATCGTCTACACCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCt                                                                                                         >  4:141557/1‑149 (MQ=255)
                                             atatCGTCTACACCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTc                                                                                                      <  3:136128/149‑1 (MQ=255)
                                                tCGTCTACACCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTg                                                                                                    >  7:71293/1‑148 (MQ=255)
                                                tCGTCTACACCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTg                                                                                                    <  8:71293/148‑1 (MQ=255)
                                                tCGTCTACACCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGa                                                                                                   <  1:318728/149‑1 (MQ=255)
                                                          cATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCtt                                                                                            >  4:275402/1‑146 (MQ=255)
                                                          cATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCtt                                                                                            <  3:275402/146‑1 (MQ=255)
                                                          cATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCt                                                                                         >  2:162564/1‑149 (MQ=255)
                                                          cATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCt                                                                                         >  8:180592/1‑149 (MQ=255)
                                                           aTCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTc                                                                                                      >  5:303122/1‑135 (MQ=255)
                                                           aTCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTc                                                                                                      <  6:303122/135‑1 (MQ=255)
                                                            tCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGGATTAATCCTTTCTGATTATCTTGCtt                                                                                        >  6:279380/1‑148 (MQ=255)
                                                              gcgcTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTaa                                                                                      <  2:8715/148‑1 (MQ=255)
                                                               cgcTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACt                                                                                    <  7:54361/149‑1 (MQ=255)
                                                               cgcTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACt                                                                                    <  6:303279/149‑1 (MQ=255)
                                                                   gACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTgcgc                                                                                >  3:206096/1‑149 (MQ=255)
                                                                    aCTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCt                                                                               <  3:3149/149‑1 (MQ=255)
                                                                      tGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTg                                                                              >  5:97133/1‑148 (MQ=255)
                                                                       gCGATTCAGTCTACACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATc                                                                                              <  8:3259/131‑1 (MQ=255)
                                                                         gATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCAt                                                                           >  7:286085/1‑148 (MQ=255)
                                                                           ttCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCaa                                                                        >  3:258207/1‑149 (MQ=255)
                                                                             cAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAAt                                                                       <  6:178833/148‑1 (MQ=255)
                                                                               gTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAGTGATTATTTCCGGtt                                                                                                                                     <  8:323600/84‑1 (MQ=38)
                                                                               gTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAGTGATTATTTCCGGtt                                                                                                                                     >  7:323600/1‑84 (MQ=38)
                                                                                  tGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTg                                                                 >  2:80392/1‑149 (MQ=255)
                                                                                               tAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTAAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCCCACtt                                                    >  7:225167/1‑149 (MQ=255)
                                                                                                 aTCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCCCACttt                                                   <  2:40850/148‑1 (MQ=255)
                                                                                                  tctcGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCCCACTTTg                                                  >  1:186273/1‑148 (MQ=255)
                                                                                                  tctcGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCCCACTTTg                                                  >  2:94034/1‑148 (MQ=255)
                                                                                                    tcGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTTATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATAGCGCCATCACACTTTGCa                                                >  5:168673/1‑148 (MQ=255)
                                                                                                      gTCTTCATCCGCAGGTTTGCGCTGCGTATATCTGTACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCCCACTTTGCAtt                                              <  7:65938/148‑2 (MQ=255)
                                                                                                      gTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCCCACTTTGCa                                                <  8:275463/146‑1 (MQ=255)
                                                                                                      gTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCCCACTTTGCa                                                <  1:80392/146‑1 (MQ=255)
                                                                                                      gTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCCCACTTTGCATAc                                             <  7:259631/149‑1 (MQ=255)
                                                                                                      gTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCCCACTTTGCATAc                                             <  5:251013/149‑1 (MQ=255)
                                                                                                      gTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCCCACTTTGCATAc                                             >  6:251013/1‑149 (MQ=255)
                                                                                                        catcatCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCCCACTTTGCATACtt                                           <  2:303224/149‑1 (MQ=255)
                                                                                                          tcatcCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCCCACTTTGCATACTTAc                                         >  2:304218/1‑149 (MQ=255)
                                                                                                            atcCGCAGCTTTGCGCTGCTGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATATGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATATCGCCATCCCACTTTGCATACTTAc                                         <  5:249411/147‑1 (MQ=255)
                                                                                                               cGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCCCACTTTGCATACTTACCACtt                                    >  5:275208/1‑149 (MQ=255)
                                                                                                                gCAGCTTTGCGCTGCGGATATCTGAACCGGAAAGAATCACTATTTCCGCTACTATATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCCCACTTTGCATACTTACCACttt                                   >  2:209837/1‑149 (MQ=255)
                                                                                                                    cTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCCCACTTTGCATACTTACCACtttgttt                               >  2:69239/1‑149 (MQ=255)
                                                                                                                         cgcTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCCCACTTTGCATACTTACCACTTTGTTTtgtg                           <  5:270793/148‑1 (MQ=255)
                                                                                                                         cgcTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCCCACTTTGCATACTTACCACTTTGTTTtgtg                           <  8:286085/148‑1 (MQ=255)
                                                                                                                         cgcTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCCCACTTTGCATACTTACCACTTTGTTTTGTGc                          <  1:304218/149‑1 (MQ=255)
                                                                                                                           cTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCCCACTTTGCATACTTACCACTTTGTTTTGTGCaa                        <  7:287565/149‑1 (MQ=255)
                                                                                                                              cGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCttt                                                                                                       <  4:298214/67‑1 (MQ=255)
                                                                                                                              cGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCttt                                                                                                       >  3:298214/1‑67 (MQ=255)
                                                                                                                              cGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCCCACTTTGCATACTTACCACTTTGTTTTGTGCAAgg                      <  2:186273/148‑1 (MQ=255)
                                                                                                                                gATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCCCACTTTGCATACTTACCACTTTGTTTTGTGCAAGGGaa                   >  8:50662/1‑149 (MQ=255)
                                                                                                                                 atatCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCCCACTTTGCATACTTACCACTTTGTTTTGTGCAAGGGAat                  <  4:206096/149‑1 (MQ=255)
                                                                                                                                      tGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCCCACTTTGCATACTTACCACTTTGTTTTGTGCAAGGGAATATTTg             >  8:132061/1‑149 (MQ=255)
                                                                                                                                      tGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCCCACTTTGCATACTTACCACTTTGTTTTGTGCAAGGGAATATTTg             >  8:1543/1‑149 (MQ=255)
                                                                                                                                       gAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGAGTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCCCACTTTGCATACTTACCCCTTTGTTTTGTGCAAGGGAATATTTgc            >  5:290963/1‑149 (MQ=255)
                                                                                                                                          ccGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCCCACTTTGCATACTTACCACTTTGTTTTGTGCAAGGGAATATTTgcgc          >  3:58733/1‑148 (MQ=255)
                                                                                                                                           cGGAAATAGTGATTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCCCACTTTGCATAc                                             <  2:278537/112‑1 (MQ=255)
                                                                                                                                           cGGAAATAGTGATTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCCCACTTTGCATAc                                             >  1:278537/1‑112 (MQ=255)
                                                                                                                                           cGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCCCACTTTGCATACTTACCACTTTGTTTTGTGCAAGGGAATATTTGCGCTa        <  6:228349/149‑1 (MQ=255)
                                                                                                                                            ggAAATACTCACTATTTCCTGTATTTTATTCTCTTAATATGCATTATTGCTTTGTGATTATCTTGCTTAACTGCGCTGCATGAATGTATTGCGCCATCCCACTTTGCATACTTACCACTTTGTTTTGTGGAAGGGAATATTTGCGCTAt       <  7:326841/149‑1 (MQ=255)
                                                                                                                                                 taatCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCCCACTTTGCATACTTACCACTTTGTTTTGTGCAAGGTAATATTTGCGCTATGTCCg  <  1:69239/149‑1 (MQ=255)
                                                                                                                                                 taatCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGAATCAATGAATTGCGCCATCCCACTTTGCATACTTACCACTTTGTTTTGTGCAAGGGAATATTTGCTCTATGTCCg  >  1:316457/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TGCGCAAGCCGGTTAACGACCTGTCCCGTGGCGCACTGGTTGACGATATCGTCTACACCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCCCACTTTGCATACTTACCACTTTGTTTTGTGCAAGGGAATATTTGCGCTATGTCCG  >  NC_000913/2416794‑2417087

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