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,4400TC99.1% 750.0 / ‑2.6 226T137A (ACC→GCC) yieCputative receptor protein
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base C (99/125);  total (101/125)
Fisher's exact test for biased strand distribution p-value = 1.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.52e-02

Only 100 of 239 total aligned reads displayed.

TTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGGAAGAAT  >  NC_000913/3899350‑3899529
                                                                                          |                                                                                         
  aaCTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCcccc                                                                                       <  1:854209/93‑1 (MQ=255)
  aaCTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCcccc                                                                                       <  1:3889034/93‑1 (MQ=255)
     tGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCaa                                                                                    >  1:549937/1‑93 (MQ=255)
      gCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAAc                                                                                   >  2:4074998/1‑93 (MQ=255)
       cAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAAcc                                                                                  <  2:4191882/93‑1 (MQ=255)
        aaaTCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCa                                                                                 <  1:3257016/93‑1 (MQ=255)
        aaaTCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCa                                                                                 >  1:2530241/1‑93 (MQ=255)
         aaTCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAg                                                                                <  1:1186889/93‑1 (MQ=255)
         aaTCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAg                                                                                >  1:4171820/1‑93 (MQ=255)
         aaTCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAg                                                                                <  2:4967812/93‑1 (MQ=255)
         aaTCAAACCAGCCGGAGGATTCGTTACCCAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGCGAACCATGAGAGGCGGCCCCCCAACCaa                                                                                >  2:336171/1‑92 (MQ=255)
          aTCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGa                                                                               <  1:480108/93‑1 (MQ=255)
           acAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGaa                                                                              <  1:5786595/92‑1 (MQ=255)
            cAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAAc                                                                             <  2:5173351/93‑1 (MQ=255)
            cAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAAc                                                                             <  2:5331230/93‑1 (MQ=255)
              aaCCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACgg                                                                           >  1:4327680/1‑93 (MQ=255)
              aaCCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACgg                                                                           >  2:4036644/1‑93 (MQ=255)
               aCCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGt                                                                          >  2:3096427/1‑93 (MQ=255)
                ccAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTa                                                                         >  2:5388004/1‑93 (MQ=255)
                  aGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAAt                                                                       >  2:1246324/1‑93 (MQ=255)
                   gCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATa                                                                      <  2:1030342/93‑1 (MQ=255)
                   gCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATa                                                                      <  2:2350042/93‑1 (MQ=255)
                    ccGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACAATGAGAGGCGGCCCCCCAACCAGAACGGTAATAg                                                                     <  1:5823786/93‑1 (MQ=255)
                     cGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCCGCACGGTAATAGc                                                                    >  1:4041729/1‑93 (MQ=255)
                     cGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGc                                                                    <  1:5942978/93‑1 (MQ=255)
                      ggAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGcc                                                                   <  1:2382962/93‑1 (MQ=255)
                      ggAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGcc                                                                   >  2:1483581/1‑93 (MQ=255)
                      ggAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGcc                                                                   <  2:5059105/93‑1 (MQ=255)
                        aGTGTTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGt                                                                 <  1:3953024/93‑1 (MQ=255)
                        aGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGt                                                                 <  1:2902467/93‑1 (MQ=255)
                        aGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGt                                                                 >  1:4364617/1‑93 (MQ=255)
                        aGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTAAGGGGAACAATGGGAGGCGGCCCCCAACCCAGAACGGAAAAAGCCGt                                                                 >  1:1777837/1‑93 (MQ=255)
                         gTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGtt                                                                <  1:2748100/93‑1 (MQ=255)
                         gTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGtt                                                                <  2:3898548/93‑1 (MQ=255)
                         gTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGtt                                                                <  1:1915617/93‑1 (MQ=255)
                          tATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTg                                                               >  2:763876/1‑93 (MQ=255)
                           aTTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGt                                                              >  1:6087682/1‑93 (MQ=255)
                               gTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACt                                                          >  1:2647247/1‑93 (MQ=255)
                               gTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACt                                                          <  1:197229/93‑1 (MQ=255)
                                ttACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTa                                                         >  1:2983431/1‑93 (MQ=255)
                                 tACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTaa                                                        <  1:4967431/93‑1 (MQ=255)
                                  aCCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTaaa                                                       <  2:5752973/93‑1 (MQ=255)
                                  aCCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTaaa                                                       >  2:4502392/1‑93 (MQ=255)
                                  aCCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTaaa                                                       >  2:1694907/1‑93 (MQ=255)
                                   ccAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAg                                                      >  1:3373915/1‑93 (MQ=255)
                                       aGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCaa                                                  <  1:1291101/93‑1 (MQ=255)
                                       aGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCaa                                                  >  2:4346631/1‑93 (MQ=255)
                                        gCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAAt                                                 <  2:5091339/93‑1 (MQ=255)
                                        gCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAAt                                                 >  2:3943394/1‑93 (MQ=255)
                                         cGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAAtt                                                <  2:4572216/93‑1 (MQ=255)
                                         cGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAAtt                                                <  1:635396/93‑1 (MQ=255)
                                          ggCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAAttt                                               <  1:2484595/93‑1 (MQ=255)
                                           gCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTc                                              >  2:4778349/1‑93 (MQ=255)
                                           gCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAAATAAAGTCAATTTc                                              >  1:5216393/1‑93 (MQ=255)
                                             ccAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCAt                                            >  2:2253671/1‑93 (MQ=255)
                                               agagAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCt                                          >  2:2811417/1‑93 (MQ=255)
                                               agagAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCt                                          <  2:239320/93‑1 (MQ=255)
                                               agagAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCt                                          >  1:6262743/1‑93 (MQ=255)
                                               agagAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCt                                          <  2:901187/93‑1 (MQ=255)
                                                gagaACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCtt                                         >  2:356504/1‑93 (MQ=255)
                                                 agaACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGACGTTGTAACTAAAGCCAATTTTATCttt                                        >  2:3303534/1‑93 (MQ=255)
                                                  gaACGAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCGCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTc                                       <  1:3157525/93‑1 (MQ=255)
                                                  gaACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTc                                       <  1:119950/93‑1 (MQ=255)
                                                  gaACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTc                                       <  1:3511146/93‑1 (MQ=255)
                                                   aaCCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTca                                      >  1:1589090/1‑93 (MQ=255)
                                                      cAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACaa                                   >  1:2623957/1‑93 (MQ=255)
                                                       aaTCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACcaa                                  >  2:1102631/1‑93 (MQ=255)
                                                        aTCGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAAt                                 <  2:1926326/93‑1 (MQ=255)
                                                          cGCCCATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAAttt                               <  1:5664711/93‑1 (MQ=255)
                                                            cccATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTAc                             <  1:4562855/93‑1 (MQ=255)
                                                            cccATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTAc                             <  1:2280665/93‑1 (MQ=255)
                                                             ccATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACt                            <  1:1749944/93‑1 (MQ=255)
                                                              cATGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTg                           <  1:3749010/93‑1 (MQ=255)
                                                               aTGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGa                          <  1:5696537/93‑1 (MQ=255)
                                                               aTGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGa                          <  2:688603/93‑1 (MQ=255)
                                                               aTGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGa                          >  1:4722634/1‑93 (MQ=255)
                                                                tGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCCCAAATTTATTGaa                         >  1:2087544/1‑93 (MQ=255)
                                                                tGATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGaa                         <  2:5526701/93‑1 (MQ=255)
                                                                 gATTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGaaa                        <  2:6383537/93‑1 (MQ=255)
                                                                  aTTTAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAAt                       >  1:5025166/1‑93 (MQ=255)
                                                                    ttAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCt                     <  2:5725228/93‑1 (MQ=255)
                                                                     tAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCtt                    <  1:722376/93‑1 (MQ=255)
                                                                     tAGGTGAACCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCtt                    >  1:2189236/1‑93 (MQ=255)
                                                                           aaCCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCagag              <  1:4107667/93‑1 (MQ=255)
                                                                           aaCCATGAGAGGCGCCCCCCCACCCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCatag              <  1:4648171/93‑1 (MQ=255)
                                                                            aCCATGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGt             <  1:3023315/93‑1 (MQ=255)
                                                                             ccATGAGAGGCGGCCCCCCAACCAGAACGGTAATGGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTc            <  1:224910/93‑1 (MQ=255)
                                                                             ccATGAGAGGCGGCCCCCCAACCAGAACGGCAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTc            >  2:3471340/1‑93 (MQ=255)
                                                                               aTGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCAt          <  2:1540077/93‑1 (MQ=255)
                                                                               aTGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCAt          >  1:5102445/1‑93 (MQ=255)
                                                                               aTGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGCTGTAACTAAAGCCACTTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCAt          >  2:2745982/1‑93 (MQ=255)
                                                                                tGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATg         >  2:742586/1‑93 (MQ=255)
                                                                                tGAGAGGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATg         <  2:6259860/93‑1 (MQ=255)
                                                                                 gagaGGCGCCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTTATgg        <  2:6160677/93‑1 (MQ=255)
                                                                                  agagGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGga       >  1:5991945/1‑93 (MQ=255)
                                                                                  agagGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGga       >  1:52983/1‑93 (MQ=255)
                                                                                  agagGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGga       <  1:5591571/93‑1 (MQ=255)
                                                                                   gagGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGgaa      <  1:2405240/93‑1 (MQ=255)
                                                                                   gagGCGGCCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGgaa      <  2:3601477/93‑1 (MQ=255)
                                                                                     ggcggcCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGgaaga    >  2:2222161/1‑93 (MQ=255)
                                                                                          |                                                                                         
TTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGGAAGAAT  >  NC_000913/3899350‑3899529

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