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,7970AC98.9% 974.5 / ‑0.4 271I121L (ATT→CTT) yjgLSopA‑central‑domain‑like hexapeptide repeat protein
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base C (139/129);  total (140/131)
Fisher's exact test for biased strand distribution p-value = 6.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.50e-01

Only 100 of 273 total aligned reads displayed.

TCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATGAAGGGA  >  NC_000913/4475701‑4475893
                                                                                                |                                                                                                
tCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACttt                                                                                               >  1:1358081/1‑100 (MQ=255)
 cTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTc                                                                                              <  1:323248/100‑1 (MQ=255)
 cTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTc                                                                                              >  2:460477/1‑100 (MQ=255)
  ttATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCa                                                                                             <  1:299878/100‑1 (MQ=255)
  ttATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCa                                                                                             >  2:2845455/1‑100 (MQ=255)
   tATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCAc                                                                                            <  1:861777/100‑1 (MQ=255)
   tATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCAc                                                                                            <  1:2890519/100‑1 (MQ=255)
   tATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCAc                                                                                            <  1:2298472/100‑1 (MQ=255)
    aTTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACt                                                                                           >  1:1244558/1‑100 (MQ=255)
      tCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACttt                                                                                         >  1:508881/1‑100 (MQ=255)
       cTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCAACCAAGGTGATTATAGAactttcactttc                                                                                        >  1:3903897/1‑100 (MQ=255)
        tAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTcc                                                                                       <  1:495946/100‑1 (MQ=255)
         aaaGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCa                                                                                      >  1:2648360/1‑100 (MQ=255)
            gATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTa                                                                                   <  1:4090071/100‑1 (MQ=255)
              tACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATc                                                                                 <  2:3881224/100‑1 (MQ=255)
               aCTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCa                                                                                >  1:205127/1‑100 (MQ=255)
               aCTAAAGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAACGTGATTATAGAACTTTCACTTTCCACTATCa                                                                                >  2:4843820/1‑100 (MQ=255)
                cTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCac                                                                               >  2:208442/1‑100 (MQ=255)
                cTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCac                                                                               <  1:2983762/100‑1 (MQ=255)
                 tAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCact                                                                              <  2:1340171/100‑1 (MQ=255)
                 tAAAGGCGAGGACAAAATGAAGTTTTCCATGAATAATGCCGGAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACGTTCCTCTCTCact                                                                              <  2:86188/100‑1 (MQ=255)
                   aTGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCactac                                                                            >  2:3395446/1‑100 (MQ=255)
                   aTGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCactac                                                                            >  2:345726/1‑100 (MQ=255)
                    tGGCAATGAAAAAATGAAGCTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCactact                                                                           <  2:4917880/100‑1 (MQ=255)
                     ggCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCActacta                                                                          <  2:1168232/100‑1 (MQ=255)
                       cAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATg                                                                        <  1:210469/100‑1 (MQ=255)
                       cAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATg                                                                        >  2:547250/1‑100 (MQ=255)
                        aaTGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATgg                                                                       <  1:597756/100‑1 (MQ=255)
                         aTGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATggg                                                                      <  1:3330707/100‑1 (MQ=255)
                         aTGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATggg                                                                      >  2:598186/1‑100 (MQ=255)
                         aTGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATggg                                                                      <  1:2762251/100‑1 (MQ=255)
                         aTGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATggg                                                                      <  1:4516787/100‑1 (MQ=255)
                         aTGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTAACCCCAGGTGATTATAGAACTTTCACTTTCCCCTATCACTACTATggg                                                                      >  2:3305605/1‑100 (MQ=255)
                         aTGAAAAAATGAAGTTGTACATGAATAAGCCTGTAGGGGAACGAGCTGTTTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATggg                                                                      <  1:4188868/100‑1 (MQ=255)
                          tGAAAAAATGAAGTGTTACATGAATAATCCTGTAGGGGAACGAGCTGATGCACCCAAGGTGATTATAGAACTTTCACTTGCCACGATCACTACTATgggg                                                                     >  2:4931521/1‑100 (MQ=255)
                           gAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGa                                                                    >  2:415289/1‑100 (MQ=255)
                           gAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGa                                                                    >  2:5067374/1‑100 (MQ=255)
                             aaaaaTGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACt                                                                  >  1:3155534/1‑100 (MQ=255)
                                aaTGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGt                                                               >  2:2963099/1‑100 (MQ=255)
                                  tGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCa                                                             <  2:5284019/100‑1 (MQ=255)
                                   gAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCaa                                                            <  2:305298/100‑1 (MQ=255)
                                    aaGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAg                                                           <  2:3134350/100‑1 (MQ=255)
                                    aaGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAg                                                           <  1:346091/100‑1 (MQ=255)
                                        tttACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACa                                                       <  2:1472247/100‑1 (MQ=255)
                                         ttACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACAt                                                      >  2:722892/1‑100 (MQ=255)
                                          tACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATa                                                     >  2:4457418/1‑100 (MQ=255)
                                           aCATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATCATAGAACTTTCACTTTCCTCTATCACTACTATGGGGACTCGTCAAGGACATAc                                                    <  1:5164502/100‑1 (MQ=255)
                                            cATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACa                                                   >  2:3169722/1‑100 (MQ=255)
                                              tGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGc                                                 >  2:1541051/1‑100 (MQ=255)
                                               gAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGcc                                                >  1:3966845/1‑100 (MQ=255)
                                                 ataatCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCAt                                              >  2:2334883/1‑100 (MQ=255)
                                                  taatCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCAtt                                             <  1:1646894/100‑1 (MQ=255)
                                                    atCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATtat                                           >  2:3230276/1‑100 (MQ=255)
                                                    atCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATtat                                           <  2:1773212/100‑1 (MQ=255)
                                                    atCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATtat                                           <  2:1043/100‑1 (MQ=255)
                                                       cTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATAtt                                        >  1:1089430/1‑100 (MQ=255)
                                                        tGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATAttt                                       >  2:3457838/1‑100 (MQ=255)
                                                        tGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATAttt                                       >  1:4118407/1‑100 (MQ=255)
                                                        tGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATAttt                                       >  2:560012/1‑100 (MQ=255)
                                                         gTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTc                                      >  2:2553435/1‑100 (MQ=255)
                                                         gTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTc                                      >  1:4264907/1‑100 (MQ=255)
                                                         gTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTc                                      <  1:469300/100‑1 (MQ=255)
                                                         gTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTc                                      <  1:5456087/100‑1 (MQ=255)
                                                          tAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTcc                                     >  1:2371075/1‑100 (MQ=255)
                                                          tAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTcc                                     <  1:559755/100‑1 (MQ=255)
                                                           aGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCca                                    >  1:1822702/1‑100 (MQ=255)
                                                             gggAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCcaca                                  <  1:3621079/100‑1 (MQ=255)
                                                               gAACGAGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAAc                                <  2:3590202/100‑1 (MQ=255)
                                                                aaCGAGCGGATCCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAAcc                               <  1:1186118/100‑1 (MQ=255)
                                                                    aGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACTTGAt                           <  1:3022397/100‑1 (MQ=255)
                                                                    aGCGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGAt                           >  2:3949200/1‑100 (MQ=255)
                                                                      cGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATgg                         <  1:3571512/100‑1 (MQ=255)
                                                                      cGGATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATgg                         <  2:1075059/100‑1 (MQ=255)
                                                                        gATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGtt                       >  1:1402086/1‑100 (MQ=255)
                                                                        gATTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGtt                       >  1:3797835/1‑100 (MQ=255)
                                                                         aTTCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTc                      <  1:553707/100‑1 (MQ=255)
                                                                          ttCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCg                     >  2:4057935/1‑100 (MQ=255)
                                                                          ttCACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCg                     <  1:776873/100‑1 (MQ=255)
                                                                          ttCACCCAAGGTGATGATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCg                     <  1:1067273/100‑1 (MQ=255)
                                                                            cACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGAc                   >  1:2925339/1‑100 (MQ=255)
                                                                            cACCCAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGAc                   <  1:212008/100‑1 (MQ=255)
                                                                            cACCCAAGGTGATTATAGAACTTTCACTTTCCACAATCCCTACTATGGGGACTCGTCAAGTACATACAGCCATTATATTTCCACAACCTGATGGTTCGAc                   >  2:4688875/1‑100 (MQ=255)
                                                                             aCCCAAGGTGATTATAGAACTTTCACTCTCCACTATCACTACTATGGGGACTCGTCAAGCACATACAGCCATTATATTTCCACAACCTGATGGTTCGACt                  <  1:4837009/100‑1 (MQ=255)
                                                                              cccAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTa                 >  1:993317/1‑100 (MQ=255)
                                                                              cccAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTa                 >  1:5332328/1‑100 (MQ=255)
                                                                              cccAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTa                 >  2:2826059/1‑100 (MQ=255)
                                                                              cccAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTa                 >  1:1668515/1‑100 (MQ=255)
                                                                                cAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAAc               <  1:1621277/100‑1 (MQ=255)
                                                                                cAAGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAAc               <  1:4220615/100‑1 (MQ=255)
                                                                                 aaGGTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAAcc              <  1:2327943/100‑1 (MQ=255)
                                                                                   ggTGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGt            >  2:1903157/1‑100 (MQ=255)
                                                                                     taaTTATAGAACTTTCACTTTCCACTATCACTTCTATGAGGAATCGTCAAGGACATACAGCCATTATATTTCCACATCTTGATGGTTCGACTAACCGTTa          <  2:1713111/98‑1 (MQ=255)
                                                                                     tGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTa          <  2:2156297/100‑1 (MQ=255)
                                                                                     tGATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTa          <  2:788677/100‑1 (MQ=255)
                                                                                      gATTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTAt         >  2:2415474/1‑100 (MQ=255)
                                                                                       aTTATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATg        <  1:149412/100‑1 (MQ=255)
                                                                                        ttATAGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATGa       <  2:4368494/100‑1 (MQ=255)
                                                                                         tataGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATGaa      >  1:2462118/1‑100 (MQ=255)
                                                                                          ataGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATGAAg     >  2:3639696/1‑100 (MQ=255)
                                                                                           taGAACTTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATGAAgg    >  2:2488624/1‑100 (MQ=255)
                                                                                                |                                                                                                
TCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATATTTCCACAACCTGATGGTTCGACTAACCGTTATGAAGGGA  >  NC_000913/4475701‑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: