Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,968,566 A→G L287L (TTA→TTG wecB → UDP‑N‑acetyl glucosamine ‑2‑epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,968,5660AG99.3% 1389.5 / ‑2.5 432L287L (TTA→TTGwecBUDP‑N‑acetyl glucosamine ‑2‑epimerase
Reads supporting (aligned to +/- strand):  ref base A (3/0);  new base G (223/205);  total (226/206)
Fisher's exact test for biased strand distribution p-value = 2.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.84e-01

Only 100 of 446 total aligned reads displayed.

GGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTCGCTGGGGAAACCTGT  >  NC_000913/3968476‑3968655
                                                                                          |                                                                                         
ggTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGcc                                                                                         <  1:4042089/93‑1 (MQ=255)
 gTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCg                                                                                        <  2:2988006/93‑1 (MQ=255)
   gCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGtt                                                                                      <  1:7137256/93‑1 (MQ=255)
    cATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGttt                                                                                     >  2:7362592/1‑93 (MQ=255)
     aTCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTg                                                                                    >  2:3755039/1‑93 (MQ=255)
     aTCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTg                                                                                    >  1:3268807/1‑93 (MQ=255)
       ctcAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGGCGTTTGTc                                                                                  <  1:6308435/93‑1 (MQ=255)
       ctcAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTc                                                                                  >  1:9989949/1‑93 (MQ=255)
       ctcAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTc                                                                                  >  1:4285390/1‑93 (MQ=255)
        tcAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCt                                                                                 <  2:568729/93‑1 (MQ=255)
            cccGAACGTCAGAGAACCGATCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCt                                                                             <  2:7566639/93‑1 (MQ=255)
             ccGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCtg                                                                            >  1:10005844/1‑93 (MQ=255)
              cGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCtga                                                                           >  2:3864161/1‑93 (MQ=255)
              cGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCtga                                                                           <  2:1949652/93‑1 (MQ=255)
              cGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCtga                                                                           <  1:7817739/93‑1 (MQ=255)
               gAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCtgat                                                                          >  2:7705531/1‑93 (MQ=255)
                aaCGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCtgatg                                                                         >  2:6560653/1‑93 (MQ=255)
                aaCGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCtgatg                                                                         <  1:9369656/93‑1 (MQ=255)
                aaCGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCtgatg                                                                         >  2:5056183/1‑93 (MQ=255)
                aaCGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGTTTGTCTGGCtgatg                                                                         >  2:1776039/1‑93 (MQ=255)
                 aCGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCtgatga                                                                        >  2:8361591/1‑93 (MQ=255)
                  ggTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGaa                                                                       <  2:8678093/92‑1 (MQ=255)
                    tCAGAGAACCGGTCAATCGTATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATGTGCCGTTTTTCTGGCTGATGAAcc                                                                     >  1:8963797/1‑93 (MQ=255)
                    tCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAAcc                                                                     <  2:8031106/93‑1 (MQ=255)
                    tCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAAcc                                                                     >  1:9448180/1‑93 (MQ=255)
                     cAGAGAACCGGTCAATCGCATTCTTGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCa                                                                    <  1:5828040/93‑1 (MQ=255)
                      agagAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCACAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCAc                                                                   >  1:1523131/1‑93 (MQ=255)
                        agaACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGc                                                                 <  1:676397/93‑1 (MQ=255)
                        agaACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGc                                                                 <  2:60038/93‑1 (MQ=255)
                        agaACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGc                                                                 <  1:3385199/93‑1 (MQ=255)
                          aaCCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCt                                                               <  2:5790528/93‑1 (MQ=255)
                          aaCCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCt                                                               <  1:7304534/93‑1 (MQ=255)
                          aaCCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCt                                                               >  1:6988564/1‑93 (MQ=255)
                            ccGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTgg                                                             <  2:3102220/93‑1 (MQ=255)
                             cGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGc                                                            >  1:993655/1‑93 (MQ=255)
                              ggTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCt                                                           >  2:7866167/1‑93 (MQ=255)
                                tCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGa                                                         >  2:6092857/1‑93 (MQ=255)
                                tCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGa                                                         >  1:4424928/1‑93 (MQ=255)
                                tCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGa                                                         >  2:1550490/1‑93 (MQ=255)
                                 cAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGAt                                                        >  2:2210831/1‑93 (MQ=255)
                                 cAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGAt                                                        >  1:2277463/1‑93 (MQ=255)
                                 cAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGAt                                                        <  2:2459789/93‑1 (MQ=255)
                                 cAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGAt                                                        >  1:1622301/1‑93 (MQ=255)
                                  aaTCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGAtt                                                       <  1:5503031/93‑1 (MQ=255)
                                  aaTCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGAtt                                                       <  1:8436009/93‑1 (MQ=255)
                                    tCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGAtttt                                                     >  1:6802324/1‑93 (MQ=255)
                                     cGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTg                                                    <  1:749305/93‑1 (MQ=255)
                                     cGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGGATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTg                                                    >  2:6755288/1‑93 (MQ=255)
                                       cATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGAc                                                  >  1:5962956/1‑93 (MQ=255)
                                         ttCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCg                                                >  1:5085079/1‑93 (MQ=255)
                                         ttCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCg                                                >  2:4178757/1‑93 (MQ=255)
                                         ttCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCg                                                <  2:7615938/93‑1 (MQ=255)
                                          tCTGGGGGATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGTTGATTTTGACCGa                                               <  2:979559/93‑1 (MQ=255)
                                          tCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGa                                               >  1:1070387/1‑93 (MQ=255)
                                            tGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACt                                             >  1:1405856/1‑93 (MQ=255)
                                               ggCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAg                                          >  1:365832/1‑93 (MQ=255)
                                               ggCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAg                                          <  2:9259392/93‑1 (MQ=255)
                                                gCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAgg                                         >  1:6042756/1‑93 (MQ=255)
                                                 cATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAggc                                        >  1:6453610/1‑93 (MQ=255)
                                                    gtgAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAggcggc                                     <  2:1410006/93‑1 (MQ=255)
                                                     tgAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCa                                    <  1:1923850/93‑1 (MQ=255)
                                                        aaaaTGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTc                                 <  1:4884513/93‑1 (MQ=255)
                                                          aaTGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGGTTTTGACCGACTCAGGAGGCATTCAg                               >  2:6362186/1‑93 (MQ=255)
                                                          aaTGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAg                               >  1:5290357/1‑93 (MQ=255)
                                                           aTGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAgg                              >  1:123728/1‑93 (MQ=255)
                                                           aTGTCATCCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAgg                              <  2:913024/93‑1 (MQ=255)
                                                            tGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGga                             >  1:218341/1‑93 (MQ=255)
                                                             gTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGgaa                            >  2:9556799/1‑93 (MQ=255)
                                                              tCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGgaag                           >  1:3018754/1‑93 (MQ=255)
                                                                aTTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCGCGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGgaagaa                         <  2:6934511/93‑1 (MQ=255)
                                                                aTTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGgaagaa                         <  2:4378459/93‑1 (MQ=255)
                                                                aTTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGgaagaa                         <  1:8697086/93‑1 (MQ=255)
                                                                aTTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGgaagaa                         >  1:3565189/1‑93 (MQ=255)
                                                                 ttCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGGATTCAGGaagaag                        >  1:9542350/1‑93 (MQ=255)
                                                                 ttCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGaagaag                        <  1:5639764/93‑1 (MQ=255)
                                                                  tCTGATCGATCCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAgc                       <  2:2686299/93‑1 (MQ=255)
                                                                     gatcgatcCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGcc                    <  1:4780551/93‑1 (MQ=255)
                                                                       tcgatcCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCtt                  >  1:3882498/1‑93 (MQ=255)
                                                                        cgatcCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTc                 >  1:8008703/1‑93 (MQ=255)
                                                                           tcCCCAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTCGCt              >  2:4695752/1‑93 (MQ=255)
                                                                            tcccAGGAGTATTTGCCGTTTGTCTGGGTGATGAACCACGCCGGGCTGATTTTGACCGGCTCAGGCGGCATTCAGGAAGAAGCGCCTTCGCTg             <  1:1728829/92‑1 (MQ=255)
                                                                             cccAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGTCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTCGCTgg            <  2:4797298/93‑1 (MQ=255)
                                                                             cccAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTCGCTgg            <  1:2513478/93‑1 (MQ=255)
                                                                             cccAGGAGTATTGGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTCGCTgg            <  2:1678966/93‑1 (MQ=255)
                                                                              ccAGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTCGCTggg           >  2:8953488/1‑93 (MQ=255)
                                                                                aGGAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCTGGAAGAAGCGCCTTCGCTGGGGa         >  2:5522056/1‑93 (MQ=255)
                                                                                 ggAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTCGCTGGGGaa        >  2:5255082/1‑93 (MQ=255)
                                                                                  gAGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTCGCTGGGGaaa       >  1:8180277/1‑93 (MQ=255)
                                                                                   aGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTCGCTGGGGAAAc      <  1:1107753/93‑1 (MQ=255)
                                                                                   aGTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTCGCTGGGGAAAc      <  2:3562321/93‑1 (MQ=255)
                                                                                    gTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTCGCTGGGGAAAcc     >  2:7196263/1‑93 (MQ=255)
                                                                                    gTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTCGCTGGGGAAAcc     <  1:2441348/93‑1 (MQ=255)
                                                                                    gTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTCGCTGGGGAAAcc     <  1:4182680/93‑1 (MQ=255)
                                                                                    gTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTCGCTGGGGAAAcc     <  2:8186183/93‑1 (MQ=255)
                                                                                    gTATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTCGCTGGGGAAAcc     >  1:4805480/1‑93 (MQ=255)
                                                                                     tATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTCGCTGGGGAAACCt    >  1:2765120/1‑93 (MQ=255)
                                                                                     tATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTCGCTGGGGAAACCt    >  2:1165276/1‑93 (MQ=255)
                                                                                     tATTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTCGCTGGGGAAACCt    >  1:9178564/1‑93 (MQ=255)
                                                                                      aTTTGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTCGCTGGGGAAACCtg   <  2:3848828/93‑1 (MQ=255)
                                                                                       tttGCCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTCGCTGGGGAAACCtgt  >  1:4073020/1‑93 (MQ=255)
                                                                                          |                                                                                         
GGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTCGCTGGGGAAACCTGT  >  NC_000913/3968476‑3968655

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: