Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 999,571 A→G Q378R (CAG→CGG)  ycbS → putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913999,5710AG97.5% 370.9 / ‑0.9 119Q378R (CAG→CGG) ycbSputative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base A (3/0);  new base G (47/69);  total (50/69)
Fisher's exact test for biased strand distribution p-value = 7.16e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

Only 100 of 129 total aligned reads displayed.

CGTTCTCCAGTGTGCCATTACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCAGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAATATGCTGAATATTAC  >  NC_000913/999481‑999659
                                                                                          |                                                                                        
cGTTCTCCAGTGTGCCATTACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCggg                                                                                        >  1:921406/1‑93 (MQ=255)
 gTTCTCCAGTGTGCCATTACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGga                                                                                       >  1:3411276/1‑93 (MQ=255)
  ttCTCCAGTGTGCCATTACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGgaa                                                                                      >  1:2086192/1‑92 (MQ=255)
      ccAGTGTGCCATTACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCa                                                                                  <  2:181395/93‑1 (MQ=255)
      ccAGTGTGCCATTACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCa                                                                                  >  2:4284392/1‑93 (MQ=255)
       cAGTGTGCCATTACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCaa                                                                                 >  2:3235507/1‑93 (MQ=255)
       cAGTGTGCCATTACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCaa                                                                                 <  1:1109465/93‑1 (MQ=255)
       cAGTGTGCCATTACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCaa                                                                                 >  2:1266055/1‑93 (MQ=255)
         gtgtGCCATTACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAAt                                                                               >  2:26323/1‑93 (MQ=255)
         gtgtGCCATTACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAAt                                                                               >  2:2638486/1‑93 (MQ=255)
          tgtgCCATTACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAAtt                                                                              <  1:1211950/93‑1 (MQ=255)
          tgtgCCATTACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAAtt                                                                              <  2:2648535/93‑1 (MQ=255)
            tgCCATTACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGTGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTg                                                                            <  2:1153701/93‑1 (MQ=255)
            tgCCATTACTCCAACGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTg                                                                            <  1:4150365/93‑1 (MQ=255)
              ccATTACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCAGGAGAGCAAATTTGcc                                                                          >  1:3296981/1‑93 (MQ=255)
               cATTACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGccc                                                                         >  1:3762842/1‑93 (MQ=255)
               cATTACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGccc                                                                         <  2:3478721/93‑1 (MQ=255)
               cATTACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGccc                                                                         <  2:2399162/93‑1 (MQ=255)
               cATTACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGccc                                                                         <  2:3792949/93‑1 (MQ=255)
                aTTACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCa                                                                        <  1:719622/93‑1 (MQ=255)
                 ttACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAg                                                                       <  2:323093/93‑1 (MQ=255)
                  tACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAgg                                                                      <  1:182246/93‑1 (MQ=255)
                  tACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAgg                                                                      >  2:3042016/1‑93 (MQ=255)
                    cTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTTCCCAGGcc                                                                    <  2:3731838/93‑1 (MQ=255)
                       cAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACg                                                                 >  2:1689285/1‑93 (MQ=255)
                        aGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGt                                                                <  2:170782/93‑1 (MQ=255)
                        aGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGt                                                                <  1:1729734/93‑1 (MQ=255)
                        aGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGt                                                                >  2:589554/1‑93 (MQ=255)
                         gCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGtt                                                               <  1:606380/93‑1 (MQ=255)
                          cGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTg                                                              <  2:2235250/93‑1 (MQ=255)
                           gTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGc                                                             >  1:1142847/1‑93 (MQ=255)
                           gTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGc                                                             >  1:485037/1‑93 (MQ=255)
                             cAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAg                                                           >  1:3659972/1‑93 (MQ=255)
                             cAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAg                                                           <  2:1054516/93‑1 (MQ=255)
                             cAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAg                                                           <  1:2261630/93‑1 (MQ=255)
                             cAGGGGCGAATCAAATACGCGGAGACACTGGCGAAATACAGAACCAATAGTAATCAACAGCGGGAGAGCAAATTTGCCCAGGCCACGttctcg                                                           >  1:3659973/1‑89 (MQ=255)
                              aGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGt                                                          <  1:204167/93‑1 (MQ=255)
                              aGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGt                                                          >  1:957536/1‑93 (MQ=255)
                               ggggCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTg                                                         >  2:2817003/1‑93 (MQ=255)
                                gggCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTgg                                                        <  2:3839608/93‑1 (MQ=255)
                                  gCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTgggg                                                      <  2:1261637/93‑1 (MQ=255)
                                      aTCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGa                                                  >  2:2413516/1‑93 (MQ=255)
                                      aTCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGa                                                  <  2:321522/93‑1 (MQ=255)
                                       tCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGAc                                                 <  2:3660265/93‑1 (MQ=255)
                                       tCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGAc                                                 >  2:2440019/1‑93 (MQ=255)
                                          aaTACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCgg                                              >  1:978172/1‑92 (MQ=255)
                                          aaTACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGt                                              <  1:1998133/93‑1 (MQ=255)
                                          aaTACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGt                                              <  1:299695/93‑1 (MQ=255)
                                          aaTACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGt                                              >  1:2361495/1‑93 (MQ=255)
                                           aTAGGCGGTGACACTGGGAAAATACAGAAACAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTg                                             <  1:938411/93‑1 (MQ=255)
                                           aTACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTg                                             >  2:2971621/1‑93 (MQ=255)
                                            tACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTgg                                            >  2:562061/1‑93 (MQ=255)
                                            tACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTgg                                            >  1:1753200/1‑93 (MQ=255)
                                            tACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTgg                                            <  2:210545/93‑1 (MQ=255)
                                               gcgGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGa                                         >  1:754953/1‑93 (MQ=255)
                                                cgGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGaa                                        >  1:3524294/1‑93 (MQ=255)
                                                 ggTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAAc                                       <  1:2825403/93‑1 (MQ=255)
                                                  gTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACg                                      >  1:642978/1‑93 (MQ=255)
                                                     acacTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGAcc                                   <  1:1720095/93‑1 (MQ=255)
                                                      cacTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCt                                  <  1:1313445/93‑1 (MQ=255)
                                                      cacTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCt                                  >  1:1858067/1‑93 (MQ=255)
                                                        cTGGCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTgg                                <  1:781618/93‑1 (MQ=255)
                                                          ggCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTa                              <  2:3576203/93‑1 (MQ=255)
                                                          ggCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTa                              <  2:4614572/93‑1 (MQ=255)
                                                          ggCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTa                              <  1:3459394/93‑1 (MQ=255)
                                                           gCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTAt                             <  1:2082514/93‑1 (MQ=255)
                                                           gCGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTAt                             <  2:2793914/93‑1 (MQ=255)
                                                            cGAAATACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTAtg                            <  2:360905/93‑1 (MQ=255)
                                                              aaaTACAGAACCAATAGTAATGAACAGAGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTAtggt                          <  2:3050040/93‑1 (MQ=255)
                                                               aaTACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTAtggtg                         >  2:1383652/1‑93 (MQ=255)
                                                                aTACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTAtggtgg                        <  1:597435/93‑1 (MQ=255)
                                                                 tACAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTAtggtggg                       >  1:396592/1‑92 (MQ=255)
                                                                  aCAGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTAtggtggtg                      <  2:3322581/93‑1 (MQ=255)
                                                                    aGAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGa                    <  2:1947164/93‑1 (MQ=255)
                                                                     gAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGAc                   <  2:3848402/93‑1 (MQ=255)
                                                                     gAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGAc                   <  2:3460969/93‑1 (MQ=255)
                                                                     gAACCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGAc                   >  1:4231957/1‑93 (MQ=255)
                                                                       aCCAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACaa                 <  1:2691067/93‑1 (MQ=255)
                                                                        ccAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGCCAat                >  2:1068760/1‑93 (MQ=255)
                                                                         cAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAata               <  1:1275520/93‑1 (MQ=255)
                                                                         cAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAata               <  1:1830822/93‑1 (MQ=255)
                                                                         cAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAata               <  2:1474078/93‑1 (MQ=255)
                                                                         cAATAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAata               <  2:1319974/93‑1 (MQ=255)
                                                                          aaTAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAatat              <  2:1790121/93‑1 (MQ=255)
                                                                          aaTAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAatat              <  1:4161390/93‑1 (MQ=255)
                                                                          aaTAGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAatat              >  1:639341/1‑93 (MQ=255)
                                                                             aGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAATATGCt           <  2:30072/93‑1 (MQ=255)
                                                                             aGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAATATGCt           <  1:1267672/93‑1 (MQ=255)
                                                                             aGTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAATATGCt           <  2:1333738/93‑1 (MQ=255)
                                                                              gTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAATATGCTg          >  2:3883009/1‑93 (MQ=255)
                                                                              gTAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAATATGCTg          <  1:1133552/93‑1 (MQ=255)
                                                                               tAATGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAATATGCTGa         <  1:1708006/93‑1 (MQ=255)
                                                                                aaTGAACAGCAGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAATACGCTGaa        >  1:4127116/1‑93 (MQ=255)
                                                                                 aTGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTTGACAATATGCTGAat       >  1:2957241/1‑93 (MQ=255)
                                                                                 aTGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAATATGCTGAat       <  2:3130635/93‑1 (MQ=255)
                                                                                 aTGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGCATGGTGGGGGACAATATGCTGAat       >  2:1813141/1‑93 (MQ=255)
                                                                                  tGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAATATGCTGAata      <  1:4517211/93‑1 (MQ=255)
                                                                                  tGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAATATGCTGAata      >  2:4306237/1‑93 (MQ=255)
                                                                                  tGAACAGCGGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAATATGCTGAata      >  2:2715591/1‑93 (MQ=255)
                                                                                   gacgagcggGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAATATGCTGAatat     <  2:4155227/89‑1 (MQ=255)
                                                                                          |                                                                                        
CGTTCTCCAGTGTGCCATTACTCCAGCGTCAGGGGCGAATCAAATACGCGGTGACACTGGCGAAATACAGAACCAATAGTAATGAACAGCAGGAGAGCAAATTTGCCCAGGCCACGTTGCAGTGGGGCGGACCGTGGGGAACGACCTGGTATGGTGGTGGACAATATGCTGAATATTAC  >  NC_000913/999481‑999659

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