Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,587,936 A→G R30R (CGT→CGC b1505 ← putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,587,9360AG99.5% 573.5 / ‑4.9 184R30R (CGT→CGCb1505putative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (90/93);  total (91/93)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.52e-01

Only 100 of 195 total aligned reads displayed.

TGACGTGTTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCACGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCTTGACGGTGTAACCACTCAT  >  NC_000913/1587846‑1588025
                                                                                          |                                                                                         
tGACGTGTTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACcgcg                                                                                         <  1:4597752/93‑1 (MQ=255)
tGACGTGTCCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACcgcg                                                                                         >  2:1427555/1‑93 (MQ=255)
  aCGTGTTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCt                                                                                       <  1:4195530/93‑1 (MQ=255)
     tgtTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTAc                                                                                    <  1:364285/93‑1 (MQ=255)
     tgtTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTAc                                                                                    <  1:689083/93‑1 (MQ=255)
      gtTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACt                                                                                   >  2:5159540/1‑93 (MQ=255)
      gtTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACt                                                                                   <  2:4592233/93‑1 (MQ=255)
        tCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGt                                                                                 <  1:4335027/93‑1 (MQ=255)
        tCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGt                                                                                 <  1:584903/93‑1 (MQ=255)
         ccAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTa                                                                                <  2:6935252/93‑1 (MQ=255)
           aGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGa                                                                              <  1:2973133/93‑1 (MQ=255)
           aGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGa                                                                              >  2:653917/1‑93 (MQ=255)
            gTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGaa                                                                             >  2:5492096/1‑93 (MQ=255)
              aaCTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACa                                                                           <  1:3906602/93‑1 (MQ=255)
                cTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGa                                                                         <  2:1236165/93‑1 (MQ=255)
                cTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGa                                                                         >  2:3116189/1‑93 (MQ=255)
                cTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGa                                                                         <  2:2622763/93‑1 (MQ=255)
                cTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGa                                                                         <  1:7532028/93‑1 (MQ=255)
                cTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGa                                                                         >  2:1490682/1‑93 (MQ=255)
                cTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGa                                                                         <  1:3573870/93‑1 (MQ=255)
                  ttGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGAtt                                                                       >  2:3281108/1‑93 (MQ=255)
                    ggCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTaa                                                                     >  2:7317012/1‑93 (MQ=255)
                      cGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTaaaa                                                                   <  1:108683/93‑1 (MQ=255)
                       gACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAaat                                                                  <  1:6179946/93‑1 (MQ=255)
                       gACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAaat                                                                  >  1:1722669/1‑93 (MQ=255)
                       gACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAaat                                                                  >  1:7891956/1‑93 (MQ=255)
                       gACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAaat                                                                  <  2:4412964/93‑1 (MQ=255)
                         cTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAaataa                                                                <  2:5094759/93‑1 (MQ=255)
                         cTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAaataa                                                                <  1:2770546/93‑1 (MQ=255)
                         cTGGCACTGAAAAATGTCGTACCATAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAaataa                                                                <  2:5649855/93‑1 (MQ=255)
                           ggCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATc                                                              >  2:1193468/1‑93 (MQ=255)
                            gCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCa                                                             >  2:2111779/1‑93 (MQ=255)
                             cACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCaa                                                            <  1:4246823/93‑1 (MQ=255)
                              aCTGAAAAATGTCGTACCGTCATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTCCTGTAGAACAGATTAAAATAATCAAt                                                           >  1:736217/1‑93 (MQ=255)
                              aCTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAAt                                                           <  2:3123450/93‑1 (MQ=255)
                               cTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATa                                                          <  2:1272225/93‑1 (MQ=255)
                               cTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATa                                                          <  2:5428610/93‑1 (MQ=255)
                                tGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATaa                                                         <  1:2861598/93‑1 (MQ=255)
                                tGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATaa                                                         >  2:7525315/1‑93 (MQ=255)
                                 gAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATaaa                                                        <  2:3265188/93‑1 (MQ=255)
                                 gAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATaaa                                                        >  1:3327750/1‑93 (MQ=255)
                                 gAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATaaa                                                        <  2:5811375/93‑1 (MQ=255)
                                 gAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATaaa                                                        <  1:398228/93‑1 (MQ=255)
                                  aaaaaTGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAAt                                                       <  2:1015844/93‑1 (MQ=255)
                                  aaaaaTGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAAt                                                       <  1:3593878/93‑1 (MQ=255)
                                   aaaaTGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAAtt                                                      <  1:6542569/93‑1 (MQ=255)
                                   aaaaTGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAAtt                                                      <  1:7384379/93‑1 (MQ=255)
                                   aaaaTGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAAtt                                                      <  1:2177949/93‑1 (MQ=255)
                                     aaTGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATtgt                                                    <  2:233607/93‑1 (MQ=255)
                                     aaTGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATtgt                                                    <  1:956003/93‑1 (MQ=255)
                                       tgttgtAACGTCCTTTGCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTgtgt                                                  <  2:6821662/93‑1 (MQ=255)
                                       tGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTgtgt                                                  >  1:5397678/1‑93 (MQ=255)
                                        gTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTc                                                 <  2:1106020/93‑1 (MQ=255)
                                          cGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTg                                               <  2:1768040/93‑1 (MQ=255)
                                          cGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTg                                               <  2:1798773/93‑1 (MQ=255)
                                            tACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCt                                             >  1:376131/1‑93 (MQ=255)
                                              ccGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCa                                           >  1:6358769/1‑93 (MQ=255)
                                               cGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCAt                                          >  2:1991469/1‑93 (MQ=255)
                                                gTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCAtt                                         <  1:1005519/93‑1 (MQ=255)
                                                gTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCAtt                                         <  2:4135249/93‑1 (MQ=255)
                                                gTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCAtt                                         <  2:4135247/93‑1 (MQ=255)
                                                gTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCAtt                                         >  2:7198822/1‑93 (MQ=255)
                                                gTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCAtt                                         >  1:7660958/1‑93 (MQ=255)
                                                gTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCAtt                                         >  1:5722453/1‑93 (MQ=255)
                                                 tAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGATCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTg                                        >  2:4908264/1‑93 (MQ=255)
                                                 tAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTg                                        >  2:3898561/1‑93 (MQ=255)
                                                   aTTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGgt                                      >  1:4950904/1‑93 (MQ=255)
                                                   aTTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGCGTCTGCTCATTGgt                                      <  1:7096701/93‑1 (MQ=255)
                                                    tttCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGgtg                                     >  1:660634/1‑93 (MQ=255)
                                                    tttCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGgtg                                     <  2:1597066/93‑1 (MQ=255)
                                                    tttCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGgtg                                     >  1:7902181/1‑93 (MQ=255)
                                                     tgccAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGgtgt                                    <  2:2940659/91‑1 (MQ=255)
                                                        cAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGtctc                                 >  1:5885652/1‑93 (MQ=255)
                                                             tgctgAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTa                            >  1:5083650/1‑93 (MQ=255)
                                                             tgctgAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTa                            <  1:417026/93‑1 (MQ=255)
                                                             cgctgAGAGATGCTTAATTGTTCCTGGCCGGGGTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTa                            <  2:1555413/92‑1 (MQ=255)
                                                              gctgAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTag                           >  2:7828978/1‑93 (MQ=255)
                                                              gctgAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTag                           >  1:5138870/1‑93 (MQ=255)
                                                                tgAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTagga                         <  2:4283244/93‑1 (MQ=255)
                                                                tgAGAGATGCTTATTTGTACCTGACCGCGCTTACTGTAGAACAGATTAAAAGAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTagga                         <  1:7390039/93‑1 (MQ=255)
                                                                   gagaTGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGc                      >  1:5608023/1‑93 (MQ=255)
                                                                   gagaTGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGc                      >  2:1855681/1‑93 (MQ=255)
                                                                   gagaTGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGc                      >  1:2726826/1‑93 (MQ=255)
                                                                   gagaTGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGc                      >  2:6346248/1‑93 (MQ=255)
                                                                    agaTGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCt                     <  1:7803640/93‑1 (MQ=255)
                                                                     gaTGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCtt                    >  2:2506315/1‑93 (MQ=255)
                                                                     gaTGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCtt                    >  2:6453871/1‑93 (MQ=255)
                                                                       tGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCTTGa                  >  2:935220/1‑93 (MQ=255)
                                                                       agCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCTTGa                  <  1:4084300/92‑1 (MQ=255)
                                                                        gCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCTTGAc                 <  2:6271190/93‑1 (MQ=255)
                                                                        gCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCTTGAc                 >  1:3076713/1‑93 (MQ=255)
                                                                        gCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCTTGAc                 >  2:5079487/1‑93 (MQ=255)
                                                                          ttATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCTTGACgg               >  2:2878889/1‑93 (MQ=255)
                                                                          ttATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCTTGACgg               <  1:1880402/93‑1 (MQ=255)
                                                                             tttGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCTTGACGgtgt            >  1:5452846/1‑93 (MQ=255)
                                                                                 ttCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCTTGACGGTGTAAcc        >  1:3096327/1‑93 (MQ=255)
                                                                                 ttCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCTTGACGGTGTAAcc        >  1:1423482/1‑93 (MQ=255)
                                                                                    cTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCTTGACGGTGTAACCACt     >  2:2726263/1‑93 (MQ=255)
                                                                                       aCCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCTTGACGGTGTAACCACTCAt  <  2:6105298/93‑1 (MQ=255)
                                                                                       aCCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCTTGACGGTGTAACCACTCAt  >  2:663626/1‑93 (MQ=255)
                                                                                          |                                                                                         
TGACGTGTTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCACGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCTTGACGGTGTAACCACTCAT  >  NC_000913/1587846‑1588025

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