Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,566,875 A→G R29R (CGT→CGC glgC ← glucose‑1‑phosphate adenylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,566,8750AG99.7% 953.3 / ‑5.1 309R29R (CGT→CGCglgCglucose‑1‑phosphate adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (218/90);  total (219/90)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.39e-01

Only 100 of 398 total aligned reads displayed.

GACAGCGCAAAGTCGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCACGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTAACCATGAC  >  NC_000913/3566786‑3566964
                                                                                          |                                                                                         
    gcgcAAAGTCGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTAGCGCGtcctc                                                                                    >  1:5514945/1‑93 (MQ=255)
    gcgcAAAGTCGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGtcctc                                                                                    >  1:3382840/1‑93 (MQ=255)
    gcgcAAAGTCGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGtcctc                                                                                    >  2:7658957/1‑93 (MQ=255)
     cgcAAAGTCGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGtcctcc                                                                                   <  1:2460061/93‑1 (MQ=255)
       cAAAGTCGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCg                                                                                 >  1:4687389/1‑93 (MQ=255)
       cAAAGTCGATAATGCGGAACTTACCACCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCg                                                                                 >  1:7983647/1‑93 (MQ=255)
        aaaGTCGATAATGCGGAACTTCCCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGc                                                                                <  2:3759634/93‑1 (MQ=255)
        aaaGTCGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGc                                                                                >  1:5809168/1‑93 (MQ=255)
         aaGTCGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGcc                                                                               >  2:7755082/1‑93 (MQ=255)
         aaGTCGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGcc                                                                               >  2:1356409/1‑93 (MQ=255)
         aaGTCGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGcc                                                                               <  1:1040927/93‑1 (MQ=255)
          aGTCGATAATGGGGACTTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTGCCGCGTCCTCCCGCCa                                                                              <  2:6178380/93‑1 (MQ=255)
           gTCGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAg                                                                             >  2:1639821/1‑93 (MQ=255)
           gTCGATAAGGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAACCCTTCAGGCGGGTACCGCGTCCTCCCGCCAg                                                                             <  2:2274650/93‑1 (MQ=255)
            tCGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGt                                                                            >  2:4635775/1‑93 (MQ=255)
            tCGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGt                                                                            >  1:8401853/1‑93 (MQ=255)
             cGATATGGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTa                                                                           <  2:2589161/93‑1 (MQ=255)
              gATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTAt                                                                          >  2:5031843/1‑93 (MQ=255)
              gATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTAt                                                                          >  1:3899918/1‑93 (MQ=255)
                tAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCa                                                                        >  1:1002440/1‑93 (MQ=255)
                caaTGCGTACCTCACCGGCGAATTGGCCGGCCGTTTTTTCTCGCTTATTGGTTAATTCCTGCAGGCGGGTACCGCGTCCTCCCGCCAGTATCa                                                                        <  2:2897330/92‑1 (MQ=255)
                 aaTGCGGAACTTCCCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAg                                                                       <  2:5428679/93‑1 (MQ=255)
                  aTGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAgg                                                                      >  2:5243940/1‑93 (MQ=255)
                  aTGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAgg                                                                      >  2:870182/1‑93 (MQ=255)
                  aTGCGGAATTTCCGCCGGAAGTGTCCGGCGGGTTTTGCTCGTTTATTGGTTAAATCTTGAGGGCGGGTACCGCGTCCTCCCGCCAGTATCAgg                                                                      <  2:63143/93‑1 (MQ=255)
                      ggACCTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCaa                                                                  <  2:5340811/93‑1 (MQ=255)
                      ggAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCaa                                                                  <  2:3127180/93‑1 (MQ=255)
                       gAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAAc                                                                 >  2:419543/1‑93 (MQ=255)
                       gAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGAACCGCGTACTCCCGCCAGTATCAGGTCAAc                                                                 >  2:2398208/1‑93 (MQ=255)
                         aCTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAg                                                               <  1:6231528/93‑1 (MQ=255)
                            tACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGAtt                                                            >  2:3838063/1‑93 (MQ=255)
                            tACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGAtt                                                            >  2:6500906/1‑93 (MQ=255)
                             cccgccgAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGAttt                                                           <  1:2309085/92‑1 (MQ=255)
                              ccgccgAGGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTCCCGCGTCCCCCCGCCAGTATCAGGGCAACAGATTTc                                                          <  1:9902951/93‑1 (MQ=255)
                              ccgccgAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTc                                                          >  1:314570/1‑93 (MQ=255)
                               cgccgAGGTTTACGGCCGTTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGTTCCCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCa                                                         <  2:3963840/93‑1 (MQ=255)
                               cgccgAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCa                                                         <  1:6899119/93‑1 (MQ=255)
                                gccgAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCaa                                                        >  2:5650590/1‑93 (MQ=255)
                                gccgAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCaa                                                        >  2:5481258/1‑93 (MQ=255)
                                gccgAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCaa                                                        >  1:2862305/1‑93 (MQ=255)
                                 ccgAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAAt                                                       >  2:2758233/1‑93 (MQ=255)
                                 ccgAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAAt                                                       <  2:3443099/93‑1 (MQ=255)
                                    gAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATgg                                                     <  1:6796009/93‑1 (MQ=255)
                                    gAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCCTCCTCCCCCCAGTATCAGGGCAACAAATTTCATTgg                                                     >  2:8433313/1‑93 (MQ=255)
                                     aaGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGc                                                    >  1:4992778/1‑93 (MQ=255)
                                     aaGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGc                                                    >  1:2515710/1‑93 (MQ=255)
                                     aaGTGTACGGCCGGTTTTGCCCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGc                                                    <  1:5562038/93‑1 (MQ=255)
                                      aGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTCCCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCa                                                   <  2:8525/93‑1 (MQ=255)
                                      aGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCa                                                   <  2:1426194/93‑1 (MQ=255)
                                       gtgtACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAg                                                  >  1:9602021/1‑93 (MQ=255)
                                        tgtACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGc                                                 >  1:6166742/1‑93 (MQ=255)
                                        tgtACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGc                                                 >  1:453064/1‑93 (MQ=255)
                                         gtACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCt                                                <  2:6749899/93‑1 (MQ=255)
                                          tACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTg                                               >  1:6379609/1‑93 (MQ=255)
                                           aCGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTgg                                              >  1:9556518/1‑93 (MQ=255)
                                              gCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGgtgc                                           >  2:4462585/1‑93 (MQ=255)
                                                 ggTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGcc                                        >  2:6830820/1‑93 (MQ=255)
                                                 ggTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGcc                                        >  1:1515913/1‑93 (MQ=255)
                                                 ggTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGGACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGcc                                        >  2:7352792/1‑93 (MQ=255)
                                                   ttttGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCaa                                      >  1:7096781/1‑93 (MQ=255)
                                                    tttGCTCGCTTATTGGTTAAATCCTTCAGGGGGGTACCGCGTACTCCCGCCAGTATCAGGGCAACAGCTTTCAATGGCAGCTTGTGCGCCaaa                                     >  2:9947065/1‑92 (MQ=255)
                                                     ttGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGACATCAATGGCAGCTGGCGCGCCAACa                                    >  2:3631723/1‑93 (MQ=255)
                                                     tcgCTCGCTTATTGGTAAAACCTTGCAGGCGGGTCCCGCTTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACa                                    <  2:978016/91‑1 (MQ=255)
                                                      tGCGCGCGTAGTGTTTAACTCTTTGAGGCGTCTACCGCCTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACAt                                   <  1:1753798/93‑1 (MQ=255)
                                                         tCGCTTATTGTTTAAATCCTTCAGGCGGGTCCCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTaa                                <  2:1221648/93‑1 (MQ=255)
                                                         tCGCTTATTGGTTAAACCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTaa                                <  2:9305313/93‑1 (MQ=255)
                                                         ccGCTTATGGGTTAAACCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTaa                                <  1:2739478/92‑1 (MQ=255)
                                                          cGCTTATTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAg                               >  1:3755546/1‑93 (MQ=255)
                                                            cTTATTGTTAAAATCCTTCAGGCGGCTCCCGCGCCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTg                             <  2:7191121/93‑1 (MQ=255)
                                                             tttttGGTAAAATCCTTCAGGCGGGTCCCGCGTCCCCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGa                            <  2:8839878/90‑1 (MQ=255)
                                                               ataggTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGGTGGCGCGCCAACATTAAGTGATc                          >  1:749543/4‑93 (MQ=255)
                                                               aTTGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATc                          >  2:2461700/1‑93 (MQ=255)
                                                                 tGGTTAAATCCTTCAGGCGGGTCCCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGt                        <  1:5251846/93‑1 (MQ=255)
                                                                 tGGTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGt                        >  1:833560/1‑93 (MQ=255)
                                                                  ggTTAAATCCTTCGGGCGGGTACCGCTCCCCCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGtt                       <  1:7385473/93‑1 (MQ=255)
                                                                  ggTTAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGtt                       >  2:1971537/1‑93 (MQ=255)
                                                                    ttAAATCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGttct                     <  1:1760191/93‑1 (MQ=255)
                                                                       aaTCCTTCAGGGGGGTCCCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTtcttct                  <  2:5130405/93‑1 (MQ=255)
                                                                       aaTCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGttcttcc                  >  1:9761930/1‑92 (MQ=255)
                                                                       aaTCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTtcttct                  >  2:8103297/1‑93 (MQ=255)
                                                                        aTCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTtctc                 >  2:5942437/1‑93 (MQ=255)
                                                                          ccTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTa               >  1:8791719/1‑93 (MQ=255)
                                                                           cTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTaa              >  2:1845693/1‑93 (MQ=255)
                                                                            ttCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTaaa             >  1:7806893/1‑93 (MQ=255)
                                                                            tgcgggCGGGTACCGCGTCCCCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTaaa             <  2:1396841/89‑1 (MQ=255)
                                                                             tCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAAc            >  2:7990305/1‑93 (MQ=255)
                                                                               aGGGGGTTACCGCGCCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTa          <  1:4331878/93‑1 (MQ=255)
                                                                               aGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCGACATTAAGTGATCGTTCTTCTCTAAACTa          >  1:8233709/1‑93 (MQ=255)
                                                                               aGGCGGGCCCCGCGCCCCCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTa          <  1:1910321/93‑1 (MQ=255)
                                                                                ggCGGGTCCCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTaa         <  1:1632988/93‑1 (MQ=255)
                                                                                ggCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTaa         >  2:337381/1‑93 (MQ=255)
                                                                                  gggttACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTAAcc       <  1:7751343/89‑1 (MQ=255)
                                                                                  ggggTACCGCGCCCCCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTAAcc       <  2:1767076/92‑1 (MQ=255)
                                                                                   gggTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTAACCa      >  1:6495982/1‑93 (MQ=255)
                                                                                   gggTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTAACCa      >  2:3468957/1‑93 (MQ=255)
                                                                                    ggTACTGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCATCTGGCGCGCCAACATTAAGTGATCGGTCTTCTCTAAACTAACCAt     >  1:8944628/1‑93 (MQ=255)
                                                                                    ggTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTAACCAt     >  1:2001156/1‑93 (MQ=255)
                                                                                    ggTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTAACCAt     >  1:1303018/1‑93 (MQ=255)
                                                                                     gTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTAACCATg    >  2:9069364/1‑93 (MQ=255)
                                                                                      tACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTAACCATGa   >  1:9239512/1‑93 (MQ=255)
                                                                                          |                                                                                         
GACAGCGCAAAGTCGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAATCCTTCAGGCGGGTACCACGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTAACCATGAC  >  NC_000913/3566786‑3566964

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: