Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,339,624 A→G K348R (AAG→AGG)  yciM → putative heat shock protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,339,6240AG99.6% 714.8 / ‑4.7 222K348R (AAG→AGG) yciMputative heat shock protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (107/114);  total (108/114)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.30e-01

Only 100 of 226 total aligned reads displayed.

CCATAAGTTAATGGATTACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAAGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTGTCGGGCCTGGT  >  NC_000913/1339535‑1339713
                                                                                         |                                                                                         
 cATAAGTTAATGGATTACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTa                                                                                       >  2:5852890/1‑93 (MQ=255)
 cATAAGTTAATGGATTACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTa                                                                                       >  1:1517658/1‑93 (MQ=255)
  aTAAGTTAATGGATTACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTAc                                                                                      >  2:2559592/1‑93 (MQ=255)
   tAAGTTAATGGATTACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACg                                                                                     >  2:6126764/1‑93 (MQ=255)
   tAAGTTAATGGATTACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACg                                                                                     >  2:5079376/1‑93 (MQ=255)
   tAAGTTAATGGATTACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACg                                                                                     >  2:4339015/1‑93 (MQ=255)
    aaGTTAATGGATTACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACgt                                                                                    >  1:1154563/1‑93 (MQ=255)
    aaGTTAATGGATTACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACgt                                                                                    >  1:4248837/1‑93 (MQ=255)
    aaGTTAATGGATTACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGACGAGAGGGTACgt                                                                                    >  1:5317112/1‑93 (MQ=255)
     aGTTAATGGATTACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACgta                                                                                   >  1:3583498/1‑93 (MQ=255)
       ttAATGGATTACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACgtagt                                                                                 <  1:705057/93‑1 (MQ=255)
         aaTGGATTACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTaa                                                                               >  1:3844600/1‑93 (MQ=255)
         aaTGGATTACCACTTAAATGAAGCGGAAGAAGGGCGGGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTaa                                                                               >  1:958576/1‑93 (MQ=255)
           tGGATTACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGc                                                                             >  2:3663042/1‑93 (MQ=255)
             gATTACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCt                                                                           <  2:5107948/93‑1 (MQ=255)
             gATTACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCt                                                                           <  1:2936696/93‑1 (MQ=255)
             gATTACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCt                                                                           >  1:3207690/1‑93 (MQ=255)
             gATTACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCt                                                                           >  2:5484852/1‑93 (MQ=255)
              aTTACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTc                                                                          >  1:2270655/1‑93 (MQ=255)
              aTTACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTc                                                                          >  2:6552644/1‑93 (MQ=255)
              aTTACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTc                                                                          <  2:445885/93‑1 (MQ=255)
               ttACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCg                                                                         <  1:1938644/93‑1 (MQ=255)
               ttACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCg                                                                         <  2:2392200/93‑1 (MQ=255)
               ttACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGAGTACGTAGTAAGCCTCg                                                                         >  1:4768313/1‑93 (MQ=255)
                tACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGt                                                                        <  1:5841930/93‑1 (MQ=255)
                 aCCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGtt                                                                       <  2:139124/93‑1 (MQ=255)
                 aCCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGta                                                                       >  2:3643002/1‑92 (MQ=255)
                  ccACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTa                                                                      >  1:1749512/1‑93 (MQ=255)
                  ccACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTa                                                                      >  2:181154/1‑93 (MQ=255)
                  ccACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTa                                                                      <  1:1424331/93‑1 (MQ=255)
                   cACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTaa                                                                     <  1:4139943/93‑2 (MQ=255)
                   cACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTAt                                                                     <  1:2003987/93‑1 (MQ=255)
                    aCTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGGTGCGTGACATGGTTGGGGAGAGGGTACGTAGTAATCCTCGTTATc                                                                    <  1:2940296/93‑1 (MQ=255)
                    aCTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATc                                                                    >  2:6634328/1‑93 (MQ=255)
                    aCTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATc                                                                    >  2:1254940/1‑93 (MQ=255)
                     cTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCg                                                                   <  2:4590301/93‑1 (MQ=255)
                     cTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCg                                                                   >  1:6357774/1‑93 (MQ=255)
                      ttAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGt                                                                  <  1:2066420/93‑1 (MQ=255)
                        aaaTGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTg                                                                >  1:489114/1‑93 (MQ=255)
                          aTGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTTGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGcc                                                              <  1:6501763/93‑1 (MQ=255)
                            gAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAg                                                            >  2:6279922/1‑93 (MQ=255)
                              aGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGaa                                                          >  1:779379/1‑93 (MQ=255)
                               gCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGaaa                                                         >  2:3889124/1‑93 (MQ=255)
                               gCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGaaa                                                         >  1:5980044/1‑93 (MQ=255)
                               gCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGaaa                                                         <  1:5519073/93‑1 (MQ=255)
                                gggAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAAt                                                        <  2:2895515/92‑1 (MQ=255)
                                cGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAAt                                                        >  2:4798309/1‑93 (MQ=255)
                                cGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAAt                                                        <  1:4566469/93‑1 (MQ=255)
                                  gaagaaGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAAtgt                                                      <  2:6581582/93‑1 (MQ=255)
                                  gaagaaGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAAtgt                                                      >  2:2919619/1‑93 (MQ=255)
                                   aagaagGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAAtgtg                                                     <  1:6216620/93‑1 (MQ=255)
                                     gaagGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGt                                                   <  1:5099960/93‑1 (MQ=255)
                                     gaagGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGt                                                   >  2:4266813/1‑93 (MQ=255)
                                      aagGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGtt                                                  <  2:3535339/93‑1 (MQ=255)
                                       agGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGttt                                                 >  2:3969489/1‑93 (MQ=255)
                                       agGGAGTGCGATAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGttt                                                 <  2:3283909/93‑1 (MQ=255)
                                        gggCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGtttt                                                <  1:4714991/93‑1 (MQ=255)
                                          gCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTAc                                              >  2:5726023/1‑93 (MQ=255)
                                               ccAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCAt                                         <  1:5683583/93‑1 (MQ=255)
                                                cAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATa                                        <  2:5188275/93‑1 (MQ=255)
                                                 aaaGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATac                                       >  1:2597096/1‑93 (MQ=255)
                                                  aaGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATaca                                      >  1:1007625/1‑93 (MQ=255)
                                                   agagagCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATacac                                     <  2:5189624/93‑1 (MQ=255)
                                                    gagagCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACAcc                                    >  2:260366/1‑93 (MQ=255)
                                                    gagagCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACAcc                                    <  1:5018587/93‑1 (MQ=255)
                                                     agagCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATAAAccc                                   >  2:5091556/1‑93 (MQ=255)
                                                     agagCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGGACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGGTTTACCGCATACAccc                                   >  1:6458394/1‑93 (MQ=255)
                                                      gagCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCct                                  >  2:1709306/1‑93 (MQ=255)
                                                        gCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCctct                                <  2:4833114/93‑1 (MQ=255)
                                                         ccTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTa                               <  1:4303153/93‑1 (MQ=255)
                                                         ccTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTa                               <  1:4144809/93‑1 (MQ=255)
                                                          cTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCTTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTAc                              >  2:5673020/1‑93 (MQ=255)
                                                           tGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACt                             <  2:5363292/93‑1 (MQ=255)
                                                           tGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACt                             >  1:6157553/1‑93 (MQ=255)
                                                           tGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACt                             <  2:6496170/93‑1 (MQ=255)
                                                           tGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACt                             >  1:5299074/1‑93 (MQ=255)
                                                           tGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACAGCATACACCCTCTACt                             >  1:1254829/1‑93 (MQ=255)
                                                           tGATGGTGCGGCCTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACt                             <  2:4179726/93‑1 (MQ=255)
                                                            gATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACCCCCTCTACTg                            <  2:6370712/93‑1 (MQ=255)
                                                            gATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTg                            <  1:2433495/93‑1 (MQ=255)
                                                             aTGGTGCTGCGTGACATGGTTGGCGAGAGGGTAAGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTgg                           >  1:1877965/1‑93 (MQ=255)
                                                              tGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGc                          >  2:336682/1‑93 (MQ=255)
                                                               ggTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCa                         <  2:4710246/93‑1 (MQ=255)
                                                                  gctgcGTGACAGGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCGCTACTGGCATTg                      <  1:353179/93‑1 (MQ=255)
                                                                   ctgcGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGt                     >  2:1161732/1‑93 (MQ=255)
                                                                    tgcGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTc                    >  2:6431604/1‑93 (MQ=255)
                                                                    tgcGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTc                    <  2:2500780/93‑1 (MQ=255)
                                                                    tgcGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTc                    <  1:2184622/93‑1 (MQ=255)
                                                                       gTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGt                 >  1:4656443/1‑93 (MQ=255)
                                                                         gACATGGTTGTCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCt               <  2:798380/93‑1 (MQ=255)
                                                                         gACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCt               <  2:6171664/93‑1 (MQ=255)
                                                                         gACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCt               <  2:1572349/93‑1 (MQ=255)
                                                                          aCATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCtt              >  2:4349236/1‑93 (MQ=255)
                                                                           cATGGTTGGCGAGAGGGTACGTCGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTg             <  2:1459626/93‑1 (MQ=255)
                                                                           cATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCTACCACCCTCTCCTGGCATTGTCCGTCTTg             >  2:1925247/1‑93 (MQ=255)
                                                                            aTGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTGt            <  2:4054048/93‑1 (MQ=255)
                                                                             tGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTGTc           <  2:4985518/93‑1 (MQ=255)
                                                                                   gCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTGTCGGGCCt     <  2:682432/93‑1 (MQ=255)
                                                                                    cGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTGTCGGGCCTg    <  2:5287700/93‑1 (MQ=255)
                                                                                     gagaGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTGTCGGGCCTgg   >  1:6641880/1‑93 (MQ=255)
                                                                                         |                                                                                         
CCATAAGTTAATGGATTACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAAGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTGTCGGGCCTGGT  >  NC_000913/1339535‑1339713

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