Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,043,629 A→G I185T (ATT→ACT)  yihG ← putative endonuclease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,043,6290AG99.2% 802.8 / ‑3.5 246I185T (ATT→ACT) yihGputative endonuclease
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (103/141);  total (104/142)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

Only 100 of 259 total aligned reads displayed.

CTGCCGCCTTTGGTGGCAACAAGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAATGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGATGACGTAATAAATAA  >  NC_000913/4043537‑4043718
                                                                                            |                                                                                         
  gccgccTTTGGTGGCAACAAGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGtg                                                                                         <  2:418388/93‑1 (MQ=255)
   ccgccTTTGGTGGCAACAAGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGtgg                                                                                        <  1:323175/93‑1 (MQ=255)
   ccgccTTTGGTGGCAACAAGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGtgg                                                                                        <  2:4544264/93‑1 (MQ=255)
    cgccTTTGGTGGCAACAAGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGtggt                                                                                       <  2:1605029/93‑1 (MQ=255)
    cgccTTTGGTGGCAACAAGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGtggt                                                                                       <  2:844137/93‑1 (MQ=255)
     gccTTTGGTGGCAACAAGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAgtggtg                                                                                      <  1:2512615/93‑1 (MQ=255)
     gccTTTGGTGGCAACAAGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAgtggtg                                                                                      >  1:1476705/1‑93 (MQ=255)
       cTTTGGTGGCAACAAGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAACTTGACAGtggtggt                                                                                    >  2:2727424/1‑93 (MQ=255)
       cTTTGGTGGCAACAAGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGtggtggt                                                                                    <  1:1333722/93‑1 (MQ=255)
       cTTTGGTGGCAACAAGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGtggtggt                                                                                    <  2:190314/93‑1 (MQ=255)
       cTTTGGTGGCAACAAGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGtggtggg                                                                                    >  2:3322053/1‑92 (MQ=255)
        tttGGTGGCAACAAGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAgtggtggtg                                                                                   <  2:3765337/93‑1 (MQ=255)
        tttGGTGGCAACAAGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAgtggtggtg                                                                                   >  2:3379545/1‑93 (MQ=255)
        tttGGTGGCAACAAGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAgtggtggtg                                                                                   >  1:2645937/1‑93 (MQ=255)
        tgtggtggCAACAAGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCGGCAACGAAAGAGACAgtggtggtg                                                                                   >  1:4219194/3‑93 (MQ=255)
         ttGGTGGCAACAAGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGtggtggtgg                                                                                  <  2:2382008/93‑1 (MQ=255)
         ttGGTGGCAACAAGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAAtggtggtgg                                                                                  <  1:4007289/93‑1 (MQ=255)
          tggtggCAACAAGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTggg                                                                                 <  1:3987298/93‑1 (MQ=255)
          tggtggCAACAAGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTggg                                                                                 <  2:3260622/93‑1 (MQ=255)
           ggtggCAACAAGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGa                                                                                <  1:1149511/93‑1 (MQ=255)
            gtggCAACAAGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTgggtt                                                                               <  2:4577584/93‑3 (MQ=255)
             tggCAACAAGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATg                                                                              <  2:2394152/93‑1 (MQ=255)
              ggCAACAAGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGc                                                                             <  2:1463376/93‑1 (MQ=255)
              ggCAACAAGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGc                                                                             <  1:102126/93‑1 (MQ=255)
               ggaacaaGTTTTGAAAAGTGGAGTGGGTTTGCTGAGGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGACGCa                                                                            <  1:1070917/91‑1 (MQ=255)
                caacaaGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAATGGTGGTGGGATGCAg                                                                           >  1:4102667/1‑93 (MQ=255)
                 aacaaGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAgg                                                                          >  2:1316262/1‑93 (MQ=255)
                 aacaaGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAgg                                                                          >  2:3825938/1‑93 (MQ=255)
                 aaaaaGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAgg                                                                          <  2:2057692/90‑1 (MQ=255)
                  acaaGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGc                                                                         >  1:1545256/1‑93 (MQ=255)
                  acaaGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGc                                                                         >  2:1981505/1‑93 (MQ=255)
                  acaaGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGc                                                                         <  2:2703288/93‑1 (MQ=255)
                   caaGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCg                                                                        <  1:18427/93‑1 (MQ=255)
                     aGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTGCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGaa                                                                      <  1:411661/93‑1 (MQ=255)
                      gTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGaaa                                                                     >  1:4310208/1‑93 (MQ=255)
                       ttttGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACTAAATTGCCAGTGGTGGTGGGATGCAGGCGaaat                                                                    >  1:463003/1‑92 (MQ=255)
                       ttttGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAAc                                                                    <  1:2871382/93‑1 (MQ=255)
                           gAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACtttt                                                                >  2:3640352/1‑93 (MQ=255)
                              aaGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTcac                                                             >  2:2979667/1‑93 (MQ=255)
                                gtggagtggGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACaa                                                           >  2:3694991/1‑93 (MQ=255)
                                 tggagtggGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAg                                                          >  1:4681968/1‑93 (MQ=255)
                                  ggagtggGTTTGCTGATGTTTTTCTTGCGTGAAACTGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGa                                                         <  2:1372059/93‑1 (MQ=255)
                                  ggagtggGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGa                                                         >  2:1010161/1‑93 (MQ=255)
                                   gagtggGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGaa                                                        <  1:2896126/93‑1 (MQ=255)
                                   gagtggGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGaa                                                        <  1:1931649/93‑1 (MQ=255)
                                   gagtggGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGaa                                                        <  2:3575022/93‑1 (MQ=255)
                                    agtggGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAAc                                                       >  1:806139/1‑93 (MQ=255)
                                    agtggGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAAc                                                       <  2:790638/93‑1 (MQ=255)
                                     gtggGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACg                                                      <  2:4284955/93‑1 (MQ=255)
                                     gtggGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACg                                                      >  1:1285002/1‑93 (MQ=255)
                                      tggGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGc                                                     <  2:2781200/93‑1 (MQ=255)
                                        ggTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCACCTACATTGCCAGTGGTGGTTGAATGCAGGCGAACCTTTTCACAACAACGCCg                                                   >  1:3320481/1‑93 (MQ=255)
                                        ggTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTTGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCg                                                   >  1:3770002/1‑93 (MQ=255)
                                         gTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGa                                                  <  2:1205827/93‑1 (MQ=255)
                                         gTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGa                                                  <  1:1582461/93‑1 (MQ=255)
                                           ttGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGt                                                >  2:3790836/1‑93 (MQ=255)
                                            tGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGATTg                                               <  2:2878057/93‑1 (MQ=255)
                                            tGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTg                                               >  1:3216220/1‑93 (MQ=255)
                                            tGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTg                                               <  2:580260/93‑1 (MQ=255)
                                             gCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCGGGCGAAACTTTTCACAAGAACGCCGAGTgg                                              >  1:941360/1‑93 (MQ=255)
                                              cTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTTGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGt                                             >  1:2081525/1‑93 (MQ=255)
                                              cTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGt                                             <  1:2431628/93‑1 (MQ=255)
                                               tGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGtt                                            <  2:3558002/93‑1 (MQ=255)
                                                gATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGttt                                           >  2:751956/1‑93 (MQ=255)
                                                 aTGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTc                                          >  1:3128286/1‑93 (MQ=255)
                                                  tGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGACGAGGGGGTTCa                                         >  1:2493992/1‑93 (MQ=255)
                                                     ttttCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACa                                      >  2:656054/1‑93 (MQ=255)
                                                       ttCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATc                                    <  1:4651921/93‑1 (MQ=255)
                                                        tCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCt                                   >  1:355241/1‑93 (MQ=255)
                                                         cTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCtt                                  <  1:3407330/93‑1 (MQ=255)
                                                           tGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTa                                <  1:4481066/93‑1 (MQ=255)
                                                              gTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCa                             <  2:274310/93‑1 (MQ=255)
                                                                gAAACGGGAGCCTTCAACGAAATTGACACAGGTGGTGGGATGCAGGCGAAACTTTCCACCAGAACGCCGAGTGGGTTCAACGTCTTTACCACg                           >  2:1550460/1‑93 (MQ=255)
                                                                 aaaCGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGt                          <  1:4603004/93‑1 (MQ=255)
                                                                 aaaCGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGt                          >  2:1728954/1‑93 (MQ=255)
                                                                   aCGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCg                        >  2:1961584/1‑93 (MQ=255)
                                                                      ggAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTc                     >  1:2638116/1‑93 (MQ=255)
                                                                      ggAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTc                     >  2:402935/1‑93 (MQ=255)
                                                                      ggAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTc                     <  1:1722565/93‑1 (MQ=255)
                                                                       gAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTcc                    <  1:4350662/93‑1 (MQ=255)
                                                                       gAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTcc                    <  2:2848273/93‑1 (MQ=255)
                                                                        aGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCg                   >  1:1635827/1‑93 (MQ=255)
                                                                         gCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCgg                  >  1:3601389/1‑93 (MQ=255)
                                                                          ccTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGa                 <  1:4486132/93‑1 (MQ=255)
                                                                           cTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGAt                <  1:765327/93‑1 (MQ=255)
                                                                           cTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGAt                >  2:222124/1‑93 (MQ=255)
                                                                            ttCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGATg               >  2:2646894/1‑93 (MQ=255)
                                                                            ttCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGATg               >  1:141389/1‑93 (MQ=255)
                                                                            ttCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGACTCGACATCTTTACCAAGTCGCTCCGGCTg               >  2:504833/1‑93 (MQ=255)
                                                                             tCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGATGa              >  2:469948/1‑93 (MQ=255)
                                                                              cAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGATGAc             <  2:2245631/93‑1 (MQ=255)
                                                                               aaCGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGATGACg            <  1:3485338/93‑1 (MQ=255)
                                                                                aCGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGATGACGt           >  2:4107211/1‑93 (MQ=255)
                                                                                aCGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGATGACGt           >  1:2226244/1‑93 (MQ=255)
                                                                                   aaaTTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGATGACGtaat        <  1:2970198/93‑1 (MQ=255)
                                                                                   aaaTTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGATGACGtaat        >  1:319827/1‑93 (MQ=255)
                                                                                      ttGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGATGACGTaataaa     >  1:2308710/1‑93 (MQ=255)
                                                                                        gACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGATGACGTaataaata   <  2:629269/93‑1 (MQ=255)
                                                                                        gACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGATGACGTaataaata   <  1:568/93‑1 (MQ=255)
                                                                                         aCAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGATGACGTAataaataa  <  2:2491158/93‑1 (MQ=255)
                                                                                            |                                                                                         
CTGCCGCCTTTGGTGGCAACAAGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAATGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGATGACGTAATAAATAA  >  NC_000913/4043537‑4043718

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