Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,883,067 T→C D121D (GAT→GAC yidC → 60 KD inner‑membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,883,0670TC99.4% 990.3 / ‑4.1 314D121D (GAT→GACyidC60 KD inner‑membrane protein
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base C (146/166);  total (146/168)
Fisher's exact test for biased strand distribution p-value = 5.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 319 total aligned reads displayed.

GAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGATAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCAAAACGAACTGCAGGTGCCGAT  >  NC_000913/3882977‑3883156
                                                                                          |                                                                                         
gAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACaa                                                                                         <  2:267449/93‑1 (MQ=255)
  aCTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAAcc                                                                                       >  2:2416922/1‑93 (MQ=255)
   ctctACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAAccc                                                                                      >  1:4946783/1‑93 (MQ=255)
    tctACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCg                                                                                     >  2:3320380/1‑93 (MQ=255)
    tctACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCg                                                                                     >  1:5671074/1‑93 (MQ=255)
    tctACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCg                                                                                     >  2:5171135/1‑93 (MQ=255)
     ctACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCgg                                                                                    <  1:2024328/93‑1 (MQ=255)
      tACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGc                                                                                   <  2:5493727/93‑1 (MQ=255)
       aCCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCt                                                                                  >  1:5712772/1‑93 (MQ=255)
       aCCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCt                                                                                  <  2:4465391/93‑1 (MQ=255)
         ccAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTaa                                                                                <  2:2919974/93‑1 (MQ=255)
          cAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAAc                                                                               <  1:4349604/93‑1 (MQ=255)
           aGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACg                                                                              >  1:87697/1‑93 (MQ=255)
            gCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACgg                                                                             <  1:1028895/93‑1 (MQ=255)
             ccGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGc                                                                            <  2:2724228/93‑1 (MQ=255)
              cGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGcc                                                                           <  2:2843855/93‑1 (MQ=255)
              cGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGcc                                                                           >  2:3571958/1‑93 (MQ=255)
               gTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGccc                                                                          >  1:2927419/1‑93 (MQ=255)
               gTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGccc                                                                          <  2:4054711/93‑1 (MQ=255)
                 tCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCcgc                                                                        <  1:2495295/93‑1 (MQ=255)
                  ccAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCcgcg                                                                       <  1:4964028/93‑1 (MQ=255)
                   cAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGt                                                                      >  2:2809568/1‑93 (MQ=255)
                    aGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTc                                                                     <  2:5381027/93‑1 (MQ=255)
                    aGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTc                                                                     >  1:773944/1‑93 (MQ=255)
                     gCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTcc                                                                    >  2:6706655/1‑93 (MQ=255)
                       tGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGc                                                                  >  2:1535097/1‑93 (MQ=255)
                        gTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCt                                                                 >  1:5855595/1‑93 (MQ=255)
                          tGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGt                                                               <  1:2128944/93‑1 (MQ=255)
                          tGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGt                                                               <  2:5594354/93‑1 (MQ=255)
                           ggAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGta                                                              <  1:4182008/93‑1 (MQ=255)
                             aaaCTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCAGGCTAACGGCCCGCGCCCGCTGcata                                                            >  2:3434749/1‑93 (MQ=255)
                              aaCTTCACCGCAGTTTATTTATCAGGCACAGAGGGGTCTGACCGGTCGTGATGGCCCGGCCAACCCGGCTAACGGCCCGCGTCCGCTGTATaa                                                           <  1:5263952/93‑1 (MQ=255)
                              aaCTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGtatta                                                           >  1:6118908/1‑91 (MQ=255)
                              aaCTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATaa                                                           >  1:530297/1‑93 (MQ=255)
                              aaCTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATaa                                                           <  1:997720/93‑1 (MQ=255)
                               aCTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAAc                                                          <  2:5555204/93‑1 (MQ=255)
                                cTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGGCAACCCGGCTAACGGCCCGCGTCCGCTGTATGACg                                                         >  1:4479367/1‑93 (MQ=255)
                                 ttCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGt                                                        >  1:6630306/1‑93 (MQ=255)
                                 ttCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGt                                                        <  2:2112587/93‑1 (MQ=255)
                                 ttCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGt                                                        <  2:1632151/93‑1 (MQ=255)
                                 ttCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGt                                                        <  1:4669488/93‑1 (MQ=255)
                                 ttCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGt                                                        <  2:5973535/93‑1 (MQ=255)
                                 ttCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGt                                                        <  1:2897271/93‑1 (MQ=255)
                                 ttCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGt                                                        <  2:5436679/93‑1 (MQ=255)
                                     ccGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGaa                                                    >  1:2301114/1‑93 (MQ=255)
                                      cGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGaaa                                                   <  1:6654953/93‑1 (MQ=255)
                                      cGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGaaa                                                   >  2:1289131/1‑93 (MQ=255)
                                       gCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGaaaa                                                  >  1:4810858/1‑93 (MQ=255)
                                       gCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGaaaa                                                  <  2:6455395/93‑1 (MQ=255)
                                          gTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGa                                               >  1:4004430/1‑93 (MQ=255)
                                            ttatttatCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGAct                                             >  2:2453281/1‑92 (MQ=255)
                                              atttatCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCt                                           <  1:5157754/93‑1 (MQ=255)
                                              atttatCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCt                                           <  2:3479616/93‑1 (MQ=255)
                                              atttatCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCt                                           >  2:3143383/1‑93 (MQ=255)
                                               tttatCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCtt                                          <  2:1339226/93‑1 (MQ=255)
                                               tttatCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCtt                                          >  2:849990/1‑93 (MQ=255)
                                                ttatCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTa                                         <  2:4631543/93‑1 (MQ=255)
                                                ttatCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTa                                         <  2:1273503/93‑1 (MQ=255)
                                                 tatCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTAt                                        <  2:6343568/93‑1 (MQ=255)
                                                 tatCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTAt                                        <  2:3864175/93‑1 (MQ=255)
                                                  atCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTAtg                                       <  1:1191804/93‑1 (MQ=255)
                                                   tCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTAtgt                                      <  1:3814625/93‑1 (MQ=255)
                                                     aGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTgc                                    <  2:690620/93‑1 (MQ=255)
                                                      ggCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTgct                                   >  1:1495504/1‑93 (MQ=255)
                                                         acaGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTgctggc                                >  1:6194166/1‑93 (MQ=255)
                                                         acaGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTgctggc                                <  2:6399395/93‑1 (MQ=255)
                                                          caGAGCGGTCTGACCGGTCGTGATGGCCCGGATAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTgctggct                               <  2:725105/93‑1 (MQ=255)
                                                           agagCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTgctggctg                              <  1:627587/93‑1 (MQ=255)
                                                            gagCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGa                             >  2:3621877/1‑93 (MQ=255)
                                                             agCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGaa                            <  2:4414506/93‑1 (MQ=255)
                                                              gCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAg                           <  1:1362630/93‑1 (MQ=255)
                                                               cGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAgg                          <  2:508573/93‑1 (MQ=255)
                                                                ggTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGt                         <  2:5700597/93‑1 (MQ=255)
                                                                 gTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTc                        <  2:5260521/93‑1 (MQ=255)
                                                                  tCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCa                       >  2:6649572/1‑93 (MQ=255)
                                                                   cTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCaa                      <  1:3375706/93‑1 (MQ=255)
                                                                    tGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCaaa                     <  2:1919545/93‑1 (MQ=255)
                                                                     gACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCaaaa                    >  2:4953447/1‑93 (MQ=255)
                                                                     gACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCaaaa                    >  1:4643908/1‑93 (MQ=255)
                                                                     gACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCaaaa                    <  2:3039265/93‑1 (MQ=255)
                                                                      aCCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCAAAAc                   >  1:3131013/1‑93 (MQ=255)
                                                                      aCCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCAAAAc                   <  1:2221262/93‑1 (MQ=255)
                                                                       ccGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAATGTCAAAAca                  >  1:4572119/1‑92 (MQ=255)
                                                                       ccGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCAAAACg                  <  2:4183218/93‑1 (MQ=255)
                                                                       ccGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCAAAACg                  >  1:1957278/1‑93 (MQ=255)
                                                                          gTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCAAAACGAAc               <  2:1826544/93‑1 (MQ=255)
                                                                           tCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCAAAAAGAACt              >  1:2704775/1‑93 (MQ=255)
                                                                            cGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCAAAACGAACTg             >  1:1605546/1‑93 (MQ=255)
                                                                            cGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCAAAACGAACTg             <  2:34131/93‑1 (MQ=255)
                                                                             gTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCAAAACGAACTGc            >  2:4828424/1‑93 (MQ=255)
                                                                              tGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCAAAACGAACTGCa           >  1:1766469/1‑93 (MQ=255)
                                                                              tGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCAAAACGAACTGCa           <  1:6540679/93‑1 (MQ=255)
                                                                               gATGGCCCGGACACCCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCAAAACGAACTGCAg          <  1:3688318/93‑1 (MQ=255)
                                                                               gATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCAAAACGAACTGCAg          >  1:4754447/1‑93 (MQ=255)
                                                                                aTGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGGATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGCCAAAACGAACTGCAgg         >  1:5951227/1‑93 (MQ=255)
                                                                                 tGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCAAAACGAACTGCAGGt        <  1:6431473/93‑1 (MQ=255)
                                                                                    cccGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCAAAACGAACTGCAGGTGcc     >  2:1183754/1‑93 (MQ=255)
                                                                                      cGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCAAAACGAACTGCAGGTGCCGa   >  2:5818682/1‑93 (MQ=255)
                                                                                       ggACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCAAAACGAACTGCAGGTGCCGAt  <  1:4700129/93‑1 (MQ=255)
                                                                                       ggACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCAAAACGAACTGCAGGTGCCGAt  >  1:5611164/1‑93 (MQ=255)
                                                                                          |                                                                                         
GAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGATAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCAAAACGAACTGCAGGTGCCGAT  >  NC_000913/3882977‑3883156

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