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,9360AG98.4% 606.9 / ‑1.3 191R30R (CGT→CGCb1505putative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base A (3/0);  new base G (91/96);  total (94/97)
Fisher's exact test for biased strand distribution p-value = 1.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

Only 100 of 201 total aligned reads displayed.

GTGACGTGTTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCACGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCTTGACGGTGTAACCACTCATT  >  NC_000913/1587845‑1588026
                                                                                           |                                                                                          
 tGACGTGTTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCTTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACcgcg                                                                                          >  1:3387354/1‑93 (MQ=255)
 tGACGTGTTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACcgcg                                                                                          >  2:4097556/1‑93 (MQ=255)
 tGACGTGTTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACcgcg                                                                                          >  1:2482901/1‑93 (MQ=255)
  gACGTGTTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCgcgc                                                                                         >  1:2736467/1‑93 (MQ=255)
   aCGTGTTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCt                                                                                        >  1:3198046/1‑93 (MQ=255)
   aCGTGTTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCt                                                                                        >  2:403369/1‑93 (MQ=255)
   aCGTGTTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCACGCt                                                                                        >  1:3797266/1‑93 (MQ=255)
     gtgtTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTa                                                                                      >  2:3318300/1‑93 (MQ=255)
       gtTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACt                                                                                    >  1:1363868/1‑93 (MQ=255)
       gtTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACt                                                                                    <  1:1114689/93‑1 (MQ=255)
       gtTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACt                                                                                    <  1:2039355/93‑1 (MQ=255)
       gtTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACt                                                                                    <  1:233963/93‑1 (MQ=255)
          ccAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTa                                                                                 <  1:1125129/93‑1 (MQ=255)
          ccAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTa                                                                                 <  2:2200554/93‑1 (MQ=255)
          ccAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTa                                                                                 <  2:1017022/93‑1 (MQ=255)
            aGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGATCGCGCTTACTGTAGa                                                                               >  1:1152257/1‑93 (MQ=255)
            aGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGa                                                                               <  1:4386224/93‑1 (MQ=255)
            aGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGa                                                                               <  1:4677979/93‑1 (MQ=255)
            aGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGa                                                                               <  2:1153060/93‑1 (MQ=255)
            aGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGa                                                                               <  1:592991/93‑1 (MQ=255)
             gTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGaa                                                                              <  2:3065770/93‑1 (MQ=255)
             gTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGaa                                                                              >  1:1461043/1‑93 (MQ=255)
               aaCTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACa                                                                            <  2:4480948/93‑1 (MQ=255)
                aCTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAg                                                                           <  1:3410992/93‑1 (MQ=255)
                 cTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGa                                                                          <  1:4415581/93‑1 (MQ=255)
                 cTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGa                                                                          >  2:1557824/1‑93 (MQ=255)
                 cTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGa                                                                          <  1:159059/93‑1 (MQ=255)
                 cTTTGGCGACTGACACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGa                                                                          <  2:1369807/93‑1 (MQ=255)
                   ttGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGAtt                                                                        <  1:2790715/93‑1 (MQ=255)
                    tGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTa                                                                       >  2:3401404/1‑93 (MQ=255)
                     ggCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTaa                                                                      <  2:3140421/93‑1 (MQ=255)
                        gACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAaat                                                                   >  2:3520009/1‑93 (MQ=255)
                          cTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACCGATTACaataa                                                                 >  1:137404/1‑93 (MQ=255)
                           tGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAaataat                                                                <  2:3658352/93‑1 (MQ=255)
                            ggCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATc                                                               >  2:2338723/1‑93 (MQ=255)
                             gCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCa                                                              <  2:4288944/93‑1 (MQ=255)
                             gCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCa                                                              <  1:3685618/93‑1 (MQ=255)
                               aCTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAAt                                                            <  2:4374996/93‑1 (MQ=255)
                                 tGAAAAATGTCGTACCGTAATTTCCAAGCTTCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATaa                                                          <  1:3351954/93‑1 (MQ=255)
                                 tGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATaa                                                          >  1:2289971/1‑93 (MQ=255)
                                 tGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATaa                                                          <  1:1526790/93‑1 (MQ=255)
                                 tGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATaa                                                          >  1:882624/1‑93 (MQ=255)
                                  gAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATaaa                                                         >  2:1732047/1‑93 (MQ=255)
                                      aaTGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATtgt                                                     <  1:2855369/93‑1 (MQ=255)
                                       aTGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATtgtg                                                    >  2:1885484/1‑93 (MQ=255)
                                       aTGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATtgtg                                                    <  2:1282648/93‑1 (MQ=255)
                                       aTGTCGGTCCGGAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATtgtg                                                    <  2:4595017/93‑1 (MQ=255)
                                         gTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTc                                                  <  2:2534020/93‑1 (MQ=255)
                                         gTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTc                                                  >  2:676983/1‑93 (MQ=255)
                                          tCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCt                                                 <  2:723247/93‑1 (MQ=255)
                                           cGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTg                                                <  2:639150/93‑1 (MQ=255)
                                            gTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGc                                               >  1:846011/1‑93 (MQ=255)
                                            gTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGc                                               <  1:2482479/93‑1 (MQ=255)
                                            gTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGc                                               >  1:3102312/1‑93 (MQ=255)
                                              tccGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTc                                             <  2:3602601/92‑1 (MQ=255)
                                              aCCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTc                                             >  1:4343854/1‑93 (MQ=255)
                                              aCCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTc                                             >  1:2816098/1‑93 (MQ=255)
                                               ccGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCa                                            <  1:527613/93‑1 (MQ=255)
                                                  tAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTg                                         >  1:346993/1‑93 (MQ=255)
                                                  tAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTg                                         <  1:3550090/93‑1 (MQ=255)
                                                   aaTTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTgg                                        >  1:2505625/1‑93 (MQ=255)
                                                   aaTTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTgg                                        <  1:329611/93‑1 (MQ=255)
                                                    aTTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGgt                                       <  1:305089/93‑1 (MQ=255)
                                                    aTTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGgt                                       <  2:2981371/93‑1 (MQ=255)
                                                    aTTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGgt                                       <  2:3835499/93‑1 (MQ=255)
                                                    aTTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGAGTCTGCTCATTGgt                                       >  2:2415703/1‑93 (MQ=255)
                                                       tCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGtc                                    <  2:1714175/93‑1 (MQ=255)
                                                         cAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGtctc                                  >  2:2745228/1‑93 (MQ=255)
                                                          aaGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTcc                                 >  1:4573587/1‑93 (MQ=255)
                                                          aaGCTGCTGAGAGATGCTTATTTGTTACTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTcc                                 >  2:2168438/1‑93 (MQ=255)
                                                            gctgctGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCgg                               <  2:3149979/93‑1 (MQ=255)
                                                            gctgctGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCgg                               >  1:2104048/1‑93 (MQ=255)
                                                              tgctgAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTa                             <  2:1299662/93‑1 (MQ=255)
                                                              tgctgAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTa                             >  2:1236226/1‑93 (MQ=255)
                                                               gctgAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTag                            >  2:279929/1‑93 (MQ=255)
                                                                 tgAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTagga                          <  1:1907621/93‑1 (MQ=255)
                                                                  gagagaTGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTaggag                         <  1:714625/93‑1 (MQ=255)
                                                                    gagaTGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGc                       <  2:1182996/93‑1 (MQ=255)
                                                                    gagaTGCTTATTTGTTCCTGACCGCGCTTACTGTACAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGc                       >  2:4695333/1‑93 (MQ=255)
                                                                      gaTGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCtt                     <  2:4746689/93‑1 (MQ=255)
                                                                      gaTGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCtt                     >  1:184512/1‑93 (MQ=255)
                                                                      gaTGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCtt                     <  1:2565199/93‑1 (MQ=255)
                                                                      gaTGCTTATTCGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCtg                     >  2:3598621/1‑92 (MQ=255)
                                                                       aTGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCTTg                    <  2:2065602/93‑1 (MQ=255)
                                                                        tGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCTTGa                   <  1:3337318/93‑1 (MQ=255)
                                                                          cTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCTTGACg                 >  2:4613516/1‑93 (MQ=255)
                                                                           ttATTTGTGCCTGACCGCGCTTACTGTAGAGCAGATTAAAATAGTCAATAAATTGTGTCTGCGCATTGGGGTCTCCGGTAGGAGGCTTGACgg                <  1:264245/93‑1 (MQ=255)
                                                                             ttttGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCTTGACGgtg              <  2:1609142/92‑1 (MQ=255)
                                                                             aTTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCTTGACGgtg              >  2:3921650/1‑93 (MQ=255)
                                                                                 gTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCTTGACGGTGTAAc          >  2:26782/1‑93 (MQ=255)
                                                                                   tCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCTTGACGGTGTAACCa        >  2:3706783/1‑93 (MQ=255)
                                                                                    ccTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGAAGGAGGCTTGACGGTGTAACCAc       >  1:3105186/1‑93 (MQ=255)
                                                                                     cTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCTTAACGGTGTAACCACt      <  2:1330345/93‑1 (MQ=255)
                                                                                      tGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCTTGACGGTGTAACCACTc     >  1:4298251/1‑93 (MQ=255)
                                                                                      tGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCTTGACGGTGTAACCACTc     >  1:2124435/1‑93 (MQ=255)
                                                                                      tGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCTTGACGGTGTAACCACTc     >  1:4146259/1‑93 (MQ=255)
                                                                                      tGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCTTGACGGTGTAACCACTc     <  1:3598000/93‑1 (MQ=255)
                                                                                        aCCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCTTGCGGGTGAAACCACTctt   >  2:914235/1‑91 (MQ=255)
                                                                                        aCCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCTTGACGGTGTAACCACTCAt   >  1:586533/1‑93 (MQ=255)
                                                                                         ccACGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCTTGACGGTGTAACCACTCAtt  <  1:4104552/93‑1 (MQ=255)
                                                                                           |                                                                                          
GTGACGTGTTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCACGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCTTGACGGTGTAACCACTCATT  >  NC_000913/1587845‑1588026

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