Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,078,438 G→A *330* (TAG→TAA yiiD → GNAT family putative N‑acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,078,4380GA97.7% 273.3 / ‑6.3 87*330* (TAG→TAAyiiDGNAT family putative N‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (37/48);  minor base T (1/1);  total (38/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCACGGGTGCAGATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAGCTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCTATAAATTTTTATTCTTCGTTGGCTCCATACGTGTCTTGCTACATCA  >  NC_000913/4078310‑4078570
                                                                                                                                |                                                                                                                                    
gCACGGGTGCAGATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGgaga                                                                                                                            >  1:84206/1‑139 (MQ=255)
   cGGGTGCAGATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAggg                                                                                                                         >  2:248735/1‑139 (MQ=255)
       tGCAGATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCgg                                                                                                                     >  1:730021/1‑139 (MQ=255)
            aTGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGtgctg                                                                                                                <  1:633655/139‑1 (MQ=255)
             tGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGtgctgc                                                                                                               <  2:1522646/139‑1 (MQ=255)
             tGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGtgctgc                                                                                                               <  2:84206/139‑1 (MQ=255)
             tGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGtgctgc                                                                                                               <  2:837427/139‑1 (MQ=255)
                 ggTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCctct                                                                                                           >  2:1191939/1‑139 (MQ=255)
                   tCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCctctct                                                                                                         >  2:914466/1‑139 (MQ=255)
                       aaTCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCAtt                                                                                                     <  2:840643/139‑1 (MQ=255)
                          cTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAg                                                                                                  <  2:607864/139‑1 (MQ=255)
                           tttGGCGACGAGACGCCTGGTGCAGTGTTTGAAGGCACGTTTATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGTAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCACTCTCATTCAgg                                                                                                 <  1:1769490/139‑1 (MQ=255)
                           tttGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAgg                                                                                                 <  1:578389/139‑1 (MQ=255)
                           tttGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAgg                                                                                                 <  2:52011/139‑1 (MQ=255)
                           tttGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAgg                                                                                                 <  1:1426636/139‑1 (MQ=255)
                           tttGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAgg                                                                                                 <  2:815078/139‑1 (MQ=255)
                           tttGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAgg                                                                                                 <  1:914466/139‑1 (MQ=255)
                               gCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCAt                                                                                             <  1:113151/139‑1 (MQ=255)
                                cgacgaGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCAtt                                                                                            <  1:502914/139‑1 (MQ=255)
                                          ccGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAg                                                                                                  <  1:897126/123‑1 (MQ=255)
                                          ccGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAg                                                                                                  >  2:897126/1‑123 (MQ=255)
                                                    tgtTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTAtt                                                                        <  1:673427/139‑1 (MQ=255)
                                                           aGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGt                                                                   <  1:493149/137‑1 (MQ=255)
                                                           aGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGt                                                                   >  2:493149/1‑137 (MQ=255)
                                                              cACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAgggg                                                                                                                        <  1:1710726/81‑1 (MQ=255)
                                                              cACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAgggg                                                                                                                        >  2:1710726/1‑81 (MQ=255)
                                                              cACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTatat                                                              <  1:170406/139‑1 (MQ=255)
                                                              cACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTatat                                                              >  2:170406/1‑139 (MQ=255)
                                                              cACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTatat                                                              <  1:1486691/139‑1 (MQ=255)
                                                                cGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAggg                                                                                                                         <  1:632718/78‑1 (MQ=255)
                                                                cGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAggg                                                                                                                         >  2:632718/1‑78 (MQ=255)
                                                                cGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTa                                                            >  2:697330/1‑139 (MQ=255)
                                                                        gTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGaa                                                    <  1:486603/139‑1 (MQ=255)
                                                                        gTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGaa                                                    >  2:486603/1‑139 (MQ=255)
                                                                         ttCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTGACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAgg                                                                                                 >  2:1077754/1‑93 (MQ=255)
                                                                         ttCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTGACTGATGGAGAGGGGCGATGCTGCCTCTCTCATTCAgg                                                                                                 <  1:1077754/93‑1 (MQ=255)
                                                                         ttCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAggg                                                                                                                         >  1:405141/1‑69 (MQ=255)
                                                                         ttCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAggg                                                                                                                         <  2:405141/69‑1 (MQ=255)
                                                                         ttCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGaa                                                    <  1:433523/138‑1 (MQ=255)
                                                                         ttCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGaa                                                    >  2:433523/1‑138 (MQ=255)
                                                                          tCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCAtt                                                                                                     <  1:1006326/88‑1 (MQ=255)
                                                                          tCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCAtt                                                                                                     >  2:1006326/1‑88 (MQ=255)
                                                                           cTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAgggg                                                                                                                        <  2:1513096/68‑1 (MQ=255)
                                                                           cTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAgggg                                                                                                                        >  1:1513096/1‑68 (MQ=255)
                                                                           cTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATAtt                                                                                      <  1:782903/102‑1 (MQ=255)
                                                                           cTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATAtt                                                                                      >  2:782903/1‑102 (MQ=255)
                                                                           cTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAActc                                                 <  1:1191939/139‑1 (MQ=255)
                                                                           cTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAActc                                                 >  2:1400338/1‑139 (MQ=255)
                                                                            tgtccGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCAt                                                                                             <  2:115089/91‑1 (MQ=255)
                                                                            tgtccGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCAt                                                                                             >  1:115089/4‑94 (MQ=255)
                                                                            tGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAgggg                                                                                                                        <  2:859870/67‑1 (MQ=255)
                                                                            tGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAgggg                                                                                                                        >  1:859870/1‑67 (MQ=255)
                                                                               ccGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCtata                                             <  1:1400338/139‑1 (MQ=255)
                                                                                  cgAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTATCTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGaa                                                    >  2:559839/1‑129 (MQ=255)
                                                                                  cgAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTATCTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGaa                                                    <  1:559839/129‑1 (MQ=255)
                                                                                   gAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAgggg                                                                                                                        <  1:1222144/60‑1 (MQ=255)
                                                                                   gAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCTATAAAtt                                         >  2:492538/1‑139 (MQ=255)
                                                                                   gAAGCCATTTGGCCCGTATAAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAgggg                                                                                                                        >  2:1222144/1‑60 (MQ=255)
                                                                                              gCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCTATAAATTTTTATTCTTCg                              <  2:730021/139‑1 (MQ=255)
                                                                                               cccGTATGACGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTtata                                                               >  2:1229952/1‑105 (MQ=255)
                                                                                               cccGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGGTATTTCGGTtata                                                               <  1:1229952/105‑1 (MQ=255)
                                                                                                            ggCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGtg                                                                                                                   <  1:33477/40‑1 (MQ=255)
                                                                                                            ggCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGtg                                                                                                                   >  2:33477/1‑40 (MQ=255)
                                                                                                            ggCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTAt                                                                         >  1:1380048/1‑82 (MQ=255)
                                                                                                            ggCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTAt                                                                         <  2:1380048/82‑1 (MQ=255)
                                                                                                             gCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAgg                                                                                                 >  1:1632169/1‑57 (MQ=255)
                                                                                                             gCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAgg                                                                                                 <  2:1632169/57‑1 (MQ=255)
                                                                                                             gCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGt                                                                   <  1:1373499/87‑1 (MQ=255)
                                                                                                             gCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGt                                                                   >  2:1373499/1‑87 (MQ=255)
                                                                                                             gCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCTATAAATTTTTATTCTTCGTTGGCTCCATACgtg               >  1:1215893/1‑139 (MQ=255)
                                                                                                              cGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCctctc                                                                                                          >  1:1312825/1‑47 (MQ=255)
                                                                                                              cGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCctctc                                                                                                          <  2:1312825/47‑1 (MQ=255)
                                                                                                              cGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAgg                                                                                                 >  2:1215890/1‑56 (MQ=255)
                                                                                                              cGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAgg                                                                                                 <  1:1215890/56‑1 (MQ=255)
                                                                                                              cGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCTATAAATTTTTATTCTTCGTTGGCTCCATACgtgt              >  2:1086639/1‑139 (MQ=255)
                                                                                                              cGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCTATAAATTTTTATTCTTCGTTGGCTCCATACgtgt              >  2:1360872/1‑139 (MQ=255)
                                                                                                                 gAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAgg                                                                                                 <  2:239573/53‑1 (MQ=255)
                                                                                                                 gAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAgg                                                                                                 >  1:239573/1‑53 (MQ=255)
                                                                                                                    cGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTc                                                                                  <  2:650310/65‑1 (MQ=255)
                                                                                                                    cGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTc                                                                                  >  1:650310/1‑65 (MQ=255)
                                                                                                                     gaagaaGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTAATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCTATAAATTTTTATTCTTCGTTGGCTCCATACGTGTCTTGCTa       <  1:697330/139‑1 (MQ=255)
                                                                                                                       agaagaGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAgg                                                                                                 >  2:342806/1‑47 (MQ=255)
                                                                                                                       agaagaGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAgg                                                                                                 <  1:342806/47‑1 (MQ=255)
                                                                                                                          agagTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCTATAAATTTTTATTCTTCGTTGGCTCCATACGTGTCTTGCTACATCa  <  2:1215893/139‑1 (MQ=255)
                                                                                                                          agagTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCTATAAATTTTTATTCTTCGTTGGCTCCATACGTGTCTTGCTACATCa  <  1:1360872/139‑1 (MQ=255)
                                                                                                                            agTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCAtt                                                                                            <  2:1445948/47‑1 (MQ=255)
                                                                                                                            agTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCAtt                                                                                            >  1:1445948/1‑47 (MQ=255)
                                                                                                                                |                                                                                                                                    
GCACGGGTGCAGATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAGCTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCTATAAATTTTTATTCTTCGTTGGCTCCATACGTGTCTTGCTACATCA  >  NC_000913/4078310‑4078570

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: