Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,475,797 A→C I121L (ATT→CTT)  yjgL → SopA‑central‑domain‑like hexapeptide repeat protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,475,7970AC96.5% 578.6 / 2.3 175I121L (ATT→CTT) yjgLSopA‑central‑domain‑like hexapeptide repeat protein
Reads supporting (aligned to +/- strand):  ref base A (5/1);  new base C (76/91);  total (83/92)
Fisher's exact test for biased strand distribution p-value = 9.95e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 177 total aligned reads displayed.

CATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATGAAGGGA  >  NC_000913/4475699‑4475893
                                                                                                  |                                                                                                
cATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATGCACTCAAGGTGATTATAGAAAt                                                                                                 >  1:2452088/1‑100 (MQ=255)
   cTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTc                                                                                              <  2:2361716/100‑1 (MQ=255)
   cTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTc                                                                                              <  1:3053984/100‑1 (MQ=255)
    ttATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCa                                                                                             <  1:2018977/100‑1 (MQ=255)
     tATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCAc                                                                                            <  2:1800762/100‑1 (MQ=255)
      aTTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACt                                                                                           <  2:204293/100‑1 (MQ=255)
        tCTAAAGATACTAATGGCAACGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACttt                                                                                         >  1:3182039/1‑100 (MQ=255)
         cTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAactttcactttc                                                                                        >  2:2442428/1‑100 (MQ=255)
         cTAAAGATACTAATAGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTc                                                                                        >  2:491315/1‑100 (MQ=255)
          tAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTcc                                                                                       >  2:2126685/1‑100 (MQ=255)
          tAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTcc                                                                                       <  1:2501669/100‑1 (MQ=255)
            aaGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCAc                                                                                     <  2:755301/100‑1 (MQ=255)
                 aCTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCa                                                                                <  1:2870831/100‑1 (MQ=255)
                  cTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCac                                                                               >  1:4063697/1‑100 (MQ=255)
                  cTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCac                                                                               <  1:539313/100‑1 (MQ=255)
                  cTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCac                                                                               >  1:936540/1‑100 (MQ=255)
                  cTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCac                                                                               >  2:3547206/1‑100 (MQ=255)
                  cTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCac                                                                               >  1:1147815/1‑100 (MQ=255)
                   tAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCact                                                                              <  2:2671213/100‑1 (MQ=255)
                   tAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTAGAGAACTTTCACTTTCCACTATCact                                                                              <  2:3185749/100‑1 (MQ=255)
                      tGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCactact                                                                           >  2:186680/1‑100 (MQ=255)
                        gCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTAt                                                                         >  2:3416920/1‑100 (MQ=255)
                        gCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTAt                                                                         >  2:3308274/1‑100 (MQ=255)
                        gCAATGAAAAAANGAAGTTTTACATCAATAATCCTGTCGGGGAACGGGCGTATTCACCCAAGGTGGTTATAGAACTTTCACTTTCCACTATCACTACTAt                                                                         <  1:3796479/100‑1 (MQ=255)
                         cAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATg                                                                        >  2:496701/1‑100 (MQ=255)
                         cAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATg                                                                        >  2:1717960/1‑100 (MQ=255)
                         cAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATg                                                                        <  2:1294189/100‑1 (MQ=255)
                           aTGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATggg                                                                      >  1:215038/1‑100 (MQ=255)
                            tGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATgggg                                                                     <  1:655157/100‑1 (MQ=255)
                            tGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATgggg                                                                     <  1:362468/100‑1 (MQ=255)
                             gAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGa                                                                    >  2:1088866/1‑100 (MQ=255)
                             gAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGa                                                                    >  2:60741/1‑100 (MQ=255)
                              aaaaaaTGAAGTTTTACATGGATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGAc                                                                   <  1:3156098/100‑1 (MQ=255)
                              aaaaaaTAAAGTTTTACATGAATAATCCTGTAGGTGAACTTGCGGATTCACCCAAGGTGATTATAGAACTTTTACTTTCCACTATCACTACTATGGGGAc                                                                   <  2:3017990/100‑1 (MQ=255)
                                 aaaTGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCg                                                                <  1:1828799/100‑1 (MQ=255)
                                   aTGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTc                                                              >  2:2542719/1‑100 (MQ=255)
                                    tGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCa                                                             <  2:3958598/100‑1 (MQ=255)
                                      aaGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAg                                                           >  1:1023636/1‑100 (MQ=255)
                                       aGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAgg                                                          <  1:454723/100‑1 (MQ=255)
                                        gTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGa                                                         >  2:3367897/1‑100 (MQ=255)
                                         ttctACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGAc                                                        <  2:2565384/97‑1 (MQ=255)
                                           ttACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACAt                                                      <  1:647577/100‑1 (MQ=255)
                                             aCATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATAc                                                    >  2:2479629/1‑100 (MQ=255)
                                              cATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACa                                                   >  1:3976203/1‑100 (MQ=255)
                                              cATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACa                                                   >  1:224437/1‑100 (MQ=255)
                                              cATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACa                                                   >  1:74394/1‑100 (MQ=255)
                                              cATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACa                                                   >  1:32885/1‑100 (MQ=255)
                                               aTGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAg                                                  <  1:1521465/100‑1 (MQ=255)
                                                tGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGc                                                 <  2:1619306/100‑1 (MQ=255)
                                                tGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGc                                                 >  1:247179/1‑100 (MQ=255)
                                                 gAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGcc                                                <  1:473315/100‑1 (MQ=255)
                                                    taatCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCAtt                                             <  1:3480943/100‑1 (MQ=255)
                                                     aatCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATta                                            <  1:61818/100‑1 (MQ=255)
                                                     aatCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATta                                            <  2:424464/100‑1 (MQ=255)
                                                      atCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCGCTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATtat                                           >  1:1628879/1‑100 (MQ=255)
                                                       tCCTGTAGGGGAACGAGCGGATTCACCCAAGGTTATTATAGAACTTTCACTTTCCACTAACACTACAATGGGGACTCGTCAAGGACATACAGCCATtata                                          <  2:1610827/100‑1 (MQ=255)
                                                           gTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTc                                      >  2:1896782/1‑100 (MQ=255)
                                                            tAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTcc                                     >  2:1251608/1‑100 (MQ=255)
                                                              ggggAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCcac                                   >  1:1708190/1‑100 (MQ=255)
                                                               gggAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCcaca                                  <  1:2717918/100‑1 (MQ=255)
                                                                  aaCGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAAcc                               <  1:1560995/100‑1 (MQ=255)
                                                                    cGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTg                             >  1:3282203/1‑100 (MQ=255)
                                                                      aGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGan                           <  1:3260824/100‑2 (MQ=255)
                                                                      aGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGAt                           >  2:1198622/1‑100 (MQ=255)
                                                                      aGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGAt                           <  1:2269573/100‑1 (MQ=255)
                                                                      aGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGAt                           <  2:3097698/100‑1 (MQ=255)
                                                                      aGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGAt                           >  1:2675921/1‑100 (MQ=255)
                                                                        cGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATgg                         <  2:2164834/100‑1 (MQ=255)
                                                                         ggATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGt                        >  2:1065674/1‑100 (MQ=255)
                                                                           aTTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTc                      <  1:3652856/100‑1 (MQ=255)
                                                                           aTTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTc                      >  2:2697955/1‑100 (MQ=255)
                                                                           aTTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTc                      <  1:508963/100‑1 (MQ=255)
                                                                            ttCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCg                     <  1:1363222/100‑1 (MQ=255)
                                                                            ttCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCg                     <  2:1002137/100‑1 (MQ=255)
                                                                              cACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGAc                   >  2:3396414/1‑100 (MQ=255)
                                                                              cACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGAc                   >  1:2947912/1‑100 (MQ=255)
                                                                              cACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGAc                   <  2:2841338/100‑1 (MQ=255)
                                                                              atctcAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGAc                   <  2:1688499/96‑1 (MQ=255)
                                                                                cccAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGTCATTATATTTCCACAACCTGATGGTTCGACTa                 >  1:3740666/1‑100 (MQ=255)
                                                                                cccAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTa                 >  2:2219212/1‑100 (MQ=255)
                                                                                cccAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTa                 >  1:2554243/1‑100 (MQ=255)
                                                                                cccAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTa                 <  1:2377559/100‑1 (MQ=255)
                                                                                  cAAGGTGATTATAGAACTTTCACTTTCCAGTATCACTACTATGGGGACTCGTCAAGGAAATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAAc               >  1:3908090/1‑100 (MQ=255)
                                                                                  cAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAAc               >  2:1522417/1‑100 (MQ=255)
                                                                                       tGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTa          >  1:839829/1‑100 (MQ=255)
                                                                                       tGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTa          <  2:2326487/100‑1 (MQ=255)
                                                                                       tGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTa          <  1:3117480/100‑1 (MQ=255)
                                                                                       tGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTa          <  1:1538405/100‑1 (MQ=255)
                                                                                       tGATTATAGAACTTTCACTTTCCACTATCACTAATATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTa          <  2:3235339/100‑1 (MQ=255)
                                                                                       tGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTa          >  2:379895/1‑100 (MQ=255)
                                                                                        gATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTAt         <  1:815940/100‑1 (MQ=255)
                                                                                        gATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTAt         <  1:3945769/100‑1 (MQ=255)
                                                                                        gATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTAt         <  1:2136459/100‑1 (MQ=255)
                                                                                        gATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTAt         <  1:1904186/100‑1 (MQ=255)
                                                                                          ttATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATGa       <  1:2768753/100‑1 (MQ=255)
                                                                                           tataGAACTTTCACTTTCCACTATCACTACTATGGGGTCTCGTCAAGAACATACAGCCATTATATTTCCCCAACCTAATGGTTCGACTAGCCGTTagtaa      >  2:3360311/1‑96 (MQ=255)
                                                                                             taGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATGAAgg    >  2:1875896/1‑100 (MQ=255)
                                                                                               gAACTTTCACTTTCCACTATCACTATTGTGGGGACTCGTGTAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATGAAGGGa  <  2:3821621/100‑1 (MQ=255)
                                                                                               gAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATGAAGGGa  >  1:4003326/1‑100 (MQ=255)
                                                                                               gAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATGAAGGGa  <  2:916862/100‑1 (MQ=255)
                                                                                                  |                                                                                                
CATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATGAAGGGA  >  NC_000913/4475699‑4475893

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