Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,023,531 T→C H318H (CAT→CAC yigC → putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,023,5310TC99.6% 850.7 / NA 263H318H (CAT→CACyigCputative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (152/110);  minor base A (1/0);  total (153/110)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 266 total aligned reads displayed.

TGGCGACCACACCGGTTACTATAATGAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCATTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTTGTGCCGATTCTGCAAAAACAGTT  >  NC_000913/4023441‑4023620
                                                                                          |                                                                                         
tGGCGACCACACCGGTTACTATAATGAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTc                                                                                         >  1:2503536/1‑93 (MQ=255)
 ggCGACCACACCGGTTACTATAATGAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTcc                                                                                        <  1:3352389/93‑1 (MQ=255)
 ggCGACCACACCGGTTACTATAATGAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTcc                                                                                        >  1:4533186/1‑93 (MQ=255)
   cGACCACACCGGTTACTATAATGAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCAc                                                                                      >  2:3476389/1‑93 (MQ=255)
   cGACCACACCGGTTACTATAATGAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCAc                                                                                      <  2:3714532/93‑1 (MQ=255)
    gACCACACCGGTTACTATAATGAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCAcc                                                                                     <  1:4226350/93‑1 (MQ=255)
         cacCGGTTACTATAATGAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATAc                                                                                >  2:3961260/1‑93 (MQ=255)
         cacCGGTTACTATAATGAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATAc                                                                                >  1:6269169/1‑93 (MQ=255)
          acCGGTTACTATAATGAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATAcc                                                                               >  2:3054533/1‑93 (MQ=255)
           ccGGTTACTATAATGAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCg                                                                              >  1:3239455/1‑93 (MQ=255)
            cGGTTACTATAATGAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCgg                                                                             <  2:3682055/93‑1 (MQ=255)
            cGGTTACTATAATGAAGTAGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCgg                                                                             <  2:5063123/93‑1 (MQ=255)
             ggTTACTATAATGAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCggg                                                                            <  1:3780354/93‑1 (MQ=255)
             ggTTACTATAATGAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCggg                                                                            >  1:5995365/1‑93 (MQ=255)
             ggTTACTATAATGAAGTCGATAGTTGCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCTCCTATACCggg                                                                            <  2:1547723/93‑1 (MQ=255)
              gTTACTATAATGAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGc                                                                           >  2:5779146/1‑93 (MQ=255)
              gTTACTATAATGAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGc                                                                           >  2:405167/1‑93 (MQ=255)
                tACTATAATGAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGt                                                                         >  1:988407/1‑93 (MQ=255)
                  cTATAATGAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTcc                                                                       >  1:5918687/1‑93 (MQ=255)
                   tataATGAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCg                                                                      >  1:2557629/1‑93 (MQ=255)
                   tataATGAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCg                                                                      >  1:682553/1‑93 (MQ=255)
                    ataATGAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGc                                                                     >  1:3722602/1‑93 (MQ=255)
                     taATGAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGcc                                                                    <  2:3820976/93‑1 (MQ=255)
                      aaTGAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCa                                                                   <  1:2744546/93‑1 (MQ=255)
                      aaTGAAGTCCATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCa                                                                   <  2:2129664/93‑1 (MQ=255)
                       aTGAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAg                                                                  >  1:6150527/1‑93 (MQ=255)
                        tGAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGa                                                                 >  1:3150795/1‑93 (MQ=255)
                         gAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGAt                                                                <  2:1457959/93‑1 (MQ=255)
                         gAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGAt                                                                >  2:2953015/1‑93 (MQ=255)
                          aaGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTCCCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATg                                                               <  2:3388145/93‑1 (MQ=255)
                          aaGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATg                                                               >  2:2251973/1‑93 (MQ=255)
                          aaGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATg                                                               >  1:3550777/1‑93 (MQ=255)
                          aaGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATg                                                               >  1:5723399/1‑93 (MQ=255)
                          aaGTCGATAGTTTCCCGGGATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACTACTCCACCTATACCGGTCGTCCGCCAGATg                                                               >  2:2413815/1‑93 (MQ=255)
                           aGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGa                                                              <  2:5154045/93‑1 (MQ=255)
                           aGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGa                                                              >  1:2493641/1‑93 (MQ=255)
                            gTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAg                                                             <  1:4411849/93‑1 (MQ=255)
                            gTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAg                                                             >  2:2191044/1‑93 (MQ=255)
                            gTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAg                                                             >  1:4061720/1‑93 (MQ=255)
                              cGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAggc                                                           >  2:889363/1‑91 (MQ=255)
                              cGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGcc                                                           <  2:5551947/93‑1 (MQ=255)
                                gtAGTTTCCCGGTATTTGCCGTGACGCATGTTCACCAGCGTAAAGATGCAATTTACCACTCCACCAATACCGAACGTCCGCCAGATAAGCCcg                                                         <  2:6305700/92‑1 (MQ=255)
                                aTAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCcg                                                         >  2:5209452/1‑93 (MQ=255)
                                     ttCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTg                                                    >  1:5613117/1‑93 (MQ=255)
                                         cGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTgg                                                >  1:159894/1‑93 (MQ=255)
                                          ggTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTggg                                               >  2:6229263/1‑93 (MQ=255)
                                          ggTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTggg                                               >  1:5323937/1‑93 (MQ=255)
                                          ggTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTggg                                               >  1:4882305/1‑93 (MQ=255)
                                           gTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGgt                                              >  2:4444526/1‑93 (MQ=255)
                                            tATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGgtg                                             >  2:675612/1‑93 (MQ=255)
                                              tttACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTc                                           <  2:6110722/93‑1 (MQ=255)
                                                tACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGc                                         >  2:2785153/1‑93 (MQ=255)
                                                 aCCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCa                                        >  1:6024868/1‑93 (MQ=255)
                                                     tGACGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGa                                    >  1:101613/1‑93 (MQ=255)
                                                       aCGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAAc                                  >  2:1254282/1‑93 (MQ=255)
                                                       aCGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAAc                                  <  2:1925238/93‑1 (MQ=255)
                                                       aCGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAAc                                  >  1:1135016/1‑93 (MQ=255)
                                                       aCGCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAAc                                  <  2:5025446/93‑1 (MQ=255)
                                                         gCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGa                                >  2:4239530/1‑93 (MQ=255)
                                                         gCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGa                                <  1:5040363/93‑1 (MQ=255)
                                                         gCATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGa                                >  2:4860884/1‑93 (MQ=255)
                                                          cATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGaa                               <  2:4226045/93‑1 (MQ=255)
                                                          cATATTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGaa                               >  2:2741061/1‑93 (MQ=255)
                                                           atatTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAg                              >  2:5286677/1‑93 (MQ=255)
                                                           atatTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAg                              >  2:431893/1‑93 (MQ=255)
                                                           atatTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAg                              <  1:1431471/93‑1 (MQ=255)
                                                            tatTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGGGTCGCACTGAACGAAgt                             >  2:3200470/1‑93 (MQ=255)
                                                             atTACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAgtg                            <  2:3231025/93‑1 (MQ=255)
                                                              ttACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAgtgt                           >  1:2941806/1‑93 (MQ=255)
                                                               tACCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGtt                          <  1:2151024/93‑1 (MQ=255)
                                                                aCCCAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGttt                         <  1:5287635/93‑1 (MQ=255)
                                                                  ccAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTtgt                       >  2:473001/1‑93 (MQ=255)
                                                                   cAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTtgtg                      >  1:4293307/1‑93 (MQ=255)
                                                                   cAGCGTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTtgtg                      >  1:277245/1‑93 (MQ=255)
                                                                       gTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTTGTGCCgt                  >  1:5223529/1‑92 (MQ=255)
                                                                       gTGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTTGTGCCGa                  <  1:5198801/93‑1 (MQ=255)
                                                                        tGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTTGTGCCGAt                 <  1:4394365/93‑1 (MQ=255)
                                                                        tGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTTGTGCCGAt                 <  1:1988247/93‑1 (MQ=255)
                                                                        tGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTTGTGCCGAt                 >  1:1639522/1‑93 (MQ=255)
                                                                        tGAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTTGTGCCGAt                 >  1:5526719/1‑93 (MQ=255)
                                                                        tGAAGATGAGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTTGTGCCGAt                 >  1:4885621/1‑93 (MQ=255)
                                                                         gAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTTGTGCCGAtt                >  2:947253/1‑93 (MQ=255)
                                                                         gAAGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTTGTGCCGAtt                <  1:1093102/93‑1 (MQ=255)
                                                                          aaGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTTGTGCCGATTc               >  1:861965/1‑93 (MQ=255)
                                                                           aGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTTGTGCCGATTCt              <  2:1513238/93‑1 (MQ=255)
                                                                           aGATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTTGTGCCGATTCt              <  1:54238/93‑1 (MQ=255)
                                                                            gATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTTGTGCCGATTCTg             >  2:3943681/1‑93 (MQ=255)
                                                                            gATGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTTGTGCCGATTCTg             <  1:3865418/93‑1 (MQ=255)
                                                                             aTGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTTGTGCCGATTCTGc            <  2:1159761/93‑1 (MQ=255)
                                                                             aTGCGATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTTGTGCCGATTCTGc            <  1:2360444/93‑1 (MQ=255)
                                                                                 gATTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTTGTGCCGATTCTGCaaaa        >  1:1935067/1‑93 (MQ=255)
                                                                                  aTTTACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTTGTGGCGATTCTGCaaaaa       >  1:1545348/1‑93 (MQ=255)
                                                                                   tttACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTTGTGCCGATTCTGCAAAAAc      <  2:5641979/93‑1 (MQ=255)
                                                                                   tttACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTTGTGCCGATTCTGCAAAAAc      <  2:6074222/93‑1 (MQ=255)
                                                                                    ttACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTTGTGCCGATTCTGCAAAAACa     >  2:3416712/1‑93 (MQ=255)
                                                                                     tACCACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTTGTGCCGATTCTGCAAAAACAg    <  2:3548154/93‑1 (MQ=255)
                                                                                      aCCACTCCACTTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTTGTGCCGATTCTGCAAAAACAGt   >  2:3858682/1‑93 (MQ=255)
                                                                                      aCCACTCCACCTATACCGGTCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTTGTGCCGATTCTGCAAATACAGt   >  1:3811771/1‑93 (MQ=255)
                                                                                       ccACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTTGTGCCGATTCTGCAAAAACAGtt  <  1:3454866/93‑1 (MQ=255)
                                                                                       ccACTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTTGTGCCGATTCTGCAAAAACAGtt  <  2:1005220/93‑1 (MQ=255)
                                                                                          |                                                                                         
TGGCGACCACACCGGTTACTATAATGAAGTCGATAGTTTCCCGGTATTTACCGTGACGCATATTACCCAGCGTGAAGATGCGATTTACCATTCCACCTATACCGGGCGTCCGCCAGATGAGCCCGCGGTGCTGGGTGTCGCACTGAACGAAGTGTTTGTGCCGATTCTGCAAAAACAGTT  >  NC_000913/4023441‑4023620

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