Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,363,167 T→G L121V (TTA→GTA)  yhcD → putative outer membrane fimbrial subunit usher protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,363,1670TG95.1% 593.3 / 8.1 187L121V (TTA→GTA) yhcDputative outer membrane fimbrial subunit usher protein
Reads supporting (aligned to +/- strand):  ref base T (8/1);  new base G (94/82);  total (102/85)
Fisher's exact test for biased strand distribution p-value = 4.33e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 187 total aligned reads displayed.

GGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGTACCACAATCACACCATTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGATCTCACAAACGCAGCAACTGGTGGAACAGCTG  >  NC_000913/3363070‑3363263
                                                                                                 |                                                                                                
 gCATTTCTGCACCAAGCCGGTTTAAAAATTCGTTCTGAAGAAAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGTACCACAATCACACCAGTAc                                                                                               <  1:2440830/100‑1 (MQ=255)
  cagttCTGCAACAAGCCGGTTTAAAAATTAAAAAAAAAAAAAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGTACCACAATCACACCAGTAcc                                                                                              <  2:3616108/97‑1 (MQ=255)
  cATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGTACCACAATCACACCAGTAcc                                                                                              <  2:300750/100‑1 (MQ=255)
   aTTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGTACCACAATCACACCAGTACCt                                                                                             >  1:1145163/1‑100 (MQ=255)
    tttCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGTACCACAATCACACCAGTACCTa                                                                                            >  1:3441135/1‑100 (MQ=255)
       cTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAGAATGGCTGTTACGACTATATATTGTCCTACCCGGGTACCACAATCACACCAGTACCTAAcc                                                                                         >  2:40614/1‑100 (MQ=255)
       cTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGTACCACAATCACACCAGTACCTAAcc                                                                                         >  2:42474/1‑100 (MQ=255)
          caacaaGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGTACCACAATCACACCAGTACCTAACCAgg                                                                                      <  1:1804230/100‑1 (MQ=255)
             taaGTCGGTGTAAAAATTCCTTCTGAAGAAAAAAAGGGGTGTTACGATTATATATTGTCCTACCCGGGTATCACAATCACACCAGTACCTAACCAGGAAg                                                                                   <  2:892519/99‑1 (MQ=255)
              aaGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGTACCACAATCACACCAGTACCTAACCAGGAAGc                                                                                  >  1:1488654/1‑100 (MQ=255)
                gCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGTACCACAATCACACCAGTACCTAACCAGGAAGCGt                                                                                >  2:2277587/1‑100 (MQ=255)
                 ccGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGTACCACAATCACACCAGGACCTAGCCAGGAAGCGtt                                                                               <  2:3352852/100‑1 (MQ=255)
                   ggTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGTACCACAATCACACCAGTACCTAACCAGGAAGCGTTAg                                                                             >  2:2126667/1‑100 (MQ=255)
                   ggTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGTACCACAATCACACCAGTACCTAACCAGGAAGCGTTAg                                                                             <  1:1504631/100‑1 (MQ=255)
                     tttAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACTATATATTGTCCTACACGGGTACCACAATCACACCAGTACCTACCCAGGACGCGTTAGat                                                                           >  2:1936903/1‑100 (MQ=255)
                      ttAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGTACCACAATCACACCAGTACCTAACCAGGAAGCGTTAGata                                                                          <  2:575140/100‑1 (MQ=255)
                      ttAAAAATTCCTACTGAAGAAAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGTACCACAATCACACCAGTACCTAACCAGGAAGCGTTAGata                                                                          >  1:3517003/1‑100 (MQ=255)
                           aaTTCCTTCTGAAGAAAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGTACCACAATCACACCAGTACCTAACCAGGAAGCGTTAGATATTATc                                                                     <  1:3170400/100‑1 (MQ=255)
                            aTTCCTTCTGAAGAAAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGTACCACAATCACACCAGTACCTAACCAGGAAGCGTTAGATATTATCg                                                                    <  2:2472837/100‑1 (MQ=255)
                               ccTTCTGAAGAAAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGTACCACAATCACACCAGTACCTAACCAGGAAGCGTTAGATATTATCGttt                                                                 >  1:1754880/1‑100 (MQ=255)
                               ccTTCTGAAGAAAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGTACCACAATCACACCAGTACCTAACCAGGAAGCGTTAGATATTATCGttt                                                                 <  2:1106616/100‑1 (MQ=255)
                                 ttCTGAAGAAAAAAATGGCTGTTGCGACTATATATTGTCCTACGCGGGTATCACAATCACACCAGTACCTGACCAGGAAGCGTTAGTTATTACCGTTTca                                                               <  2:3299384/100‑1 (MQ=255)
                                 ttCTGAAGAAAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGTACCACAATCACACCAGTACCTAACCAGGAAGCGGTAGATATTATCGTTTca                                                               >  1:190041/1‑100 (MQ=255)
                                  tCTGAAGAAAAAAATGGCTGTTACGACTATATATTGTACTACCCGGGTACCTCAATCACCCCAGTACCTAACCAGGAAGCGTTAGATATTATCGTTTcac                                                              <  2:3628582/100‑1 (MQ=255)
                                   cTGAAGAAAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGTACCACAATCACACCAGTACCTAACCAGGAAGCGTTAGATATTATCGTTTcacc                                                             >  1:1954260/1‑100 (MQ=255)
                                    tGAAGAAAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGTACCACAATCACACCAGTACCTAACCAGGAAGCGTTAGATATTATCGTTTcacca                                                            <  2:143491/100‑1 (MQ=255)
                                     gaagaaAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGTACCACAATCACACCAGTACCTAACCAGGAAGCGTTAGATATTATCGTTTcaccac                                                           >  1:3826003/1‑100 (MQ=255)
                                      aagaaAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGTACCACAATCACACCATTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACcaca                                                          >  2:1244953/1‑100 (MQ=255)
                                       agaaAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGTACCACAATCACACCAGTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAg                                                         >  1:3334705/1‑100 (MQ=255)
                                        gaaAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGTACCACAATCACACCAGTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAgg                                                        <  2:3245591/100‑1 (MQ=255)
                                        gaaAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGTACCACAATCACACCAGTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAgg                                                        >  1:410800/1‑100 (MQ=255)
                                        gaaAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGTACCACAAGCACACCAGTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAgg                                                        >  1:705370/1‑100 (MQ=255)
                                         aaaaaaaTGGCTGTTACGACTATATATTGTCCTACCCGGGTACCACAATGACACCACTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGc                                                       <  1:3955001/100‑1 (MQ=255)
                                         aaaaaaaTGGCTGTTACGACTATATATTGTCCTACCCGGGTACCACAATCACACCAGTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAgcc                                                       >  2:225441/1‑98 (MQ=255)
                                         aaaaaaaTGGCTGTTACGACTATATATTGTCCTACCCGGGTACCACAATCACACCAGTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGc                                                       <  1:398506/100‑1 (MQ=255)
                                         aaaaaaaTGGCTGTTACGACTATATATTGTCCTACCCGGGTACCACAATCACACCAGTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGc                                                       <  2:293747/100‑1 (MQ=255)
                                         aaaaaaaTGGCTGTTACGACTATATATTGTCCTACCCGGGTACCACAATCACACCAGTACCTAACCAGGAAGCGTTAAATATTATCGTTTCACCACAGGc                                                       <  2:1095007/100‑1 (MQ=255)
                                         aaaaaaaTGGCTGTTACGACTATACAGTGTCCTACCCGGGAACCACAACCCCACCAGCACCTAACCAGGAAGCGATAGATATTATCGTTTCACCACAGGc                                                       >  2:1774808/1‑100 (MQ=255)
                                              aaTGGCTGTTACGACTATATATTGTCCTACCCGGGTACCACAATCACACCAGTACCTAACCAGGAAGCGTTAGGTATTATCGTTTCACCACAGGCGATCa                                                  <  2:601123/100‑1 (MQ=255)
                                               cgatcTGTTACGACTATATATTGTCCTACCCGGGTACCACAATCACACCAGTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCAt                                                 <  2:2035724/96‑1 (MQ=255)
                                               aTGGCTGTTACGACTATATATTGTCCTACCCGGGTACCACAATCACACCAGTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCAt                                                 >  1:1603216/1‑100 (MQ=255)
                                                  gCTGTTACGACTATATATTGTCCTACCCGGGTACCACAATCACACCAGTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTcc                                              <  1:3393494/100‑1 (MQ=255)
                                                   cTGTTACGACTATATATTGTCCTACCCGGGTACCACAATCACACCAGTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTccc                                             >  1:2609446/1‑100 (MQ=255)
                                                     gTTACGACTATATATTGTCCTACCCGGGTACCACAATCACACCAGTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCAt                                           >  1:322287/1‑100 (MQ=255)
                                                     gTTACGACTATATATTGTCCTACCCGGGTACCACAATCACACCAGTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCAt                                           >  1:550007/1‑100 (MQ=255)
                                                      ttACGACTATATATTGTCCTACCCGGGTACCACAATCACACCAGTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATa                                          <  2:3884506/100‑1 (MQ=255)
                                                        aCGACTATATATTGTCCTACCCGGGTACCACAATCACACCAGTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAgg                                        <  1:2582479/100‑1 (MQ=255)
                                                           aCTATATATTGTCCTACCCGGGTACCACAATCACACCAGTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGtt                                     >  2:2834838/1‑100 (MQ=255)
                                                            cTATATATTGTCCTACCCGGGTACCACAATCACACCAGTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTg                                    <  2:1301436/100‑1 (MQ=255)
                                                             tatataTTGTCCTACCCGGGTACCACAATCACACCAGTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTgg                                   <  1:993220/100‑1 (MQ=255)
                                                               tatatTGTCCTACCCGGGTACCACAATCACACCAGTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGAt                                 >  1:66403/1‑100 (MQ=255)
                                                                atatTGTCCTACCCGGGTACCACAATCACACCAGTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGAtc                                <  1:2156775/100‑1 (MQ=255)
                                                                 tatTGTCCTACCCGGGTACCACAATCACACCATTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGAtct                               >  1:1137896/1‑100 (MQ=255)
                                                                 tatTGTCCTACCCGGGTACCACAATCACACCAGTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGAtgt                               >  1:2696500/1‑98 (MQ=255)
                                                                 tatTGTCCTACCCGGGTACCACAATCACACCAGTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGAtct                               >  1:1444690/1‑100 (MQ=255)
                                                                 tatTGTCCTACCCGGGTACCACAATCACACCAGTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGAtct                               >  2:2594421/1‑100 (MQ=255)
                                                                 tatTGTCCTACCCGGGTACCACAATCACACCAGTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGAtct                               <  2:1923908/100‑1 (MQ=255)
                                                                 tatTGTCCTACCCGGGTACCACAATCACACCAGTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGAtct                               >  1:3718916/1‑100 (MQ=255)
                                                                  atTGTCCTACCCGGGTACCACAATCACACCAGTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGAtctc                              >  1:1062305/1‑100 (MQ=255)
                                                                  atTGTCCTACCCGGGTACCACAATCACACCAGTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGAtctc                              >  1:2392649/1‑100 (MQ=255)
                                                                  atTGTCCTACCCGGGTACCACAATCACACCAGTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGAtctc                              <  1:2720665/100‑1 (MQ=255)
                                                                   gtGTCCTACCCGGGTACCACAATCACACCAGTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGATCTca                             <  2:690599/99‑1 (MQ=255)
                                                                    tGTCCTACCCGGGTACCACAATCACACCAGTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGATCTcac                            <  1:910835/100‑1 (MQ=255)
                                                                    tGTCCTACCCGGGTACCACAAAAAAAAAAAAAAAAAAAAAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGATCTcac                            <  2:1559491/100‑1 (MQ=255)
                                                                     gTCCTACCCGGGTACCACAATCACACCAGTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGATCTcaca                           >  2:1421448/1‑100 (MQ=255)
                                                                       ccTACCCGGGTACCACAATCACACCAGTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGATCTCACaaa                         <  2:717250/100‑1 (MQ=255)
                                                                       ccTACCCGGGTACCACAATCACACCAGTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGATCTCACaaa                         >  1:1288911/1‑100 (MQ=255)
                                                                        cTACCCGGGTACCACAATCACACCATTACCTACCCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGATCTCACAAAc                        <  2:957823/100‑1 (MQ=255)
                                                                        cTACCCGGGTACCACAATCACACCAGTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGATCTCACAAAc                        <  2:2789052/100‑1 (MQ=255)
                                                                         tACCCGGGTACCACAATCACACCAGTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGATCTCACAAACg                       >  1:2274406/1‑100 (MQ=255)
                                                                         tACCCGGGTACCACAATCACACCAGTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGATCTCACAAACg                       >  1:36199/1‑100 (MQ=255)
                                                                          aCCCGGGTACCACAATCACACCAGTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGATCTCACAAACgc                      <  1:1591071/100‑1 (MQ=255)
                                                                           cccGGGTACCACAATCACACCATTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGATCTCACAAACgca                     >  2:2009375/1‑100 (MQ=255)
                                                                           cccGGGTACCACAATCACACCAGTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGATCTCACAAACgca                     <  1:2904226/100‑1 (MQ=255)
                                                                            ccGGGTACCACAATCACACCAGTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGATCTCACAAACgcag                    >  1:1113164/1‑100 (MQ=255)
                                                                            ccGGGTACCACAATCACACCAGTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGATCTCACAAACgcag                    <  1:3997922/100‑1 (MQ=255)
                                                                              gggTACCACAATCACACCAGTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGATCTCACAAACgcagca                  <  1:759361/100‑1 (MQ=255)
                                                                                gTACCACAATCACACCAGTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGATCTCACAAACGCAGCAAc                <  2:1488958/100‑1 (MQ=255)
                                                                                 tACCACAATCACACCAGTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGATCTCACAAACGCAGCAACt               <  2:838387/100‑1 (MQ=255)
                                                                                 tACCACAATCACACCAGTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGATCTCACAAACGCAGCAACt               <  1:1566625/100‑1 (MQ=255)
                                                                                   ccACAATCACACCAGTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCCTTCCCATAGGGTTGGATCTCCCAAACGCAGCAACtgg             >  2:3082955/1‑100 (MQ=255)
                                                                                   ccACAATCACACCAGTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGATCTCACAAACGCAGCAACtgg             >  2:2449205/1‑100 (MQ=255)
                                                                                   ccACAATCACACCAGTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGATCTCACAAACGCAGCAACtgg             >  1:2458212/1‑100 (MQ=255)
                                                                                    cacaATCACACCAGTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGATCTCACAAACGCAGCAACtggt            >  2:3824160/1‑100 (MQ=255)
                                                                                    cacaATCACACCAGTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGATCTCACAAACGCAGCAACtggt            >  1:4069987/1‑100 (MQ=255)
                                                                                     acaATCACACCAGTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGATCTCACAAACGCAGCAACtggtg           >  1:3032021/1‑100 (MQ=255)
                                                                                     acaATCACACCAGTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGATCTCACAAACGCAGCAACtggtg           <  2:656656/100‑1 (MQ=255)
                                                                                      caATCACACCAGTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGATCTCACAAACGCAGCAACtggtgg          >  1:2900146/1‑100 (MQ=255)
                                                                                      caATCACACCAGTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGATCTCACAAACGCAGCAACtggtgg          <  1:1979422/100‑1 (MQ=255)
                                                                                        aTCACACCATTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGATCTCACAAACGCAGCAACTGGTGGaa        >  2:303845/1‑100 (MQ=255)
                                                                                         tCACACCAGTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGATCTCACAAACGCAGCAACTGGTGGAAc       <  1:229959/100‑1 (MQ=255)
                                                                                          cacaCCAGTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGATCTCACAAACGCAGCAACTGGTGGAACa      >  2:2524799/1‑100 (MQ=255)
                                                                                           ccacCAGTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGATCTCACAAACGCAGCAACTGGTGGAACAg     <  1:868883/99‑1 (MQ=255)
                                                                                           acacCATTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGATCTCACAAACGCAGCAACTGGTGGAACAg     >  2:3681429/1‑100 (MQ=255)
                                                                                           acacCAGTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGATCTCACAAACGCAGCAACTGGTGGAACAg     >  1:1267883/1‑100 (MQ=255)
                                                                                            cacCAGTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGATCTCACAAACGCAGCAACTGGTGGAACAgc    >  2:3674537/1‑100 (MQ=255)
                                                                                             acCATTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGATCTCACAAACGCAGCAACTGGTGGAACAgct   >  1:278350/1‑100 (MQ=255)
                                                                                             acCAGTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGATCTCACAAACGCAGCAACTGGTGGAACAgct   <  2:2355907/100‑1 (MQ=255)
                                                                                              ccAGTACCTAACCAGGAAGCGTTAGATGTTATCGTTTCACCACAGGCGTTCATTACCATAGGGTTGGATCTCACAAACGCAGCAACTGGTGGAACAgcgg  >  2:3984949/1‑98 (MQ=255)
                                                                                              ccAGTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGATCTCACAAACGCAGCAACTGGTGGAACAgctg  >  2:3860943/1‑100 (MQ=255)
                                                                                                 |                                                                                                
GGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAAAAAATGGCTGTTACGACTATATATTGTCCTACCCGGGTACCACAATCACACCATTACCTAACCAGGAAGCGTTAGATATTATCGTTTCACCACAGGCGATCATTCCCATAGGGTTGGATCTCACAAACGCAGCAACTGGTGGAACAGCTG  >  NC_000913/3363070‑3363263

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