Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,146,436 Δ1 bp 7.5% coding (179/879 nt) pflC → putative [formate‑C‑acetyltransferase 2]‑activating enzyme; pyruvate formate‑lyase 1‑activating enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,146,4360T. 7.5% 290.8 / 16.5 80coding (179/879 nt)pflCputative [formate‑C‑acetyltransferase 2]‑activating enzyme; pyruvate formate‑lyase 1‑activating enzyme
Reads supporting (aligned to +/- strand):  ref base T (38/36);  new base . (3/3);  total (41/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.72e-01

GTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGTCTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGT  >  NC_000913/4146289‑4146577
                                                                                                                                                    |                                                                                                                                             
gTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGTCTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGAt                                                                                                                                               <  1:565455/148‑1 (MQ=255)
            aaaCGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGTCTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAAt                                                                                                                                   >  4:350613/1‑148 (MQ=255)
                        atgatgTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGTCTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACg                                                                                                                       >  4:652861/1‑148 (MQ=255)
                         tgatgTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGTCTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGATCGACTCCgga                                                                                                                                        >  1:111911/1‑129 (MQ=255)
                         tgatgTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGTCTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGATCGACTCCgga                                                                                                                                        <  2:111911/130‑2 (MQ=255)
                           atgTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGTCTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCag                                                                                                                    <  3:652861/148‑1 (MQ=255)
                                   cgcATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGTCTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGATCGACTCCggatt                                                                                                                                      >  4:958390/1‑119 (MQ=255)
                                   cgcATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGTCTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGATCGACTCCggatt                                                                                                                                      <  3:958390/122‑4 (MQ=255)
                                     cATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGTCTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGaaa                                                                                                          >  2:473380/1‑148 (MQ=255)
                                         ttCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGTCTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTc                                                                                                      <  3:350613/148‑1 (MQ=255)
                                             aCATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGTCTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAg                                                                                                                               >  1:607086/1‑119 (MQ=255)
                                             aCATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGTCTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAg                                                                                                                               <  2:607086/119‑1 (MQ=255)
                                              cATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGTCTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGATCGACTCCggatt                                                                                                                                      >  4:901748/1‑108 (MQ=255)
                                              cATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGTCTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGATCGACTCCggatt                                                                                                                                      <  3:901748/111‑4 (MQ=255)
                                              cATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGTCTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAg                                                                                                                               >  1:608385/1‑118 (MQ=255)
                                              cATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGTCTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAg                                                                                                                               <  2:608385/118‑1 (MQ=255)
                                                ttCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGTCTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGTTACGCagagag                                                                                                                >  2:264600/1‑131 (MQ=255)
                                                ttCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGTCTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCagagag                                                                                                                <  1:264600/131‑1 (MQ=255)
                                                ttCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGTCTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCag                                                                                                                    >  4:810288/1‑127 (MQ=255)
                                                ttCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGTCTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCag                                                                                                                    <  3:810288/127‑1 (MQ=255)
                                                 tCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGTCTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACtgt                                                                                              >  1:256959/1‑148 (MQ=255)
                                                        tATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGTCTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAgg                                                                                                               <  2:834217/124‑1 (MQ=255)
                                                        tATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGTCTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAgg                                                                                                               >  1:834217/1‑124 (MQ=255)
                                                        tATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGTCTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAAt                                                                                       >  3:923552/1‑148 (MQ=255)
                                                              cTGAATGACGGCGAGGGCATTCGTACGGTGGTCTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGaaa                                                                                                          >  4:881831/1‑123 (MQ=255)
                                                              cTGAATGACGGCGAGGGCATTCGTACGGTGGTCTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGaaa                                                                                                          <  3:881831/123‑1 (MQ=255)
                                                              cTGAATGACGGCGAGGGCATTCGTACGGTGGTCTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATg                                                                                                        <  1:709887/125‑1 (MQ=255)
                                                              cTGAATGACGGCGAGGGCATTCGTACGGTGGTCTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATg                                                                                                        >  2:709887/1‑125 (MQ=255)
                                                                gAATGACGGCGAGGGCATTCGTACGGTGGTCTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCg                                                                                <  4:47329/147‑1 (MQ=255)
                                                                gAATGACGGCGAGGGCATTCGTACGGTGGTCTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCg                                                                                >  3:47329/1‑147 (MQ=255)
                                                                    gACGGCGAGGGCATTCGTACGGTGGTCTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATg                                                                           >  2:318141/1‑148 (MQ=255)
                                                                            gggCATTCGTACGGTGGTCTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGaa                                                                   <  1:473380/148‑1 (MQ=255)
                                                                              gCATTCGTACGGTGGTCTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATg                                                                 <  2:256959/148‑1 (MQ=255)
                                                                                           ggtCTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGaaa                                                                                        >  1:444081/1‑112 (MQ=255)
                                                                                           ggtCTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGaaa                                                                                        <  2:444081/112‑1 (MQ=255)
                                                                                           ggtCTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCgg                                                                        <  4:872050/128‑1 (MQ=255)
                                                                                           ggtCTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCgg                                                                        >  3:872050/1‑128 (MQ=255)
                                                                                           ggtCTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCg                                                    <  2:1028173/148‑1 (MQ=255)
                                                                                              cTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCag                                                                                                                    >  4:714578/1‑81 (MQ=255)
                                                                                              cTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCag                                                                                                                    <  3:714578/81‑1 (MQ=255)
                                                                                                tttttAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGttt                                                                                   <  2:373417/112‑1 (MQ=255)
                                                                                                tttttAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGttt                                                                                   >  1:373417/1‑112 (MQ=255)
                                                                                                tttttAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGa                                               >  3:957233/1‑148 (MQ=255)
                                                                                                 ttttAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGaa                                              <  4:957233/148‑1 (MQ=255)
                                                                                                          cTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGaaa                                                                                                          >  1:1090722/1‑79 (MQ=255)
                                                                                                          cTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGaaa                                                                                                          <  2:1090722/79‑1 (MQ=255)
                                                                                                          cTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGAtt                                        >  1:452475/1‑145 (MQ=255)
                                                                                                          cTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGAtt                                        <  2:452475/145‑1 (MQ=255)
                                                                                                               cGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGtt                                                                                    <  4:963/96‑1 (MQ=255)
                                                                                                               cGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGtt                                                                                    >  3:963/1‑96 (MQ=255)
                                                                                                                gCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATg                                                                 <  4:592384/114‑1 (MQ=255)
                                                                                                                gCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATg                                                                 >  3:592384/1‑114 (MQ=255)
                                                                                                                 cATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACt                                                                                                >  2:885731/1‑82 (MQ=255)
                                                                                                                 cATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACt                                                                                                <  1:885731/82‑1 (MQ=255)
                                                                                                                  aTCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGatat                             <  4:431565/148‑1 (MQ=255)
                                                                                                                  aTCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGatat                             >  3:431565/1‑148 (MQ=255)
                                                                                                                   tCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATc                            <  1:318141/148‑1 (MQ=255)
                                                                                                                     tttGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGa                                                                       >  3:26021/1‑103 (MQ=255)
                                                                                                                     tttGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGa                                                                       <  4:26021/103‑1 (MQ=255)
                                                                                                                         cccGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATc                                                                                                                                  >  2:857598/1‑40 (MQ=255)
                                                                                                                         cccGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATc                                                                                                                                  <  1:857598/40‑1 (MQ=255)
                                                                                                                                     tAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCgg                                                        <  3:1100112/102‑1 (MQ=255)
                                                                                                                                     tAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCgg                                                        >  4:1100112/1‑102 (MQ=255)
                                                                                                                                     tAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCgg                                                        <  1:516672/102‑1 (MQ=255)
                                                                                                                                     tAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCgg                                                        >  2:516672/1‑102 (MQ=255)
                                                                                                                                     tAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGAt                                         >  3:132354/1‑117 (MQ=255)
                                                                                                                                     tAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGAt                                         <  4:132354/117‑1 (MQ=255)
                                                                                                                                      aaTCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGtt                                                                                    >  2:880235/1‑73 (MQ=255)
                                                                                                                                      aaTCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGtt                                                                                    <  1:880235/73‑1 (MQ=255)
                                                                                                                                      aaTCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATg                                                                 >  1:605451/1‑92 (MQ=255)
                                                                                                                                      aaTCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATg                                                                 <  2:605451/92‑1 (MQ=255)
                                                                                                                                        tCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGa                                                                       <  2:809273/84‑1 (MQ=255)
                                                                                                                                        tCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGa                                                                       >  1:809273/1‑84 (MQ=255)
                                                                                                                                         ccGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACggg      >  3:1125595/1‑148 (MQ=255)
                                                                                                                                         ccGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACggg      >  1:290456/1‑148 (MQ=255)
                                                                                                                                                tCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCgg                                                                        >  1:432727/1‑76 (MQ=255)
                                                                                                                                                tCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCgg                                                                        <  2:432727/76‑1 (MQ=255)
                                                                                                                                                   aTCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGc                                                                <  1:894876/81‑1 (MQ=255)
                                                                                                                                                   aTCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGc                                                                >  2:894876/1‑81 (MQ=255)
                                                                                                                                                   aTCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGt  <  4:532896/143‑1 (MQ=255)
                                                                                                                                                   aTCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGt  >  3:532896/1‑143 (MQ=255)
                                                                                                                                                    |                                                                                                                                             
GTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGTCTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGT  >  NC_000913/4146289‑4146577

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