Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,899,440 T→C T137A (ACC→GCC)  yieC ← putative receptor protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,899,4400TC83.8% 677.2 / 127.0 315T137A (ACC→GCC) yieCputative receptor protein
Reads supporting (aligned to +/- strand):  ref base T (26/25);  new base C (117/146);  total (143/172)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.54e-04

Only 100 of 323 total aligned reads displayed.

TAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGGAAGAATAAA  >  NC_000913/3899351‑3899532
                                                                                          |                                                                                            
 aaCTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTcccc                                                                                          >  1:4346364/1‑93 (MQ=255)
    tGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCaa                                                                                       <  2:2045681/93‑1 (MQ=255)
    tGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCaa                                                                                       >  1:6829728/1‑93 (MQ=255)
     gCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAAc                                                                                      <  1:1304096/93‑1 (MQ=255)
     gCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAAc                                                                                      <  2:6128038/93‑1 (MQ=255)
      cAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAAcc                                                                                     <  2:5892191/93‑1 (MQ=255)
       aaaTCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCa                                                                                    >  2:567556/1‑93 (MQ=255)
       aaaTCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCa                                                                                    >  1:4638087/1‑93 (MQ=255)
        aaTCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAg                                                                                   >  2:4642501/1‑93 (MQ=255)
        aaTCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAg                                                                                   >  2:411402/1‑93 (MQ=255)
          tCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGaa                                                                                 <  2:2890836/93‑1 (MQ=255)
          tCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGaa                                                                                 <  1:7381855/93‑1 (MQ=255)
          tCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGaa                                                                                 >  2:1831751/1‑93 (MQ=255)
          acacaCCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGaa                                                                                 <  1:2587790/89‑1 (MQ=255)
            aaaCCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACg                                                                               >  2:6338241/1‑93 (MQ=255)
            aaaCCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACg                                                                               <  1:2056929/93‑1 (MQ=255)
             aaCCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACgg                                                                              >  2:3122319/1‑93 (MQ=255)
             aaCCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACgg                                                                              >  2:7472859/1‑93 (MQ=255)
              aCCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGt                                                                             >  2:7029839/1‑93 (MQ=255)
              aCCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGt                                                                             >  1:3589554/1‑93 (MQ=255)
               ccAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTa                                                                            <  1:7041696/93‑1 (MQ=255)
                cAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTaa                                                                           <  2:1749099/93‑1 (MQ=255)
                  gCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATa                                                                         >  1:6322870/1‑93 (MQ=255)
                  gCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATa                                                                         <  1:6288653/93‑1 (MQ=255)
                    cGGAGTGTTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGc                                                                       >  1:1957970/1‑93 (MQ=255)
                     ggAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGcc                                                                      <  1:5819268/93‑1 (MQ=255)
                     ggAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGcc                                                                      >  2:6730601/1‑93 (MQ=255)
                       aGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGt                                                                    <  1:7221040/93‑1 (MQ=255)
                        gTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGtt                                                                   >  2:862948/1‑93 (MQ=255)
                        gTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGtt                                                                   <  1:668483/93‑1 (MQ=255)
                         tATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTg                                                                  <  2:5768220/93‑1 (MQ=255)
                          aTTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGt                                                                 <  1:7448668/93‑1 (MQ=255)
                           ttCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTa                                                                >  2:1019213/1‑93 (MQ=255)
                           ttCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTa                                                                >  1:3044047/1‑93 (MQ=255)
                           ttCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTa                                                                >  1:7272395/1‑93 (MQ=255)
                            tCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTaa                                                               >  2:3616699/1‑93 (MQ=255)
                            tCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTaa                                                               >  1:2952607/1‑93 (MQ=255)
                            tCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTaa                                                               <  1:2516794/93‑1 (MQ=255)
                              gTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACt                                                             <  1:3156447/93‑1 (MQ=255)
                               ttACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTa                                                            >  2:5390479/1‑93 (MQ=255)
                               ttACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTa                                                            >  2:2109568/1‑93 (MQ=255)
                                tACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTaa                                                           <  2:4499958/93‑1 (MQ=255)
                                tACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTaa                                                           <  1:6795508/93‑1 (MQ=255)
                                 aCCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTaaa                                                          >  1:3010763/1‑93 (MQ=255)
                                   cAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGc                                                        <  2:7322347/93‑1 (MQ=255)
                                   cAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGc                                                        <  1:3756464/93‑1 (MQ=255)
                                   cAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGc                                                        >  2:5493539/1‑93 (MQ=255)
                                     aaGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCa                                                      >  2:681087/1‑93 (MQ=255)
                                      aGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCaa                                                     >  1:5104146/1‑93 (MQ=255)
                                       gCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAAt                                                    <  2:3537302/93‑1 (MQ=255)
                                       gCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAAt                                                    <  1:3314695/93‑1 (MQ=255)
                                         ggCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAAttt                                                  <  1:1140797/93‑1 (MQ=255)
                                         ggCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAAttt                                                  >  2:6534023/1‑93 (MQ=255)
                                          gCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTc                                                 >  2:1189860/1‑93 (MQ=255)
                                          gCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTc                                                 >  2:4911411/1‑93 (MQ=255)
                                             cAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCAt                                               >  1:2697189/1‑93 (MQ=255)
                                             cAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCGAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATc                                              <  2:299118/93‑1 (MQ=255)
                                             cAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATc                                              <  1:2044559/93‑1 (MQ=255)
                                             cAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATc                                              <  2:4334217/93‑1 (MQ=255)
                                             cAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATc                                              <  2:370570/93‑1 (MQ=255)
                                              agagAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCt                                             >  2:2402400/1‑93 (MQ=255)
                                               gagaACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCtt                                            >  2:1520257/1‑93 (MQ=255)
                                               gagaACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCtt                                            >  1:5550168/1‑93 (MQ=255)
                                                agaACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCttt                                           >  2:5605952/1‑93 (MQ=255)
                                                 gaACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTc                                          <  2:2282484/93‑1 (MQ=255)
                                                   aCCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTcac                                        >  1:3981564/1‑93 (MQ=255)
                                                   aCCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTcac                                        <  1:651932/93‑1 (MQ=255)
                                                    ccAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTcaca                                       >  1:1387272/1‑93 (MQ=255)
                                                       aTCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAAt                                    >  2:2649640/1‑93 (MQ=255)
                                                        tCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAAtt                                   >  2:7136772/1‑93 (MQ=255)
                                                         cGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAAttt                                  >  2:1123888/1‑93 (MQ=255)
                                                           cccATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTAc                                <  1:1546729/93‑1 (MQ=255)
                                                           cccATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTAc                                <  1:5496415/93‑1 (MQ=255)
                                                            ccATGATTTAGGTGAACCATGAGAGGCGGCCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTAc                                <  1:5198781/93‑1 (MQ=255)
                                                            ccATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACt                               <  2:442641/93‑1 (MQ=255)
                                                             cATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTg                              >  1:1773036/1‑93 (MQ=255)
                                                             cATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTg                              >  1:5723150/1‑93 (MQ=255)
                                                              aTGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGa                             >  2:2747792/1‑93 (MQ=255)
                                                              aTGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGa                             <  1:2843350/93‑1 (MQ=255)
                                                              aTGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGa                             <  2:7231847/93‑1 (MQ=255)
                                                                gATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGaaa                           <  2:231848/93‑1 (MQ=255)
                                                                  tttAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATc                         >  1:4049799/1‑93 (MQ=255)
                                                                   ttAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCt                        <  1:1046802/93‑1 (MQ=255)
                                                                   ttAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCt                        <  1:5008306/93‑1 (MQ=255)
                                                                   ttAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCt                        <  2:3740450/93‑1 (MQ=255)
                                                                    tAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCtt                       >  1:3703456/1‑93 (MQ=255)
                                                                     aGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCttt                      <  1:5907480/93‑1 (MQ=255)
                                                                      ggTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTc                     <  2:7178849/93‑1 (MQ=255)
                                                                      ggTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTc                     <  1:4505945/93‑1 (MQ=255)
                                                                       gcgAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCa                    <  1:6200549/91‑1 (MQ=255)
                                                                          aaCCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCagag                 <  2:3284595/93‑1 (MQ=255)
                                                                          aaCCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCagag                 <  2:3997546/93‑1 (MQ=255)
                                                                           aCCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGt                >  2:7059962/1‑93 (MQ=255)
                                                                              aTGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCAt             <  2:6268145/93‑1 (MQ=255)
                                                                               tGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATg            <  2:4567670/93‑1 (MQ=255)
                                                                                  gagGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGgaa         >  2:483541/1‑93 (MQ=255)
                                                                                  gagGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGgaa         <  1:2797518/93‑1 (MQ=255)
                                                                                    ggCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGgaaga       <  2:3395108/93‑1 (MQ=255)
                                                                                       ggTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGGAAGAATa    >  2:5048741/1‑93 (MQ=255)
                                                                                          tCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGGAAGAATaaa  <  1:4765293/93‑1 (MQ=255)
                                                                                          |                                                                                            
TAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGGAAGAATAAA  >  NC_000913/3899351‑3899532

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