Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,146,431:1 +A 7.8% coding (174/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,4311.A 7.8% 320.8 / 24.0 77coding (174/879 nt)pflCputative [formate‑C‑acetyltransferase 2]‑activating enzyme; pyruvate formate‑lyase 1‑activating enzyme
Reads supporting (aligned to +/- strand):  ref base . (36/35);  new base A (3/3);  total (39/38)
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.50e-01

GTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGTCTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGG  >  NC_000913/4146289‑4146573
                                                                                                                                               |                                                                                                                                              
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                                                                                                                           <  2:607086/119‑1 (MQ=255)
                                             aCATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGTCTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAg                                                                                                                           >  1:607086/1‑119 (MQ=255)
                                              cATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGTCTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGATCGACTCCggatt                                                                                                                                  <  3:901748/111‑4 (MQ=255)
                                              cATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGTCTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGATCGACTCCggatt                                                                                                                                  >  4:901748/1‑108 (MQ=255)
                                              cATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGTCTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAg                                                                                                                           <  2:608385/118‑1 (MQ=255)
                                              cATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGTCTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAg                                                                                                                           >  1:608385/1‑118 (MQ=255)
                                                ttCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGTCTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGTTACGCagagag                                                                                                            >  2:264600/1‑131 (MQ=255)
                                                ttCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGTCTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCagagag                                                                                                            <  1:264600/131‑1 (MQ=255)
                                                ttCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGTCTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCag                                                                                                                <  3:810288/127‑1 (MQ=255)
                                                ttCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGTCTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCag                                                                                                                >  4:810288/1‑127 (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                                                                            >  3:47329/1‑147 (MQ=255)
                                                                gAATGACGGCGAGGGCATTCGTACGGTGGTCTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCg                                                                            <  4:47329/147‑1 (MQ=255)
                                                                    gACGGCGAGGGCATTCGTACGGTGGTCTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATg                                                                       >  2:318141/1‑148 (MQ=255)
                                                                            gggCATTCGTACGGTGGTCTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGaa                                                               <  1:473380/148‑1 (MQ=255)
                                                                              gCATTCGTACGGTGGTCTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATg                                                             <  2:256959/148‑1 (MQ=255)
                                                                                           ggtCTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGaaa                                                                                    <  2:444081/112‑1 (MQ=255)
                                                                                           ggtCTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGaaa                                                                                    >  1:444081/1‑112 (MQ=255)
                                                                                           ggtCTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCgg                                                                    >  3:872050/1‑128 (MQ=255)
                                                                                           ggtCTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCgg                                                                    <  4:872050/128‑1 (MQ=255)
                                                                                           ggtCTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCg                                                <  2:1028173/148‑1 (MQ=255)
                                                                                              cTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAg                                                                                                                                                 <  4:4246/49‑1 (MQ=255)
                                                                                              cTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAg                                                                                                                                                 >  3:4246/1‑49 (MQ=255)
                                                                                              cTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCag                                                                                                                >  4:714578/1‑81 (MQ=255)
                                                                                              cTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCag                                                                                                                <  3:714578/81‑1 (MQ=255)
                                                                                                tttttAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCg                                                                                                                                             >  4:59108/1‑50 (MQ=255)
                                                                                                tttttAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCg                                                                                                                                             <  3:59108/50‑1 (MQ=255)
                                                                                                tttttAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGttt                                                                               >  1:373417/1‑112 (MQ=255)
                                                                                                tttttAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGttt                                                                               <  2:373417/112‑1 (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                                                                                >  3:963/1‑96 (MQ=255)
                                                                                                               cGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGtt                                                                                <  4:963/96‑1 (MQ=255)
                                                                                                                gCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATg                                                             >  3:592384/1‑114 (MQ=255)
                                                                                                                gCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATg                                                             <  4:592384/114‑1 (MQ=255)
                                                                                                                 cATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACt                                                                                            <  1:885731/82‑1 (MQ=255)
                                                                                                                 cATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACt                                                                                            >  2:885731/1‑82 (MQ=255)
                                                                                                                  aTCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGatat                         >  3:431565/1‑148 (MQ=255)
                                                                                                                  aTCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGatat                         <  4:431565/148‑1 (MQ=255)
                                                                                                                   tCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATc                        <  1:318141/148‑1 (MQ=255)
                                                                                                                     tttGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGa                                                                   <  4:26021/103‑1 (MQ=255)
                                                                                                                     tttGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGa                                                                   >  3:26021/1‑103 (MQ=255)
                                                                                                                         cccGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATc                                                                                                                              >  2:857598/1‑40 (MQ=255)
                                                                                                                         cccGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATc                                                                                                                              <  1:857598/40‑1 (MQ=255)
                                                                                                                                     tAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCgg                                                    >  4:1100112/1‑102 (MQ=255)
                                                                                                                                     tAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCgg                                                    <  1:516672/102‑1 (MQ=255)
                                                                                                                                     tAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCgg                                                    >  2:516672/1‑102 (MQ=255)
                                                                                                                                     tAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCgg                                                    <  3:1100112/102‑1 (MQ=255)
                                                                                                                                     tAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGAt                                     >  3:132354/1‑117 (MQ=255)
                                                                                                                                     tAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGAt                                     <  4:132354/117‑1 (MQ=255)
                                                                                                                                      aaTCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGtt                                                                                <  1:880235/73‑1 (MQ=255)
                                                                                                                                      aaTCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGtt                                                                                >  2:880235/1‑73 (MQ=255)
                                                                                                                                      aaTCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATg                                                             <  2:605451/92‑1 (MQ=255)
                                                                                                                                      aaTCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATg                                                             >  1:605451/1‑92 (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)
                                                                                                                                               |                                                                                                                                              
GTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGTCTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGG  >  NC_000913/4146289‑4146573

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