Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,362,543 T→C V701A (GTA→GCA)  yhcD → putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,362,5430TC99.6% 813.6 / ‑4.6 230V701A (GTA→GCA) yhcDputative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (120/109);  total (120/110)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.30e-01

Only 100 of 240 total aligned reads displayed.

CATCAGGAGTTGAGATCGATACCGCCAGCTTGCCAAAAAATGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCCGTGGTGCAGTAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGACGGCAAAAAACTGCCCAGAGGCGTTG  >  NC_000913/3362453‑3362632
                                                                                          |                                                                                         
cATCAGGAGTTGAGATCGATACCGCCAGCTTGCCAAAAAATGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCCGTGGTGcagcag                                                                                         >  2:3191740/1‑93 (MQ=255)
  tCAGGAGTTGAGATCGATACCGCCAGCTTGCCAAAAAATGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCCGTGGTGCAGCAggg                                                                                       <  1:3725317/93‑1 (MQ=255)
  tCAGGAGTTGAGATCGATACCGCCAGCTTGCCAAAAAATGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCCGTGGTGCAGCAggg                                                                                       >  2:3617202/1‑93 (MQ=255)
   cAGGAGTTGAGATCGATACCGCCAGCTTGCCAAAAAATGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCCGTGGTGCAGCAGGGa                                                                                      >  1:5760694/1‑93 (MQ=255)
   cAGGAGTTGAGATCGATACCGCCAGCTTGCCAAAAAATGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCCGTGGTGCAGCAGGGa                                                                                      <  1:3896444/93‑1 (MQ=255)
    aGGAGTTGAGATCGATACCGCCAGCTTGCCAAAAAATGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCCGTGGTGCAGCAGGGaa                                                                                     <  2:5224352/93‑1 (MQ=255)
      gAGTTGAGATCGATACCGCCAGCTTGCCAAAAAATGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCCGTGGTGCAGCAGGGAAAg                                                                                   <  2:6348081/93‑1 (MQ=255)
       aGTTGAGATCGATCCCGCCAGCTTGCCAAAAAATGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGt                                                                                  <  2:5913684/93‑1 (MQ=255)
       aGTTGAGATCGATACCGCCAGCTTGCCAAAAAATGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGt                                                                                  <  1:5312933/93‑1 (MQ=255)
        gTTGAGATCGATACCGCCAGCTTGCCAAAAAATGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTc                                                                                 >  1:3992840/1‑93 (MQ=255)
           gagaTCGATACCGCCAGCTTGCCAAAAAATGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGa                                                                              <  2:2241240/93‑1 (MQ=255)
            agaTCGATACCGCCAGCTTGCCAAAAAATGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGAt                                                                             >  1:5188270/1‑93 (MQ=255)
             gaTCGATACCGCCAGCTTGCCAAAAAATGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGAtt                                                                            >  1:6031695/1‑93 (MQ=255)
             gaTCGATACCGCCAGCTTGCCAAAAAATGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGAtt                                                                            >  2:831569/1‑93 (MQ=255)
             gaTCGATACCGCCAGCTTGCCAAAAAATGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGAtt                                                                            <  2:930228/93‑1 (MQ=255)
             gaTCGATACCGCCAGCTTGCCAAAAAATGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGAtt                                                                            >  1:3450310/1‑93 (MQ=255)
              aTCGATACCGCCAGCTTGCCAAAAAATGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTc                                                                           <  1:283124/93‑1 (MQ=255)
              aTCGATACCGCCAGCTTGCCAAAAAATGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTc                                                                           >  2:2914710/1‑93 (MQ=255)
               tCGATACCGCCAGCTTGCCAAAAAATGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCa                                                                          <  1:673080/93‑1 (MQ=255)
               tCGATACCGCCAGCTTGCCAAAAAATGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCa                                                                          >  2:5118198/1‑93 (MQ=255)
                cGATACCGCCAGCTTGCCAAAAAATGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAg                                                                         >  2:3650802/1‑93 (MQ=255)
                 gATACCGCCAGCTTGCCAAAAAATGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGt                                                                        >  1:4041467/1‑93 (MQ=255)
                 gATACCGCCAGCTTGCCAAAAAATGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGt                                                                        <  2:5777839/93‑1 (MQ=255)
                 gATACCGCCAGCTTGCCAAAAAATGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCAGATTCAGt                                                                        >  2:3927954/1‑93 (MQ=255)
                  aTACCGCCAGCTTGCCAAAAAATGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTg                                                                       >  1:5447421/1‑93 (MQ=255)
                  aTACCGCCAGCTTGCCAAAAAATGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTg                                                                       >  1:3099679/1‑93 (MQ=255)
                   tACCGCCAGCTTGCCAAAAAATGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGc                                                                      >  2:5386480/1‑93 (MQ=255)
                    aCCGCCAGCTTGCCAAAAAATGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCg                                                                     >  2:4653991/1‑93 (MQ=255)
                     ccGCCAGCTTGCCAAAAAATGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGa                                                                    >  1:1755689/1‑93 (MQ=255)
                      cGCCAGCTTGCCAAAAAATGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGAt                                                                   <  1:1692228/93‑1 (MQ=255)
                      cGCCAGCTTGCCAAAAAATGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGAt                                                                   <  2:2412760/93‑1 (MQ=255)
                       gCCAGCTTGCCAAAAAATGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATa                                                                  <  2:792736/93‑1 (MQ=255)
                        ccAGCTTGCCAAAAAATGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATaa                                                                 <  1:759473/93‑1 (MQ=255)
                        ccAGCTTGCCAAAAAATGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATaa                                                                 <  1:2214922/93‑1 (MQ=255)
                         cAGCTTGCCAAAAAATGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAac                                                                >  1:5912021/1‑93 (MQ=255)
                         cAGCTTGCCAAAAAATGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAac                                                                <  1:1095304/93‑1 (MQ=255)
                           gCTTGCCAAAAAATGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAacac                                                              <  2:2749329/93‑1 (MQ=255)
                           gCTTGCCAAAAAATGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAacac                                                              >  2:3015526/1‑93 (MQ=255)
                              tGCCAAAAAATGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAAc                                                           <  1:5420574/93‑1 (MQ=255)
                              tGCCAAAAAATGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAAc                                                           <  2:4991847/93‑1 (MQ=255)
                              tGCCAAAAAATGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAAc                                                           >  1:3554239/1‑93 (MQ=255)
                                ccAAAAAATGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGc                                                         <  1:5733336/93‑1 (MQ=255)
                                ccAAAAAATGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGACGGATTCAGTGCGATAACACAACGc                                                         >  2:6241476/1‑93 (MQ=255)
                                 cAAAAAATGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGcc                                                        <  2:1423806/93‑1 (MQ=255)
                                 cAAAAAATGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGcc                                                        >  1:3872244/1‑93 (MQ=255)
                                  aaaaaaTGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCg                                                       <  1:4389229/93‑1 (MQ=255)
                                   aaaaaTGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGt                                                      <  1:6375515/93‑1 (MQ=255)
                                    aaaaTGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTg                                                     >  2:3988929/1‑93 (MQ=255)
                                      aaTGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCCGTGGTGCAGTAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCa                                                   <  1:3600847/93‑1 (MQ=255)
                                      aaTGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCa                                                   >  2:1760927/1‑93 (MQ=255)
                                       aTGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCAt                                                  <  1:5962103/93‑1 (MQ=255)
                                       aTGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCAt                                                  <  1:4305573/93‑1 (MQ=255)
                                       aTGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCAt                                                  >  1:3266843/1‑93 (MQ=255)
                                        tGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCAtt                                                 <  2:4513491/93‑1 (MQ=255)
                                            gATATCGGTAACGGCACAAAAATGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTc                                             <  1:3516498/93‑1 (MQ=255)
                                              tatCGGTAACGGCACAAAAATGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCaa                                           >  2:2298263/1‑93 (MQ=255)
                                               atCGGTAACGGCACAAAAATGATCGAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAat                                          <  1:1000242/93‑1 (MQ=255)
                                                tCGGTAACGGCACAAAAATGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAata                                         <  1:5592579/93‑1 (MQ=255)
                                                 cGGTAACGGCACAAAAATGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAatat                                        <  2:3036191/93‑1 (MQ=255)
                                                   gTAACGGCACAAAAATGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATca                                      >  2:1325198/1‑93 (MQ=255)
                                                   gTAACGGCACAAAAATGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATca                                      >  2:4456532/1‑93 (MQ=255)
                                                    tAACGGCACAAAAATGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATcac                                     <  1:1451687/93‑1 (MQ=255)
                                                     aaCGGCACAAAAATGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATcaca                                    >  2:4833257/1‑93 (MQ=255)
                                                      aCGGCACAAAAATGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCacac                                   >  2:799298/1‑93 (MQ=255)
                                                       cGGCACAAAAATGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACt                                  <  1:2599736/93‑1 (MQ=255)
                                                        ggCACAAAAATGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACtt                                 <  1:1887819/93‑1 (MQ=255)
                                                        ggCACAAAAATGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACtt                                 <  1:897520/93‑1 (MQ=255)
                                                         gCACAAAAATGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACttt                                <  2:5603858/93‑1 (MQ=255)
                                                            caAAAATGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCg                             >  2:1903004/1‑93 (MQ=255)
                                                              aaaaTGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCggc                           >  1:3309823/1‑91 (MQ=255)
                                                              aaaaTGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCggc                           >  2:6210267/1‑91 (MQ=255)
                                                              aaaaTGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCggc                           <  2:1978827/93‑3 (MQ=255)
                                                              aaaaTGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCggc                           <  1:1171135/93‑3 (MQ=255)
                                                                aaTGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCcggcgg                         <  2:3069077/93‑1 (MQ=255)
                                                                aaTGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCcggcgg                         >  2:3412785/1‑93 (MQ=255)
                                                                 aTGATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCggcggc                        >  2:244787/1‑93 (MQ=255)
                                                                   gATCAAACAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCaa                      >  1:1573261/1‑93 (MQ=255)
                                                                       aaaCAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCaaaaaa                  >  1:4611845/1‑93 (MQ=255)
                                                                        aaCAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAAc                 >  2:1087685/1‑93 (MQ=255)
                                                                         aCAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACt                <  1:4124438/93‑1 (MQ=255)
                                                                         aCAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACt                <  1:2055427/93‑1 (MQ=255)
                                                                          cAAGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTg               <  1:2660189/93‑1 (MQ=255)
                                                                           aaGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGc              >  1:2467171/1‑93 (MQ=255)
                                                                            aGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGcc             >  2:1214544/1‑93 (MQ=255)
                                                                            aGGCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGcc             >  2:1270203/1‑93 (MQ=255)
                                                                             ggCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGccc            >  1:2579528/1‑93 (MQ=255)
                                                                              gCCGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCa           <  2:4202988/93‑1 (MQ=255)
                                                                               ccGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCag          >  2:1553367/1‑93 (MQ=255)
                                                                               ccGTGGTGCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCag          >  1:4831278/1‑93 (MQ=255)
                                                                                 gtggtgCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCagag        <  1:520999/93‑1 (MQ=255)
                                                                                 gtggtgCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCagag        >  1:5677050/1‑93 (MQ=255)
                                                                                 gtggtgCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCagag        <  1:6176786/93‑1 (MQ=255)
                                                                                  tggtgCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAgg       >  1:2131514/1‑93 (MQ=255)
                                                                                  tggcgCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAgg       >  1:3260441/1‑93 (MQ=255)
                                                                                   ggtgCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGc      <  2:237313/93‑1 (MQ=255)
                                                                                     tgCAGCAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGt    >  2:5961488/1‑93 (MQ=255)
                                                                                      gcagcaGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGtt   >  2:4950023/1‑93 (MQ=255)
                                                                                      gcagcaGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGtt   <  1:6125166/93‑1 (MQ=255)
                                                                                       cagcagGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTg  >  1:5672107/1‑93 (MQ=255)
                                                                                       cagcagGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGGCGGCAAAAAACTGCCCAGAGGCGTTg  <  1:4366577/93‑1 (MQ=255)
                                                                                          |                                                                                         
CATCAGGAGTTGAGATCGATACCGCCAGCTTGCCAAAAAATGTCGATATCGGTAACGGCACAAAAATGATCAAACAAGGCCGTGGTGCAGTAGGGAAAGTCGGATTCAGTGCGATAACACAACGCCGTGCATTACTCAATATCACACTTTCCGACGGCAAAAAACTGCCCAGAGGCGTTG  >  NC_000913/3362453‑3362632

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