Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 293 A→G intergenic (+38/‑44) thrL → / → thrA thr operon leader peptide/aspartokinase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132930AG99.6% 735.0 / ‑4.9 238intergenic (+38/‑44)thrL/thrAthr operon leader peptide/aspartokinase I
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (98/139);  total (99/139)
Fisher's exact test for biased strand distribution p-value = 4.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.91e-01

Only 100 of 262 total aligned reads displayed.

GCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAAC  >  NC_000913/203‑382
                                                                                           |                                                                                          
gCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGt                                                                                          <  2:2040388/93‑1 (MQ=255)
 caccacCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTg                                                                                         <  1:400342/93‑1 (MQ=255)
 caccacCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTg                                                                                         <  1:2225457/93‑1 (MQ=255)
  accaccATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGc                                                                                        >  2:5131445/1‑93 (MQ=255)
   ccaccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCg                                                                                       <  1:3562965/93‑1 (MQ=255)
   ccaccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCg                                                                                       >  1:1365291/1‑93 (MQ=255)
   ccaccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCg                                                                                       >  2:3781340/1‑93 (MQ=255)
    caccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCgg                                                                                      >  2:5991927/1‑93 (MQ=255)
     cccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCggg                                                                                     <  1:2354554/92‑1 (MQ=255)
       caTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCt                                                                                  >  1:5087173/1‑93 (MQ=255)
        aTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCtt                                                                                 <  1:7717204/93‑1 (MQ=255)
         ttACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAACAAATCAAGCACCTGACGGCGCGGGCttt                                                                                >  1:3086817/1‑93 (MQ=255)
         ttACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCttt                                                                                <  2:2415184/93‑1 (MQ=255)
          tACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCtttt                                                                               >  2:2290379/1‑93 (MQ=255)
          tACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCtttt                                                                               >  1:3862779/1‑93 (MQ=255)
            ccaccaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCtttttt                                                                             >  1:4854630/1‑93 (MQ=255)
             caccaccaTCACCATTACCACAGGTAACGGTGCGGGCTGGCGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCttttttt                                                                            >  1:7877444/1‑93 (MQ=255)
               ccaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCttttttttt                                                                          >  1:6564902/1‑93 (MQ=255)
                  ccatcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCgc                                                                       >  1:6079442/1‑92 (MQ=255)
                     tcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCa                                                                    <  1:4673897/93‑1 (MQ=255)
                     tcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCa                                                                    >  1:6384848/1‑93 (MQ=255)
                       accatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCaaa                                                                  <  2:1676386/93‑1 (MQ=255)
                        ccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAg                                                                 <  1:593105/93‑1 (MQ=255)
                        ccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAg                                                                 >  1:1082199/1‑93 (MQ=255)
                        ccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAg                                                                 >  2:728682/1‑93 (MQ=255)
                           ttACCACAGGTAACGGTGCGGGCTGGCGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTa                                                              <  1:1923487/93‑1 (MQ=255)
                           ttACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTTGACCAAAGGTa                                                              >  1:3344357/1‑93 (MQ=255)
                           ttACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTa                                                              >  2:1035501/1‑93 (MQ=255)
                           ttACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTa                                                              >  2:343075/1‑93 (MQ=255)
                           ttACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTa                                                              >  1:4961146/1‑93 (MQ=255)
                           ttACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTa                                                              >  1:1787068/1‑93 (MQ=255)
                            tACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTaa                                                             <  2:4747111/93‑1 (MQ=255)
                            tACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTaa                                                             <  2:5536170/93‑1 (MQ=255)
                             aCCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAAc                                                            <  2:5207932/93‑1 (MQ=255)
                               cacaGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGa                                                          <  2:2789842/93‑1 (MQ=255)
                               caagggTAACGGGGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGa                                                          <  2:1660034/89‑1 (MQ=255)
                                acaGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAg                                                         >  1:903207/1‑93 (MQ=255)
                                acaGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAg                                                         >  2:7652203/1‑93 (MQ=255)
                                 caGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAgg                                                        <  1:5255002/93‑1 (MQ=255)
                                   ggTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTa                                                      >  2:4705264/1‑93 (MQ=255)
                                   ggTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTa                                                      <  1:2611857/93‑1 (MQ=255)
                                     tAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTaac                                                    <  2:951655/93‑1 (MQ=255)
                                       aCGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCCGCACCCACCGGGGGGGGGTTTTTTTTTTTGCCAAAGGGAACAAGGGaaaa                                                   >  2:3012544/1‑91 (MQ=255)
                                       aCGGTGCGGCCTGACGCGTACAGAAAACACAAAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTaacaa                                                  <  2:3325405/93‑1 (MQ=255)
                                         ggTGCGGGCTGACGCGTACAGAAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAAcc                                                <  2:6772434/93‑1 (MQ=255)
                                          gTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCa                                               <  1:2719374/93‑1 (MQ=255)
                                          gTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCa                                               <  2:563113/93‑1 (MQ=255)
                                             cGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGc                                            <  2:7683189/93‑1 (MQ=255)
                                                gCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAg                                         >  2:4891098/1‑93 (MQ=255)
                                                gCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAg                                         >  2:3453888/1‑93 (MQ=255)
                                                 cTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGGGCGGGCTTTTTTTTTCGACCCAAAGGAACGAGGTAACAAACATGCGAgg                                        >  1:4224959/1‑92 (MQ=255)
                                                    aCGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGtt                                     >  1:534367/1‑93 (MQ=255)
                                                      gcgTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGa                                   <  2:1313159/93‑1 (MQ=255)
                                                       cgTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGaa                                  <  1:3349240/93‑1 (MQ=255)
                                                       cgTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGaa                                  <  1:5679828/93‑1 (MQ=255)
                                                       cgTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGaa                                  <  2:1384771/93‑1 (MQ=255)
                                                        gTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGaa                                  <  1:4339005/93‑1 (MQ=255)
                                                        gTACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAg                                 >  1:7520709/1‑93 (MQ=255)
                                                         tACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGt                                >  2:5401282/1‑93 (MQ=255)
                                                         tACAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGt                                <  1:7248719/93‑1 (MQ=255)
                                                           cAGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTc                              <  1:5991928/93‑1 (MQ=255)
                                                            aGGAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcg                             <  2:7333441/93‑1 (MQ=255)
                                                             ggAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGCCCAAAGGGAACCAGTCAACAACCATCAGAGTGTTGACGTTcgg                            >  1:2084478/1‑93 (MQ=255)
                                                             ggAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcgg                            >  1:7390524/1‑93 (MQ=255)
                                                             ggAAACACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcgg                            <  1:3732930/93‑1 (MQ=255)
                                                               aaaCACAGAAAAAAGCCCGCTCCTGCCGGTGCGGGCTTTTTTTTTAGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcggcg                          <  2:6153939/93‑1 (MQ=255)
                                                                aaCACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcggcgg                         >  2:6735760/1‑93 (MQ=255)
                                                                aaCACAGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcggcgg                         >  2:6525867/1‑93 (MQ=255)
                                                                    caGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACa                     >  2:3634050/1‑93 (MQ=255)
                                                                    caGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACa                     <  1:4599303/93‑1 (MQ=255)
                                                                    caGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGGTCGGCCGTACa                     >  1:6289112/1‑93 (MQ=255)
                                                                    caGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGGGTTGAAGTTCGGCGGTACa                     >  1:2477640/1‑93 (MQ=255)
                                                                     aGAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTTGACCCAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGGGGGACAt                    >  1:5830389/1‑93 (MQ=255)
                                                                      gAAAAAAGCCCGCCCCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATc                   <  1:6199482/93‑1 (MQ=255)
                                                                      gAAAAAAGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATc                   <  1:341243/93‑1 (MQ=255)
                                                                      aaaaaaaGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATc                   <  2:3966113/92‑1 (MQ=255)
                                                                      aaaaaaaGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATc                   <  2:6007848/92‑1 (MQ=255)
                                                                      aaaaaaaCCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATc                   <  2:7132161/92‑1 (MQ=255)
                                                                       aaaaaaGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCa                  <  1:6533752/93‑1 (MQ=255)
                                                                       aaaaaaCCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCa                  <  2:6636805/93‑1 (MQ=255)
                                                                        aaaaaGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAg                 <  1:7847283/93‑1 (MQ=255)
                                                                        aaaaaGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAg                 >  2:5455103/1‑93 (MQ=255)
                                                                         aaaaCCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGt                <  2:4531441/93‑1 (MQ=255)
                                                                          aaaGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTg               <  1:1024502/93‑1 (MQ=255)
                                                                          aaaGCCCGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTg               >  2:6461401/1‑93 (MQ=255)
                                                                                ccGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCaaa          <  1:4636258/93‑1 (MQ=255)
                                                                                 cGCCCCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAAt         <  2:1424900/93‑1 (MQ=255)
                                                                                 cGCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGGAACGAGGTAACAACCATGCGAGTGTTGAAGGTCGGCGGTACATCAGTGGCAAAt         >  1:5162809/1‑93 (MQ=255)
                                                                                 cGCACCTGACGGTGCGGGCTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATg        <  1:7151594/93‑1 (MQ=255)
                                                                                  gCACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATg        >  1:1689205/1‑93 (MQ=255)
                                                                                   cACCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGc       <  1:6777545/93‑1 (MQ=255)
                                                                                    cccggACGGTGCGGGGTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCa      <  1:6427235/89‑1 (MQ=255)
                                                                                    aCCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCa      <  2:415484/93‑1 (MQ=255)
                                                                                    aCCTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCa      <  2:7480489/93‑1 (MQ=255)
                                                                                      cTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGa    >  2:4486460/1‑93 (MQ=255)
                                                                                      cTGACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGa    >  2:4029912/1‑93 (MQ=255)
                                                                                       tGACGGTGCGGGCTTTTTTTTTCGACCAAAGGGAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCaaaa   >  2:6566689/1‑93 (MQ=255)
                                                                                       tGACGGTGCGGGCTTTTTTTTTCGACCAAAGGGAACGAGGTAACAACCATGCGAGTGTTGAAATTCGGCGGTACATCAGTGGCAACGGCAGaa   >  2:3850436/1‑93 (MQ=255)
                                                                                        gACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAAc  <  1:3677153/93‑1 (MQ=255)
                                                                                        gACGGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAAc  <  1:7199436/93‑1 (MQ=255)
                                                                                           |                                                                                          
GCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAAC  >  NC_000913/203‑382

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