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,4380GA98.8% 254.6 / NA 80*330* (TAG→TAAyiiDGNAT family putative N‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (35/44);  minor base C (0/1);  total (35/45)
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

GACGAAAAGCACGGGTGCAGATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAGCTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCTATAAATTTTTATTCTTCGTTGGCTCCATACGTGTCTTGCTACAT  >  NC_000913/4078302‑4078568
                                                                                                                                        |                                                                                                                                  
gACGAAAAGCACGGGTGCAGATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACt                                                                                                                                  >  1:171396/1‑139 (MQ=255)
    aaaaGCACGGGTGCAGATGCATGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATg                                                                                                                              <  1:822802/139‑1 (MQ=255)
       aGCACGGGTGCAGATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGgag                                                                                                                           <  2:334292/139‑1 (MQ=255)
                     tGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGtgctgt                                                                                                             <  1:221506/139‑2 (MQ=255)
                     tGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGtgctgc                                                                                                             <  1:1044024/139‑1 (MQ=255)
                     tGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGtgctgc                                                                                                             <  1:35142/139‑1 (MQ=255)
                        aGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCctc                                                                                                          >  1:52155/1‑139 (MQ=255)
                        aGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCctc                                                                                                          >  2:871986/1‑139 (MQ=255)
                         ggTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCctct                                                                                                         <  2:171396/139‑1 (MQ=255)
                         ggTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCctct                                                                                                         >  1:366302/1‑139 (MQ=255)
                                   tttGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAgg                                                                                               <  2:52155/139‑1 (MQ=255)
                                   tttGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAgg                                                                                               <  2:876177/139‑1 (MQ=255)
                                     tGGCGACGAGACGCGGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGTTGCCTCTCTCATTCAGGTc                                                                                             <  2:845578/139‑1 (MQ=255)
                                     tGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGTCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAgggg                                                                                                                      <  2:514238/114‑1 (MQ=255)
                                     tGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGTCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAgggg                                                                                                                      >  1:514238/1‑114 (MQ=255)
                                     tGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTc                                                                                             <  1:810342/139‑1 (MQ=255)
                                     tGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTc                                                                                             <  1:252494/139‑1 (MQ=255)
                                       gCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCAt                                                                                           <  1:265807/139‑1 (MQ=255)
                                       gCGACGAGACGCCGGGTGCAGTGGTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCAt                                                                                           <  2:617252/139‑1 (MQ=255)
                                                  ccGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCAt                                                                                           <  2:692232/128‑1 (MQ=255)
                                                  ccGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCAt                                                                                           >  1:692232/1‑128 (MQ=255)
                                                        gCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGc                                                                          <  1:1040383/139‑1 (MQ=255)
                                                             gtTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTAttt                                                                     <  1:777464/139‑1 (MQ=255)
                                                                tGAAGGCACGTATATCGTTCTGCCCGCGAAGCAATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCgg                                                                  >  2:639490/1‑139 (MQ=255)
                                                                 gAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCAtt                                                                                          <  2:174006/114‑1 (MQ=255)
                                                                 gAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCAtt                                                                                          >  1:174006/1‑114 (MQ=255)
                                                                  aaGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGtt                                                                <  2:84135/139‑1 (MQ=255)
                                                                      cACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCAtt                                                                                          >  1:199517/1‑109 (MQ=255)
                                                                      cACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCAtt                                                                                          <  2:199517/109‑1 (MQ=255)
                                                                      cACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTatat                                                            >  2:353051/1‑139 (MQ=255)
                                                                       aCGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCAtt                                                                                                   >  1:15107/1‑99 (MQ=255)
                                                                       aCGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCAtt                                                                                                   <  2:15107/99‑1 (MQ=255)
                                                                              tCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAg                                                    <  2:366302/139‑1 (MQ=255)
                                                                                 ttCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCAtt                                                                                          >  2:193892/1‑98 (MQ=255)
                                                                                 ttCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCAtt                                                                                          <  1:193892/98‑1 (MQ=255)
                                                                                   cTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTc                                                                                                  >  1:641929/1‑88 (MQ=255)
                                                                                   cTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTc                                                                                                  <  2:641929/88‑1 (MQ=255)
                                                                                    tGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTc                                                                                             <  2:256199/92‑1 (MQ=255)
                                                                                    tGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTc                                                                                             >  1:256199/1‑92 (MQ=255)
                                                                                    tGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCgg                                                                  >  1:939423/1‑119 (MQ=255)
                                                                                    tGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCgg                                                                  <  2:939423/119‑1 (MQ=255)
                                                                                     gCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAgg                                                                                               <  1:344597/89‑1 (MQ=255)
                                                                                     gCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAgg                                                                                               >  2:344597/1‑89 (MQ=255)
                                                                                      cccGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAgggg                                                                                                                      <  2:984092/65‑1 (MQ=255)
                                                                                      cccGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAgggg                                                                                                                      >  1:984092/1‑65 (MQ=255)
                                                                                      cccGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTatat                                                                                     <  2:389497/98‑1 (MQ=255)
                                                                                      cccGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTatat                                                                                     >  1:389497/1‑98 (MQ=255)
                                                                                      cccGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCtat                                            >  1:948265/1‑139 (MQ=255)
                                                                                      cccGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCtat                                            >  2:413842/1‑139 (MQ=255)
                                                                                         gcgAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCAtt                                                                                          >  1:698056/1‑90 (MQ=255)
                                                                                         gcgAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCAtt                                                                                          <  2:698056/90‑1 (MQ=255)
                                                                                          cgAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCTATAAAt                                        >  1:663795/1‑139 (MQ=255)
                                                                                          cgAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCTATAAAt                                        >  2:206006/1‑139 (MQ=255)
                                                                                                 aTTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCTATAAATTTTTAtt                                 <  1:482819/139‑1 (MQ=255)
                                                                                                        ccGTATGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCTATAAATTTTTATTCTTCGtt                          <  2:948265/139‑1 (MQ=255)
                                                                                                             tGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAAct                                                <  2:102271/112‑1 (MQ=255)
                                                                                                             tGAAGAGGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAAct                                                >  1:102271/1‑112 (MQ=255)
                                                                                                                agagGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCgg                                                                  >  2:766612/1‑91 (MQ=255)
                                                                                                                agagGGCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCgg                                                                  <  1:766612/91‑1 (MQ=255)
                                                                                                                    ggCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCgg                                                                  <  2:395737/87‑1 (MQ=255)
                                                                                                                    ggCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCgg                                                                  >  1:395737/1‑87 (MQ=255)
                                                                                                                     gCGGGAACGAAGAAGAGTAACTGATGGTGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCTATAAATTTTTATTCTTCGTTGGCt                      <  2:399984/130‑1 (MQ=255)
                                                                                                                     gCGGGAACGAAGAAGAGTAACTGATGGTGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCTATAAATTTTTATTCTTCGTTGGCt                      >  1:399984/1‑130 (MQ=255)
                                                                                                                     gCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCTATAAATTTTTATTCTTCGTTGGCTCCATACgtg             >  2:289524/1‑139 (MQ=255)
                                                                                                                     gCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCTATAAATTTTTATTCTTCGTTGGCTCCATACgtg             >  1:179502/1‑139 (MQ=255)
                                                                                                                     gCGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCTATAAATTTTTATTCTTCGTTGGCTCCATACgtg             >  1:745383/1‑139 (MQ=255)
                                                                                                                      cgtgAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAg                                                                                                <  1:724922/52‑1 (MQ=255)
                                                                                                                      cGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCAtt                                                                                                   <  1:414111/52‑1 (MQ=255)
                                                                                                                      cGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCAtt                                                                                                   >  2:414111/1‑52 (MQ=255)
                                                                                                                      cGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAg                                                                                                >  2:724922/1‑55 (MQ=255)
                                                                                                                      cGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCTATAAATTTTTATTCTTCGTTGGc                       <  2:1037573/128‑1 (MQ=255)
                                                                                                                      cGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCTATAAATTTTTATTCTTCGTTGGc                       >  1:1037573/1‑128 (MQ=255)
                                                                                                                      cGGGAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCTATAAATTTTTATTCTTCGTTGGCTCCATACgtgt            <  2:179502/139‑1 (MQ=255)
                                                                                                                        ggAACGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCTATAAATTTTTATTCTTCGTTGGCTCCATACGTGTCt          <  2:571327/139‑1 (MQ=255)
                                                                                                                            cGAAGAAGAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCTATAAATTTTTATTCTTCGTTGGCTCCATACGTGTCTTGCt      <  2:745383/139‑1 (MQ=255)
                                                                                                                              aagaagAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAg                                                                                                >  2:121284/1‑47 (MQ=255)
                                                                                                                              aagaagAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAg                                                                                                <  1:121284/47‑1 (MQ=255)
                                                                                                                              aagaagAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTAt                                                                       <  2:315931/72‑1 (MQ=255)
                                                                                                                              aagaagAGTAACTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTAt                                                                       >  1:315931/1‑72 (MQ=255)
                                                                                                                                gaagaGTACCTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTCATATTAAATAAGAACTCTATAAATTTTTATTCTTCGTTGGCTCCATACGTGTCTTGCTACAt  <  1:639490/139‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                  
GACGAAAAGCACGGGTGCAGATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAGCTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGGTCATTTATATTAGTCCTTGGCTATTTCGGTTATATTAAATAAGAACTCTATAAATTTTTATTCTTCGTTGGCTCCATACGTGTCTTGCTACAT  >  NC_000913/4078302‑4078568

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: