Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,800,657 T→C E10G (GAG→GGG)  rfaI ← UDP‑D‑galactose:(glucosyl)lipopolysaccharide‑ alpha‑1,3‑D‑galactosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,800,6570TC99.0% 634.7 / ‑3.2 201E10G (GAG→GGG) rfaIUDP‑D‑galactose:(glucosyl)lipopolysaccharide‑ alpha‑1,3‑D‑galactosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base C (115/84);  total (115/86)
Fisher's exact test for biased strand distribution p-value = 1.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.67e-01

Only 100 of 205 total aligned reads displayed.

AGAAAGTTTTTGTCAGTTCCATAAGCGATGTCCAGACAAAGATTTTCTGTTTCTACTTTATGATCATAATCAATAACTGAATTAAGAAACTCAGTTTCCTGGAAAAACACCTGCTGCATTATTTCCTCCGGGAAATAGCAGAAATAATTACCTTATGTAATCTATCGTAATACTTCGATT  >  NC_000913/3800567‑3800746
                                                                                          |                                                                                         
aGAAAGTTTTTGTCAGTTCCATAAGCGATGTCCAGACAAAGATTTTCTGTTTCTACTTTATGATCATAATCAATAACTGAATTAAGAAACCCa                                                                                         <  2:1851457/93‑1 (MQ=255)
 gAAAGTTTTTGTCAGTTCCATAAGCGATGTCCAGACAAAGATTTTCTGTTTCTACTTTATGATCATAATCAATAACTGAATTAAGAAACCCAg                                                                                        <  1:1014286/93‑1 (MQ=255)
 gAAAGTTTTTGTCAGTTCCATAAGCGATGTCCAGACAAAGATTTTCTGTTTCTACTTTATGATCATAATCAATAACTGAATTAAGAAACCCAg                                                                                        <  2:4400405/93‑1 (MQ=255)
   aaGTTTTTGTCAGTTCCATAAGCGATGTCCAGACAAAGATTTTCTGTTTCTACTTTATGATCATAATCAATAACTGAATTAAGAAACCCAGtt                                                                                      <  2:1818958/93‑1 (MQ=255)
   aaGTTTTTGTCAGTTCCATAAGCGATGTCCAGACAAAGATTTTCTGTTTCTACTTTATGATCATAATCAATAACTGAATTAAGAAACCCAGtt                                                                                      <  2:2391616/93‑1 (MQ=255)
     gTTTTTGTCAGTTCCATAAGCGATGTCCAGACAAAGATTTTCTGTTTCTACTTTATGATCATAATCAATAACTGAATTAAGAAACCCAGTTTc                                                                                    >  2:1383873/1‑93 (MQ=255)
      tttttGTCAGTTCCATAAGCGATGTCCAGACAAAGATTTTCTGTTTCTACTTTATGATCATAATCAATAACTGAATTAAGAAACCCAGTTTcc                                                                                   <  2:547780/93‑1 (MQ=255)
        tttGTCAGTTCCATAAGCGATGTCCAGACAAAGATTTTCTGTTTCTACTTTATGATCATAATCAATAACTGAATTAAGAAACCCAGTTTCCTg                                                                                 <  1:4120060/93‑1 (MQ=255)
        tttGTCAGTTCCATAAGCGATGTCCAGACAAAGATTTTCTGTTTCTACTTTATGATCATAATCAATAACTGAATTAAGAAACCCAGTTTCCTg                                                                                 <  2:793857/93‑1 (MQ=255)
        tttGTCAGTTCCATAAGCGATGTCCAGACAAAGATTTTCTGTTTCTACTTTATGATCATAATCAATAACTGAATTAAGAAACCCAGTTTCCTg                                                                                 <  1:1600/93‑1 (MQ=255)
         ttGTCAGTTCCATAAGCGATGTCCAGACAAAGATTTTCTGTTTCTACTTTATGATCATAATCAATAACTGAATTAAGAAACCCAGTTTCCTgg                                                                                >  1:4136063/1‑93 (MQ=255)
          tGTCAGTTCCATAAGCGATGTCCAGACAAAGATTTTCTGTTTCTACTTTATGATCATAATCAATAACTGAATTAAGAAACCCAGTTTCCTGGa                                                                               <  1:3411093/93‑1 (MQ=255)
            tCAGTTCCATAAGCGATGTCCAGACAAAGATTTTCTGTTTCTACTTTATGATCATAATCAATAACTGAATTAAGAAACCCAGTTTCCTGGaaa                                                                             >  2:2132044/1‑93 (MQ=255)
               gTTCCATAAGCGATGTCCAGACAAAGATTTTCTGTTTCTACTTTATGATCATAATCAATAACTGAATTAAGAAACCCAGTTTCCTGGaaaaaa                                                                          >  2:2469239/1‑92 (MQ=255)
               gTTCCATAAGCGATGTCCAGACAAAGATTTTCTGTTTCTACTTTATGATCATAATCAATAACTGAATTAAGAAACCCAGTTTCCTGGAAAAac                                                                          >  1:2541374/1‑93 (MQ=255)
                  ccATAAGCGATGTCCAGACAAAGATTTTCTGTTTCTACTTTATGATCATAATCAATAACTGAATTAAGAAACCCAGTTTCCTGGAAAAACAcc                                                                       <  1:2646611/93‑1 (MQ=255)
                  ccATAAGCGATGTCCAGACAAAGATTTTCTGTTTCTACTTTATGATCATAATCAATAACTGAATTAAGAAACCCAGTTTCCTGGAAAAACAcc                                                                       >  1:3459143/1‑93 (MQ=255)
                   cATAAGCGATGTCCAGACAAAGATTTTCTGTTTCTACTTTATGATCATAATCAATAACTGAATTAAGAAACCCAGTTTCCTGGAAAAACACct                                                                      >  1:3739508/1‑93 (MQ=255)
                   cATAAGCGATGTCCAGACAAAGATTTTCTGTTTCTACTTTATGATCATAATCAATAACTGAATTAAGAAACCCAGTTTCCTGGAAAAACACct                                                                      >  1:2982885/1‑93 (MQ=255)
                   cATAAGCGATGTCCAGACAAAGATTTTCTGTTTCTACTTTATGATCATAATCAATAACTGAATTAAGAAACCCAGTTTCCTGGAAAAACACct                                                                      <  2:3335876/93‑1 (MQ=255)
                    aTAAGCGATGTCCAGACAAAGATTTTCTGTTTCTACTTTATGATCATAATCAATAACTGAATTAAGAAACCCAGTTTCCTGGAAAAACACctg                                                                     >  1:1416579/1‑93 (MQ=255)
                     tAAGCGATGTCCAGACAAAGATTTTCTGTTTCTACTTTATGATCATAATCAATAACTGAATTAAGAAACCCAGTTTCCTGGAAAAACACctgc                                                                    <  1:2133895/93‑1 (MQ=255)
                     tAAGCGATGTCCAGACAAAGATTTTCTGTTTCTACTTTATGATCATAATCAATAACTGAATTAAGAAACCCAGTTTCCTGGAAAAACACctgc                                                                    <  2:1486508/93‑1 (MQ=255)
                       aGCGATGTCCAGACAAAGATTTTCTGTTTCTACTTTATGATCATAATCAATAACTGAATTAAGAAACCCAGTTTCCTGGAAAAACACctgctg                                                                  >  2:2808475/1‑93 (MQ=255)
                        gCGATGTCCAGACAAAGATTTTCTGTTTCTACTTTATGATCATAATCAATAACTGAATTAAGAAACCCAGTTTCCTGGAAAAACACCtgctgc                                                                 >  1:3894377/1‑93 (MQ=255)
                           aTGTCCAGACAAAGATTTTCTGTTTCTACTTTATGATCATAATCAATAACTGAATTAAGAAACCCAGTTTCCTGGAAAAACACCTGCTGCatt                                                              >  2:1023375/1‑93 (MQ=255)
                           aTGTCCAGACAAAGATTTTCTGTTTCTACTTTATGATCATAATCAATAACTGAATTAAGAAACCCAGTTTCCTGGAAAAACACCTGCTGCatt                                                              >  2:3141458/1‑93 (MQ=255)
                            tGTCCAGACAAAGATTTTCTGTTTCTACTTTATGATCATAATCAATAACTGAATTAAGAAACCCAGTTTCCTGGAAAAACACCTGCTGCatta                                                             >  2:416681/1‑93 (MQ=255)
                            tGTCCAGACAAAGATTTTCTGTTTCTACTTTATGATCATAATCAATAACTGAATTAAGAAACCCAGTTTCCTGGAAAAACACCTGCTGCatta                                                             <  1:2238606/93‑1 (MQ=255)
                             gTCCAGACAAAGATTTTCTGTTTCTACTTTATGATCATAATCAATAACTGAATTAAGAAACCCAGTTTCCTGGAAAAACACCTGCTGCattat                                                            <  1:4280372/93‑1 (MQ=255)
                              tCCAGACAAAGATTTTCTGTTTCTACTTTATGATCATAATCAATAACTGAATTAAGAAACCCAGTTTCCTGGAAAAACACCTGCTGCattatt                                                           <  1:3029170/93‑1 (MQ=255)
                                cAGACAAAGATTTTCTGTTTCTACTTTATGATCATAATCAATAACTGAATTAAGAAACCCAGTTTCCTGGAAAAACACCTGCTGCATTATTtc                                                         >  1:3588409/1‑93 (MQ=255)
                                 aGACAAAGATTTTCTGTTTCTACTTTATGATCATAATCAATAACTGAATTAAGAAACCCAGTTTCCTGGAAAAACACCTGCTGCATTATTtcc                                                        >  2:1222972/1‑93 (MQ=255)
                                 aGACAAAGATTTTCTGTTTCTACTTTATGATCATAATCAATAACTGAATTAAGAAACCCAGTTTCCTGGAAAAACACCTGCTGCATTATTtcc                                                        <  1:2378856/93‑1 (MQ=255)
                                 aGACAAAGATTTTCTGTTTCTACTTTATGATCATAATCAATAACTGAATTAAGAAACCCAGTTTCCTGGAAAAACACCTGCTGCATTATTtcc                                                        >  2:2825958/1‑93 (MQ=255)
                                  gACAAAGATTTTCTGTTTCTACTTTATGATCATAATCAATAACTGAATTAAGAAACCCAGTTTCCTGGAAAAACACCTGCTGCATTATTtcct                                                       >  1:1230811/1‑93 (MQ=255)
                                  gACAAAGATTTTCTGTTTCTACTTTATGATCATAATCAATAACTGAATTAAGAAACCCAGTTTCCTGGAAAAACACCTGCTGCATTATTtcct                                                       >  1:369694/1‑93 (MQ=255)
                                   aCAAAGATTTTCTGTTTCTACTTTATGATCATAATCAATAACTGAATTAAGAAACCCAGTTTCCTGGAAAAACACCTGCTGCATTATTtcctc                                                      <  2:3465932/93‑1 (MQ=255)
                                   aCAAAGATTTTCTGTTTCTACTTTATGATCATAATCAATAACTGAATTAAGAAACCCAGTTTCCTGGAAAAACACCTGCTGCATTATTtcctc                                                      <  2:2108230/93‑1 (MQ=255)
                                    cAAAGATTTTCTGTTTCTACTTTATGATCATAATCAATAACTGAATTAAGAAACCCAGTTTCCTGGAAAAACACCTGCTGCATTATTtcctcc                                                     <  2:1552588/93‑1 (MQ=255)
                                    cAAAGATTTTCTGTTTCTACTTTATGATCATAATCAATAACTGAATTAAGAAACCCAGTTTCCTGGAAAAACACCTGCTGCATTATTtcctcc                                                     <  1:2770939/93‑1 (MQ=255)
                                     aaaGATTTTCTGTTTCTACTTTATGATCATAATCAATAACTGAATTAAGAAACCCAGTTTCCTGGAAAAACACCTGCTGCATTATTTCCTCCg                                                    >  1:4233217/1‑93 (MQ=255)
                                     aaaGATTTTCTGTTTCTACTTTATGATCATAATCAATAACTGAATTAAGAAACCCAGTTTCCTGGAAAAACACCTGCTGCATTATTTCCTCCg                                                    >  2:28831/1‑93 (MQ=255)
                                       aGATTTTCTGTTTCTACTTTATGATCATAATCAATAACTGAATTAAGAAACCCAGTTTCCTGGAAAAACACCTGCTGCATTATTTCCTCCggg                                                  >  1:3597052/1‑93 (MQ=255)
                                        gATTTTCTGTTTCTACTTTATGATCATAATCAATAACTGAATTAAGAAACCCAGTTTCCTGGAAAAACACCTGCTGCATTATTTCCTCCGGGa                                                 <  1:2034056/93‑1 (MQ=255)
                                        gATTTTCTGTTTCTACTTTATGATCATAATCAATAACTGAATTAAGAAACCCAGTTTCCTGGAAAAACACCTGCTGCATTATTTCCTCCGGGa                                                 <  2:4303504/93‑1 (MQ=255)
                                         aTTTTCTGTTTCTACTTTATGATCATAATCAATAACTGAATTAAGAAACCCAGTTTCCTGGAAAAACACCTGCTGCATTATTTCCTCCGGGaa                                                >  2:4620014/1‑93 (MQ=255)
                                          ttttCTGTTTCTACTTTATGATCATAATCAATAACTGAATTAAGAAACCCAGTTTCCTGGAAAAACACCTGCTGCATTATTTCCTCCGGGaaa                                               >  1:1852471/1‑93 (MQ=255)
                                          ttttCTGTTTCTACTTTATGATCATAATCAATAACTGAATTAAGAAACCCAGTTTCCTGGAAAAACACCTGCTGCATTATTTCCTCCGGGaaa                                               >  1:532692/1‑93 (MQ=255)
                                          ttttCTGTTTCTACTTTATGATCATAATCAATAACTGAATTAAGAAACCCAGTTTCCTGGAAAAACACCTGCTGCATTATTTCCTCCGGGaaa                                               >  2:1537264/1‑93 (MQ=255)
                                           tttCTGTTTCTACTTTATGATCATAATCAATAACTGAATTAAGAAACCCAGTTTCCTGGAAAAACACCTGCTGCATTATTTCCTCCGGGAAAt                                              <  2:3567174/93‑1 (MQ=255)
                                                gTTTCTACTTTATGATCATAATCAATAACTGAATTAAGAAACCCAGTTTCCTGGAAAAACACCTGCTGCATTATTTCCTCCGGGAAATAGCAg                                         >  1:2694899/1‑93 (MQ=255)
                                                 tttCTACTTTATGATCATAATCAATAACTGAATTAAGAAACCCAGTTTCCTGGAAAAACACCTGCTGCATTATTTCCTCCGGGAAATAGCAGa                                        >  2:4052158/1‑93 (MQ=255)
                                                     tACTTTATGATCATAATCAATAACTGAATTAAGAACCCCAGTTTCCTGGAAAAACACCTGCTGCATTATTTCCTCCGGGAAATAGCAGAaata                                    <  1:3223110/93‑1 (MQ=255)
                                                     tACTTTATGATCATAATCAATAACTGAATTAAGAAACCCAGTTTCCTGGAAAAACACCTGCTGCATTATTTCCTCCGGGAAATAGCAGAaata                                    >  1:4597766/1‑93 (MQ=255)
                                                      aCTTTATGATCATAATCAATAACTGAATTAAGAAACCCAGTTTCCTGGAAAAACACCTGCTGCATTATTTCCTCCGGGAAATAGCAGAaataa                                   <  1:2811668/93‑1 (MQ=255)
                                                      aCTTTATGATCATAATCAATAACTGAATTAAGAAACCCAGTTTCCTGGAAAAACACCTGCTGCATTATTTCCTCCGGGAAATAGCAGAaataa                                   <  2:2580037/93‑1 (MQ=255)
                                                      aCTTTATGATCATAATCAATAACTGAATTAAGAAACCCAGTTTCCTGGAAAAACACCTGCTGCATTATTTCCTCCGGGAAATAGCAGAaataa                                   >  2:2710449/1‑93 (MQ=255)
                                                       cTTTATGATCATAATCAATAACTGAATTAAGAAACCCAGTTTCCTGGAAAAGCACCTGCTGCATTATTTCCTCCGGGAAATAGCAGAaataat                                  <  2:2919399/93‑1 (MQ=255)
                                                       cTTTATGATCATAATCAATAACTGAATTAAGAAACCCAGTTTCCTGGAAAAACACCTGCTGCATTATTTCCTCCGGGAAATAGCAGAaataat                                  >  2:3663315/1‑93 (MQ=255)
                                                       cTTTATGATCATAATCAATAACTGAATTAAGAAACCCAGTTTCCTGGAAAAACACCTGCTGCATTATTTCCTCCGGGAAATAGCAGAaataat                                  >  1:2480724/1‑93 (MQ=255)
                                                       cTTTATGATCATAATCAATAACTGAATTAAGAAACCCAGTTTCCTGGAAAAACACCTGCTGCATTATTTCCTCCGGGAAATAGCAGAaataat                                  <  1:1392728/93‑1 (MQ=255)
                                                        tttATGATCATAATCAATAACTGAATTAAGAAACCCAGTTTCCTGGAAAAACACCTGCTGCATTATTTCCTCCGGGAAATAGCAGAAATAAtt                                 >  2:3100225/1‑93 (MQ=255)
                                                        tttATGATCATAATCAATAACTGAATTAAGAAACCCAGTTTCCTGGAAAAACACCTGCTGCATTATTTCCTCCGGGAAATAGCAGAAATAAtt                                 <  2:46099/93‑1 (MQ=255)
                                                         ttATGATCATAATCAATAACTGAATTAAGAAACCCAGTTTCCTGGAAAAACACCTGCTGCATTATTTCCTCCGGGAAATAGCAGAAATAATTa                                <  2:3641187/93‑1 (MQ=255)
                                                         ttATGATCATAATCAATAACTGAATTAAGAAACCCAGTTTCCTGGAAAAACACCTGCTGCATTATTCCCTCCGGGAAATAGCAGAAATAATTa                                <  2:2315200/93‑1 (MQ=255)
                                                          tATGATCATAATCAATAACTGAATTAAGAAACCCAGTTTCCTGGAAAAACACCTGCTGCATTATTTCCTCCGGGAAATAGCAGAAATAATTAc                               >  2:4686394/1‑93 (MQ=255)
                                                           aTGATCATAATCAATAACTGAATTAAGAAACCCAGTTTCCTGGAAAAACACCTGCTGCATTATTTCCTCCGGGAAATAGCAGAAATAATTAcc                              <  1:1382432/93‑1 (MQ=255)
                                                             gATCATAATCAATAACTGAATTAAGAAACCCAGTTTCCTGGAAAAACACCTGCTGCATTATTTCCTCCGGGAAATAGCAGAAATAATTACCtt                            >  1:976871/1‑93 (MQ=255)
                                                             gATCATAATCAATAACTGAATTAAGAAACCCAGTTTCCTGGAAAAACACCTGCTGCATTATTTCCTCCGGGAAATAGCAGAAATAATTACCtt                            >  1:4414460/1‑93 (MQ=255)
                                                              acattAATCAATAACTGAATTAAGAACCCAAGTTTCCTGGAAAACCCCCGCTTGAATTATTCCCTCGGGGAAATAGCAGAAATAATTACCTTa                           <  1:738380/89‑1 (MQ=255)
                                                              aTCATAATCAATAACTGAATTAAGAAACCCAGTTTCCTGGAAAAACACCTGCTGCATTATTTCCTCCGGGAAATAGCAGAAATAATTACCTTa                           <  2:4277996/93‑1 (MQ=255)
                                                                cATAATCAATAACTGAATTAAGAAACCCAGTTTCCTGGAAAAACACCTGCTGCATTATTTCCTCCGGGAAATAGCAGAAATAATTACCTTATg                         <  1:4353360/93‑1 (MQ=255)
                                                                 aTAATCAATAACTGAATTAAGAAACCCAGTTTCCTGGAAAAACACCTGCTGCATTATTTCCTCCGGGAAATAGCAGAAATAATTACCTTATGt                        <  2:1390817/93‑1 (MQ=255)
                                                                  tAATCAATAACTGAATTAAGAAACCCAGTTTCCTGGAAAAACACCTGCTGCATTATTTCCTCCGGGAAATAGCAGAAATAATTACCTTATGTa                       <  1:1888044/93‑1 (MQ=255)
                                                                   aaTCAATAACTGAATTAAGAAACCCAGTTTCCTGGAAAAACACCTGCTGCATTATTTCCTCCGGGAAATAGCAGAAATAATTACCTTATGTaa                      <  2:3233668/93‑1 (MQ=255)
                                                                   aaTCAATAACTGAATTAAGAAACCCAGTTTCCTGGAAAAACACCTGCTGCATTATTTCCTCCGGGAAATAGCAGAAATAATTACCTTATGTaa                      >  2:4474657/1‑93 (MQ=255)
                                                                   aaTCAATAACTGAATTAAGAAACCCAGTTTCCTGGAAAAACACCTGCTGCATTATTTCCTCCGGGAAATAGCAGAAATAATTACCTTATGTaa                      >  1:3925318/1‑93 (MQ=255)
                                                                   aaTCAATAACTGAATTAAGAAACCCAGTTTCCTGGAAAAACACCTGCTGCATTATTTCCTCCGGGAAATAGCAGAAATAATTACCTTATGTaa                      >  2:3985563/1‑93 (MQ=255)
                                                                    aTCAATAACTGAATTAAGAAACCCAGTTTCCTGGAAAAACACCTGCTGCATTATTTCCTCCGGGAAATAGCAGAAATAATTACCTTATGTAAt                     >  1:166663/1‑93 (MQ=255)
                                                                     tCAATAACTGAATTAAGAAACCCAGTTTCCTGGAAAAACACCTGCTGCATTATTTCCTCCGGGAAATAGCAGAAATAATTACCTTATGTAATc                    <  2:3872886/93‑1 (MQ=255)
                                                                        aTAACTGAATTAAGAAACCCAGTTTCCTGGAAAAACACCTGCTGCATTATTTCCTCCGGGAAATAGCAGAAATAATTACCTTATGTAATCTAt                 >  1:3643538/1‑93 (MQ=255)
                                                                          aaCTGAATTAAGAAACCCAGTTTCCTGGAAAAACACCTGCTGCATTATTTCCTCCGGGAAATAGCAGAAATAATTACCTTATGTAATCTATCg               <  2:3533445/93‑1 (MQ=255)
                                                                            cTGAATTAAGAAACCCAGTTTCCTGGAAAAACACCTGCTGCATTATTTCCTCCGGGAAATAGCAGAAATAATTACCTTATGTAATCTATCGTa             <  2:119705/93‑1 (MQ=255)
                                                                            cTGAATTAAGAAACCCAGTTTCCTGGAAAAACACCTGCTGCATTATTTCCTCCGGGAAATAGCAGAAATAATTACCTTATGTAATCTATCGTa             >  1:2190890/1‑93 (MQ=255)
                                                                             tGAATTAAGAAACCCAGTTTCCTGGAAAAACACCTGCTGCATTATTTCCTCCGGGAAATAGCAGAAATAATTACCTTATGTAATCTATCGTaa            <  2:2288909/93‑1 (MQ=255)
                                                                              gAATTAAGAAACCCAGTTTCCTGGAAAAACACCTGCTGCATTATTTCCTCCGGGAAATAGCAGAAATAATTACCTTATGTAATCTATCGTAAt           <  2:1732776/93‑1 (MQ=255)
                                                                               aaTTAAGAAACCCAGTTTCCTGGAAAAACACCTGCTGCATTATTTCCTCCGGGAAATAGCAGAAATAATTACCTTATGTAATCTATCGTAATa          <  1:709471/93‑1 (MQ=255)
                                                                               aaTTAAGAAACCCAGTTTCCTGGAAAAACACCTGCTGCATTATTTCCTCCGGGAAATAGCAGAAATAATTACCTTATGTAATCTATCGTAATa          >  1:813854/1‑93 (MQ=255)
                                                                               aaTTAAGAAACCCAGTTTCCTGGAAAAACACCTGCTGCATTATTTCCTCCGGGAAATAGCAGAAATAATTACCTTATGTAATCTATCGTAATa          >  1:4047596/1‑93 (MQ=255)
                                                                                aTTAAGAAACCCAGTTTCCTGGAAAAACACCTGCTGCATTATTTCCTCCGGGAAATAGCAGAAATAATTACCTTATGTAATCTATCGTAATAc         >  2:3263266/1‑93 (MQ=255)
                                                                                 ttAAGAAACCCAGTTTCCTGGAAAAACACCTGCTGCATTATTTCCTCCGGGAAATAGCAGAAATAATTACCTTATGTAATCTATCGTAATACt        >  1:1533743/1‑93 (MQ=255)
                                                                                 ttAAGAAACCCAGTTTCCTGGAAAAACACCTGCTGCATTATTTCCTCCGGGAAATAGCAGAAATAATTACCTTATGTAATCTATCGTAATACt        >  1:3297174/1‑93 (MQ=255)
                                                                                 ttAAGAAACCCAGTTTCCTGGAAAAACACCTGCTGCATTATTTCCTCCGGGAAATAGCAGAAATAATTACCTTATGTAATCTATCGTAATACt        >  2:2224549/1‑93 (MQ=255)
                                                                                 ttAAGAAACCCAGTTTCCTGGAAAAACACCTGCTGCATTATTTCCTCCGGGAAATAGCAGAAATAATTACCTTATGTAATCTATCGTAATACt        >  2:2797740/1‑93 (MQ=255)
                                                                                  tAAGAAACCCAGTTTCCTGGAAAAACACCTGCTGCATTATTTCCTCCGGGAAATAGCAGAAATAATTACCTTATGTAATCTATCGTAATACtt       <  2:2725470/93‑1 (MQ=255)
                                                                                   aaGAAACCCAGTTTCCTGGAAAAACACCTGCTGCATTATTTCCTCCGGGAAATAGCAGAAATAATTACCTTATGTAATCTATCGTAATACTTc      <  1:4673332/93‑1 (MQ=255)
                                                                                     gAAACCCAGTTTCCTGGAAAAACACCTGCTGCATTATTTCCTCCGGGAAATAGCAGAAATAATTACCTTATGTAATCTATCGTAATACTTCGa    >  1:1707249/1‑93 (MQ=255)
                                                                                       aaCCCAGTTTCCTGGAAAAACACCTGCTGCATTATTTCCTCCGGGAAATAGCAGAAATAATTACCTTATGTAATCTATCGTAATACTTCGAtt  >  2:4456816/1‑93 (MQ=255)
                                                                                       aaCCCAGTTTCCTGGAAAAACACCTGCTGCATTATTTCCTCCGGGAAATAGCAGAAATAATTACCTTATGTAATCTATCGTAATACTTCGAtt  >  2:619323/1‑93 (MQ=255)
                                                                                          |                                                                                         
AGAAAGTTTTTGTCAGTTCCATAAGCGATGTCCAGACAAAGATTTTCTGTTTCTACTTTATGATCATAATCAATAACTGAATTAAGAAACTCAGTTTCCTGGAAAAACACCTGCTGCATTATTTCCTCCGGGAAATAGCAGAAATAATTACCTTATGTAATCTATCGTAATACTTCGATT  >  NC_000913/3800567‑3800746

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