Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,483,536 T→C intergenic (‑17/‑250) pepA ← / → yjgP aminopeptidase A/I/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,483,5360TC99.6% 789.4 / ‑5.1 239intergenic (‑17/‑250)pepA/yjgPaminopeptidase A/I/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (128/110);  total (128/111)
Fisher's exact test for biased strand distribution p-value = 4.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

Only 100 of 245 total aligned reads displayed.

GGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAATCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTCCGCTGTTGCGATGACTTCGTGTTAATCTTAAC  >  NC_000913/4483446‑4483624
                                                                                          |                                                                                        
ggTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCt                                                                                        >  2:645811/1‑93 (MQ=255)
 gTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCtt                                                                                       <  1:718191/93‑1 (MQ=255)
  ttCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTa                                                                                      >  1:4320155/1‑93 (MQ=255)
   tCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTaa                                                                                     <  2:2477250/93‑1 (MQ=255)
   tCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTaa                                                                                     <  1:1015498/93‑1 (MQ=255)
      aaGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGa                                                                                  >  1:1176332/1‑93 (MQ=255)
       aGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGaa                                                                                 >  1:3057736/1‑92 (MQ=255)
       aGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGAc                                                                                 <  1:1453492/93‑1 (MQ=255)
          cGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAAc                                                                              >  1:515628/1‑93 (MQ=255)
          cGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAAc                                                                              <  1:3884168/93‑1 (MQ=255)
          cGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAAc                                                                              <  2:112975/93‑1 (MQ=255)
            cccACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAAcgg                                                                            >  2:3586648/1‑93 (MQ=255)
             ccACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAAcggc                                                                           <  2:3242076/93‑1 (MQ=255)
              cACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAAcggcg                                                                          >  1:4272156/1‑93 (MQ=255)
               acgacgATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAAcggcgg                                                                         >  1:427234/1‑93 (MQ=255)
               acgacgATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAAcggcgg                                                                         <  1:3633820/93‑1 (MQ=255)
                cgacgaTGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGt                                                                        <  1:735590/93‑1 (MQ=255)
                  acgaTGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGgggg                                                                      >  2:1758692/1‑93 (MQ=255)
                    gaTGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCt                                                                    >  1:4061700/1‑93 (MQ=255)
                    gaTGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCt                                                                    >  1:3328614/1‑93 (MQ=255)
                      tGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTAc                                                                  >  1:3094192/1‑93 (MQ=255)
                       gCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACa                                                                 <  2:3012330/93‑1 (MQ=255)
                         aGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGa                                                               >  2:3198427/1‑93 (MQ=255)
                         aGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGa                                                               >  2:540838/1‑93 (MQ=255)
                          ggCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGAt                                                              >  2:3304783/1‑93 (MQ=255)
                           gctcTCCACTTTTCCTCCGGGATACCGCCTTCTACACTAAACTCCATGCACTACGCTCCAGAACCTTAAAGACAACGGCGGTGGCTACAGATa                                                             <  2:2448540/90‑1 (MQ=255)
                           gCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATa                                                             <  2:1527381/93‑1 (MQ=255)
                           gCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATa                                                             <  2:3162612/93‑1 (MQ=255)
                            cACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAg                                                            <  1:1858816/93‑1 (MQ=255)
                             aCTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGa                                                           <  2:13182/93‑1 (MQ=255)
                              cTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGaa                                                          >  1:2547107/1‑93 (MQ=255)
                               tCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAAt                                                         >  2:4188766/1‑93 (MQ=255)
                               tCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAAt                                                         >  1:3264381/1‑93 (MQ=255)
                               tCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAAt                                                         <  1:2087303/93‑1 (MQ=255)
                                  gCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGc                                                      <  2:38461/93‑1 (MQ=255)
                                   cTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCa                                                     <  2:2140984/93‑1 (MQ=255)
                                    tGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCaa                                                    <  2:824534/93‑1 (MQ=255)
                                     gTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAg                                                   >  1:1625125/1‑93 (MQ=255)
                                     gTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAg                                                   >  2:16654/1‑93 (MQ=255)
                                     gTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAg                                                   >  2:4318281/1‑93 (MQ=255)
                                      tttCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGc                                                  <  2:1239751/93‑1 (MQ=255)
                                       ttCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCt                                                 <  1:1293789/93‑1 (MQ=255)
                                        tctcCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCtt                                                <  2:4572510/93‑1 (MQ=255)
                                         ctcCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACATTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCttt                                               <  1:285182/93‑1 (MQ=255)
                                          tcCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTc                                              <  1:176737/93‑1 (MQ=255)
                                           ccGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCg                                             <  1:885214/93‑1 (MQ=255)
                                            cGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGt                                            <  1:1077719/93‑1 (MQ=255)
                                            cGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGt                                            >  1:4133753/1‑93 (MQ=255)
                                              ggCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTaa                                          <  2:1839656/93‑1 (MQ=255)
                                              ggCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTaa                                          >  1:468833/1‑93 (MQ=255)
                                              ggCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTaa                                          >  1:3694841/1‑93 (MQ=255)
                                               gCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAAc                                         <  1:321884/93‑1 (MQ=255)
                                                 tACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTc                                       >  1:2679697/1‑93 (MQ=255)
                                                 tACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTc                                       <  1:203869/93‑1 (MQ=255)
                                                  aCCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCa                                      >  1:3441982/1‑93 (MQ=255)
                                                  aCCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCa                                      >  1:3005097/1‑93 (MQ=255)
                                                  aCCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTCCGTAACTCa                                      >  1:4190817/1‑93 (MQ=255)
                                                   ccGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCAt                                     <  1:1530013/93‑1 (MQ=255)
                                                     gCTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGt                                   >  2:4393073/1‑93 (MQ=255)
                                                      cTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTc                                  <  1:3814664/93‑1 (MQ=255)
                                                      cTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTc                                  <  2:1977439/93‑1 (MQ=255)
                                                      cTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTc                                  >  1:1737720/1‑93 (MQ=255)
                                                      cTTTTTACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTc                                  >  2:1323546/1‑93 (MQ=255)
                                                         tttACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTCCGc                               >  2:3090669/1‑93 (MQ=255)
                                                          ttACACTAAACTCCATGCACTACGCTCCTGAATCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTCCGCt                              <  1:3322891/93‑1 (MQ=255)
                                                          ttACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTCCGCt                              >  1:446701/1‑93 (MQ=255)
                                                          ttACACTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTCCGCt                              <  1:3547964/93‑1 (MQ=255)
                                                            acacTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTCCGCTGt                            >  1:923432/1‑93 (MQ=255)
                                                            acacTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTCCGCTGt                            >  2:4329048/1‑93 (MQ=255)
                                                            acacTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTCCGCTGt                            <  1:2232757/93‑1 (MQ=255)
                                                               cTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTCCGCTGTTGc                         <  1:2919542/93‑1 (MQ=255)
                                                               cTAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTCCGCTGTTGc                         <  2:4743845/93‑1 (MQ=255)
                                                                tAAACTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAAATGCAAGCTTTCGTAACTCATGTCCGCTGTTGCg                        >  1:4217577/1‑93 (MQ=255)
                                                                  aaCTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTCCGCTGTTGCGAt                      <  1:2194278/93‑1 (MQ=255)
                                                                  aaCTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTCCGCTGTTGCGAt                      >  2:2009749/1‑93 (MQ=255)
                                                                   aCTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTCCGCTGTTGCGATg                     >  2:3498813/1‑93 (MQ=255)
                                                                   aCTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTCCGCTGTTGCGATg                     >  2:4399726/1‑93 (MQ=255)
                                                                   aCTCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAACTTTCGTAACTCATGTCCGCTGTTGCGATg                     <  2:2630445/93‑1 (MQ=255)
                                                                     tCCATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTCCGCTGTTGCGATGat                   >  2:3442239/1‑92 (MQ=255)
                                                                      ccATGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTCCGCCGTTGCGATGACt                  >  1:49050/1‑93 (MQ=255)
                                                                       cATGCACTACGCTCCTGAACCTTAAAGACAACGGCCGTGGCTACAGATAGAATTGCAAGCTTTCGCAACGAATGTCAGGTGTGGCGATgcagt                 >  2:1199891/1‑89 (MQ=255)
                                                                        aTGCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTCCGCTGTTGCGATGACTTc                >  1:887734/1‑93 (MQ=255)
                                                                          gCACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTCCGCTGTTGCGATGACTTCgt              <  2:1862682/93‑1 (MQ=255)
                                                                           cACTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTCCGCTGTTGCGATGACTTCgtg             >  1:4021940/1‑93 (MQ=255)
                                                                            aCTACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTCCGCTGTTGCGATGACTTCgtgt            >  1:2743352/1‑93 (MQ=255)
                                                                            aCTACGCACCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTCCGCTGTTGCGATGACTTCgtgt            <  2:3350741/93‑1 (MQ=255)
                                                                              tACGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTCCGCTGTTGCGATGACTTCGTGTTa          >  2:3980187/1‑93 (MQ=255)
                                                                                cGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCCAGCTTTCGTAACTCATGTTCGCTGTTGCGATGACTTCGTGTTaaa        >  1:2806212/1‑92 (MQ=255)
                                                                                cGCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTCCGCTGTTGCGATGACTTCGTGTTAAt        <  2:828056/93‑1 (MQ=255)
                                                                                 gCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTCCGCTGTTGCGATGACTTCGTGTTAATc       <  2:1870490/93‑1 (MQ=255)
                                                                                 gCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTCCGCTGTTGCGATGACTTCGTGTTAATc       <  2:3235173/93‑1 (MQ=255)
                                                                                 gCTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCCTAACTCATGTCCGCTGTTGCGATGACTTCGTGTTAATc       <  1:4557597/93‑1 (MQ=255)
                                                                                  cTCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTCCGCTGTTGCGATGACTTCGTGTTAATCt      >  1:2590210/1‑93 (MQ=255)
                                                                                   tCCTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTCCGCTGTTGCGATGACTTCGTGTTAATCtt     >  1:2407492/1‑93 (MQ=255)
                                                                                     cTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTCCGCTGTTGCGATGACTTCGTGTTAATCTTaa   <  1:393843/93‑1 (MQ=255)
                                                                                     cTGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTCCGCTGTTGCGATGACTTCGTGTTAATCTTaa   >  1:3366084/1‑93 (MQ=255)
                                                                                     cTGAACCTGAAAGACAACGGCGCTGGCTACAGATAGAATTTCAAGCTTTCGTAACTCATGTCCGCTGTTGCGATGACTTCGTGTTAATCTTaa   <  1:934449/93‑1 (MQ=255)
                                                                                      tGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTCCGCTGTTGCGATGACTTCGTGTTAATCTTAAc  >  1:4384514/1‑93 (MQ=255)
                                                                                      tGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTCCGCTGTTGCGATGACTTCGTGTTAATCTTAAc  >  1:4158621/1‑93 (MQ=255)
                                                                                      tGAACCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTCCGCTGTTGCGATGACTTCGTGTTAATCTTAAc  >  2:2750693/1‑93 (MQ=255)
                                                                                          |                                                                                        
GGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAATCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAACTCATGTCCGCTGTTGCGATGACTTCGTGTTAATCTTAAC  >  NC_000913/4483446‑4483624

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