Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,375,376 A→G intergenic (‑76/+13) blc ← / ← ampC outer membrane lipoprotein/beta‑lactamase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,375,3760AG87.0% 755.6 / 94.9 308intergenic (‑76/+13)blc/ampCouter membrane lipoprotein/beta‑lactamase
Reads supporting (aligned to +/- strand):  ref base A (15/25);  new base G (124/144);  total (139/169)
Fisher's exact test for biased strand distribution p-value = 3.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.02e-01

Only 100 of 309 total aligned reads displayed.

GAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGA  >  NC_000913/4375285‑4375468
                                                                                            |                                                                                            
    gAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGGCCCGa                                                                                         <  2:6221293/93‑1 (MQ=255)
      gCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGGCCCGATg                                                                                       <  2:7152955/93‑1 (MQ=255)
       cAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATgg                                                                                      <  2:3325626/93‑1 (MQ=255)
         ggCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGGCCCGATGGaa                                                                                    >  1:6548693/1‑93 (MQ=255)
         ggCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGGCCCGATGGaa                                                                                    <  2:1395365/93‑1 (MQ=255)
         ggCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGGCCCGATGGaa                                                                                    >  2:2921525/1‑93 (MQ=255)
         ggCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGGCCCGATGGaa                                                                                    <  1:5099906/93‑1 (MQ=255)
         ggCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGGCCCGATGGaa                                                                                    >  2:2739698/1‑93 (MQ=255)
          gcgcATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGGCCCGATGGAAt                                                                                   >  1:3635430/1‑93 (MQ=255)
           cgcATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGGCCCGATGGAAtt                                                                                  <  2:1219019/93‑1 (MQ=255)
            gcATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGGCCCGATGGAAttt                                                                                 <  2:3542627/93‑1 (MQ=255)
            gcATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGGCCCGATGGAAttt                                                                                 <  1:6490065/93‑1 (MQ=255)
            gcATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCAGGCCCGATGGAAttt                                                                                 >  2:4092674/1‑93 (MQ=255)
             cATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGGCCCGATGGAAtttt                                                                                <  2:1083681/93‑1 (MQ=255)
               tAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGGCCCGATGGAATTTTAc                                                                              <  2:6391348/93‑1 (MQ=255)
               tAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGGCCCGATGGAATTTTAc                                                                              >  1:1443489/1‑93 (MQ=255)
               tAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTAc                                                                              <  2:1836872/93‑1 (MQ=255)
                aaaTGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGGCCCGATGGAATTTTACt                                                                             <  2:2822685/93‑1 (MQ=255)
                 aaTGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGGCCCGATGGAATTTTACTg                                                                            <  2:5293284/93‑1 (MQ=255)
                   tGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGGCCCGATGGAATTTTACTGTa                                                                          >  2:494495/1‑93 (MQ=255)
                    gTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGGCCCGATGGAATTTTACTGTag                                                                         <  2:3959032/93‑1 (MQ=255)
                     gtttcTTTCTGGTTTTTTTCCAAATGTTGGAGACAAGAAAGAAAAGCCGAAGAAAAGTCCGAAAAATTCGGGCCCGATGGAATTTTACTGTag                                                                         <  2:3147242/89‑1 (MQ=255)
                          ttACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGGCCCGATGGAATTTTACTGTAGAGCGtt                                                                   >  1:6788624/1‑93 (MQ=255)
                          ttACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGGCCCGATGGAATTTTACTGTAGAGCGtt                                                                   >  1:4330662/1‑93 (MQ=255)
                          ttACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGtt                                                                   <  2:5080626/93‑1 (MQ=255)
                           tACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGGCCCGATGGAATTTTACTGTAGAGCGTTa                                                                  >  1:996383/1‑93 (MQ=255)
                           tACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGGCCCGATGGAATTTTACTGTAGAGCGTTa                                                                  <  2:5195267/93‑1 (MQ=255)
                            aCTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGGCCCGATGGAATTTTACTGTAGAGCGTTaa                                                                 <  2:3076132/93‑1 (MQ=255)
                            aCTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTaa                                                                 <  2:2600220/93‑1 (MQ=255)
                             cTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGGCCCGATGGAATTTTACTGTAGAGCGTTAAg                                                                >  1:2173580/1‑93 (MQ=255)
                              tGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGGCCCGATGGAATTTTACTGTAGAGCGTTAAGa                                                               <  2:745146/93‑1 (MQ=255)
                              tGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGGCCCGATGGAATTTTACTGTAGAGCGTTAAGa                                                               <  1:5332019/93‑1 (MQ=255)
                                 tttttttCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGGCCCGATGGAATTTTACTGTAGAGCGTTAAGAATc                                                            <  1:7040640/93‑1 (MQ=255)
                                 tttttttCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGGCCCGATGGAATTTTACTGTAGAGCGTTAAGAATc                                                            >  1:5681002/1‑93 (MQ=255)
                                  ttttttCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGGCCCGATGGAATTTTACTGTAGAGCGTTAAGAATCt                                                           >  1:3377269/1‑93 (MQ=255)
                                  ttttttCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGGCCCGATGGAATTTTACTGTAGAGCGTTAAGAATCt                                                           <  1:2630101/93‑1 (MQ=255)
                                   tttttCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGGCCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTg                                                          <  1:3170876/93‑1 (MQ=255)
                                   tttttCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGGCCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTg                                                          <  2:379009/93‑1 (MQ=255)
                                    ttttCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGGCCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGc                                                         <  2:2169926/93‑1 (MQ=255)
                                    ttttCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGGCCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGc                                                         >  2:5467261/1‑93 (MQ=255)
                                     tttCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGGCCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGcc                                                        <  1:6389269/93‑1 (MQ=255)
                                      ttCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGGCCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCa                                                       <  1:4186691/93‑1 (MQ=255)
                                      ttCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGGCCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCa                                                       >  2:6638634/1‑93 (MQ=255)
                                       tCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGGCCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAg                                                      <  1:1256911/93‑1 (MQ=255)
                                       tCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGGCCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAg                                                      <  2:2226488/93‑1 (MQ=255)
                                        cTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGGCCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAgg                                                     <  2:3889794/93‑1 (MQ=255)
                                        cTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGGCCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAgg                                                     <  1:3692112/93‑1 (MQ=255)
                                        cTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGGCCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAgg                                                     <  1:4006265/93‑1 (MQ=255)
                                        cTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAgg                                                     >  1:6660660/1‑93 (MQ=255)
                                         tAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGGCCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAggc                                                    <  1:3093206/93‑1 (MQ=255)
                                         tAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGGCCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAggc                                                    >  2:157167/1‑93 (MQ=255)
                                           aGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGGCCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAggcgg                                                  >  2:4475027/1‑93 (MQ=255)
                                            gtgtAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGGCCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAggcggc                                                 >  2:1587700/1‑93 (MQ=255)
                                            gtgtAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAggcggc                                                 <  1:792909/93‑1 (MQ=255)
                                             tgtAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGGCCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGgcggcg                                                >  1:570813/1‑93 (MQ=255)
                                              gtAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGGCCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGt                                               <  2:5819758/93‑1 (MQ=255)
                                              gtAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGGCCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGt                                               <  2:4585802/93‑1 (MQ=255)
                                               tAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTc                                              <  1:5933970/93‑1 (MQ=255)
                                                 gATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGGCCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGa                                            >  2:3262750/1‑93 (MQ=255)
                                                  aTGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGGCCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGAc                                           <  1:3890542/93‑1 (MQ=255)
                                                  aTGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGGCCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGAc                                           <  2:6313371/93‑1 (MQ=255)
                                                    gACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGGCCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGActc                                         <  1:3045673/93‑1 (MQ=255)
                                                     aCAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGActct                                        >  1:183700/1‑93 (MQ=255)
                                                      cAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACtctc                                       <  1:4093881/93‑1 (MQ=255)
                                                       aGCAAGGAAAAGCGGAGAAAAGGTCCAAAAATTAGGGCCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCg                                      >  1:4983454/1‑93 (MQ=255)
                                                        gCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGGCCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGc                                     >  2:5496601/1‑93 (MQ=255)
                                                         cAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGGCCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCt                                    <  1:3340307/93‑1 (MQ=255)
                                                          aaGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGGCCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTg                                   <  1:1076478/93‑1 (MQ=255)
                                                            ggAAAAGCGGAGAAAAGGTCCGAAAATTCGGGCCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTggc                                 >  1:2420711/1‑92 (MQ=255)
                                                             gAAAAGCGGAGAAAAGGTCCGAAAATTCGGGCCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGAt                                <  1:905441/93‑1 (MQ=255)
                                                             gAAAAGCGGAGAAAAGGTCCGAAAATTCGGGCCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGAt                                <  1:3980991/93‑1 (MQ=255)
                                                             gAAAAGCGGAGAAAAGGTCCGAAAATTCGGGCCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGAt                                >  1:132756/1‑93 (MQ=255)
                                                              aaaaGCGGAGAAAAGGTCCGAAAATTCGGGCCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGAtt                               >  1:2220805/1‑93 (MQ=255)
                                                              aaaaGCGGAGAAAAGGTCCGAAAATTCGGGCCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGAtt                               <  2:2128609/93‑1 (MQ=255)
                                                               aaaGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTg                              >  1:1499137/1‑93 (MQ=255)
                                                                aaGCGGAGAAAAGGTCCGAAAATTCGGGCCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTgg                             >  2:564826/1‑93 (MQ=255)
                                                                aaGCGGAGAAAAGGTCCGAAAATTCGGGCCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTgg                             <  1:175536/93‑1 (MQ=255)
                                                                aaGCGGAGAAAAGGTCCGAAAATTCGGGCCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTgg                             <  1:5218566/93‑1 (MQ=255)
                                                                  gCGGAGAAAAGGTCCGAAAATTCGGGCCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGa                           <  1:478384/93‑1 (MQ=255)
                                                                    ggAGAAAAGGTCCGAAAATTCGGGCCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATa                         <  1:7143014/93‑1 (MQ=255)
                                                                       gaAAAGGTCCGAAAATTCGGGCCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGtt                      <  2:808386/93‑1 (MQ=255)
                                                                       gaAAAGGTCCGAAAATTCGGGCCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGtt                      <  1:5789315/93‑1 (MQ=255)
                                                                          aaGGTCCGAAAATTCGGGCCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGttttt                   >  2:4236324/1‑93 (MQ=255)
                                                                             gTCCGAAAATTCGGGCCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTtttgtt                <  2:4866876/93‑1 (MQ=255)
                                                                              tCCGAAAATTCGGGCCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTtttgttt               >  1:5610644/1‑93 (MQ=255)
                                                                              tCCGAAAATTCGGGCCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTtttgttt               >  1:4244388/1‑93 (MQ=255)
                                                                              tCCGAAAATTCGGGCCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTtttgttt               >  2:4513378/1‑93 (MQ=255)
                                                                              tCCGAAAATTCGGGCCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTtttgttt               >  1:623701/1‑93 (MQ=255)
                                                                                cGAAAATTCGGGCCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGc             >  1:611604/1‑93 (MQ=255)
                                                                                 gAAAATTCGGGCCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGcc            >  2:2007682/1‑93 (MQ=255)
                                                                                   aaaaTTCGGGCCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCa           >  2:1168656/1‑93 (MQ=255)
                                                                                   aaaaTTCGGGCCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCa           <  2:7048746/93‑1 (MQ=255)
                                                                                    aaaTTCGGGCCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAg          >  1:5510108/1‑93 (MQ=255)
                                                                                    aaaTTCGGGCCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAg          >  1:7510598/1‑93 (MQ=255)
                                                                                     aaTTCGGGCCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGc         <  2:4646209/93‑1 (MQ=255)
                                                                                      aTTCGGGCCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCa        <  1:1394652/93‑1 (MQ=255)
                                                                                       ttCGGGCCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCAt       >  2:6014754/1‑93 (MQ=255)
                                                                                       ttCGGGCCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCAt       >  2:2050988/1‑93 (MQ=255)
                                                                                        tCGGGCCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATc      <  1:7027832/93‑1 (MQ=255)
                                                                                         cGGGCCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCa     <  1:4620302/93‑1 (MQ=255)
                                                                                            |                                                                                            
GAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGA  >  NC_000913/4375285‑4375468

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