Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 988,555 A→G L342P (CTG→CCG)  pncB ← nicotinate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913988,5550AG99.3% 486.0 / ‑5.5 151L342P (CTG→CCG) pncBnicotinate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (74/76);  total (75/76)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 158 total aligned reads displayed.

CTGTCAGAAAGTTTCGCCACCGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCAGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAG  >  NC_000913/988465‑988644
                                                                                          |                                                                                          
cTGTCAGAAAGTTTCGCCACCGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCggg                                                                                          <  2:305636/93‑1 (MQ=255)
 tGTCAGAAAGTTTCGCCACCGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGc                                                                                         <  2:4641263/93‑1 (MQ=255)
 tGTCAGAAAGTTTCGCCACCGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGc                                                                                         >  1:4561587/1‑93 (MQ=255)
 tGTCAGAAAGTTTCGCCACCGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGc                                                                                         <  1:1359673/93‑1 (MQ=255)
      gAAAGTTTCGCCACCGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTGGGGCGAGTc                                                                                   >  1:3900720/1‑93 (MQ=255)
      gAAAGTTTCGCCACCGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTc                                                                                   >  2:2018972/1‑93 (MQ=255)
        aaGTTTCGCCACCGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTccc                                                                                 <  1:1928412/93‑1 (MQ=255)
          gTTTCGCCACCGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCaa                                                                               <  1:15913/93‑1 (MQ=255)
          gTTTCGCCACCGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCaa                                                                               <  1:8207/93‑1 (MQ=255)
            ttCGCCACCGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATa                                                                             <  2:1830491/93‑1 (MQ=255)
                 cACCGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGTTACGCAAGGCGGGGGAAGTCCCAATACCaaa                                                                        >  1:2323802/1‑93 (MQ=255)
                 cACCGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGTGATCGCAGGCCGGGCGAGTCCCAATACCaaa                                                                        >  1:2323804/1‑93 (MQ=255)
                 cACCGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCaaa                                                                        >  2:147570/1‑93 (MQ=255)
                 cACCGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCaaa                                                                        >  1:2001680/1‑93 (MQ=255)
                  cccGGTTTACCGTTACACTCTACTAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCaaaa                                                                       <  2:1185344/92‑1 (MQ=255)
                  aCCGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCaaaa                                                                       >  1:3285878/1‑93 (MQ=255)
                   ccGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATTTCGCAGGTCGGGCGAGTCCCAATACCAAAAc                                                                      >  1:868738/1‑93 (MQ=255)
                   ccGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAAc                                                                      <  1:273738/93‑1 (MQ=255)
                   ccGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAAc                                                                      >  2:3591773/1‑93 (MQ=255)
                    cGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACt                                                                     >  1:3144917/1‑93 (MQ=255)
                     ggTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACtt                                                                    <  1:710179/93‑1 (MQ=255)
                      gTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACGAAAACTTa                                                                   >  1:1784378/1‑93 (MQ=255)
                      gTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTa                                                                   <  1:3882266/93‑1 (MQ=255)
                         tACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAAtt                                                                <  1:3476810/93‑1 (MQ=255)
                         tACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAAtt                                                                <  2:541397/93‑1 (MQ=255)
                           ccGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGc                                                              >  2:4375709/1‑93 (MQ=255)
                           ccGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGc                                                              >  1:3378247/1‑93 (MQ=255)
                            cGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCa                                                             <  1:1307568/93‑1 (MQ=255)
                             gTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCAc                                                            <  2:687293/93‑1 (MQ=255)
                               tACACTCTACCAACTTAATGACAATATTCAGGTGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCAcgc                                                          <  1:3847387/93‑1 (MQ=255)
                               tACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCAcgc                                                          <  2:3739784/93‑1 (MQ=255)
                                acacTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCAcgcg                                                         <  1:2855790/93‑1 (MQ=255)
                                acacTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCAcgcg                                                         >  1:215896/1‑93 (MQ=255)
                                   ctctACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGa                                                      <  1:1451671/93‑1 (MQ=255)
                                   ctctACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGa                                                      <  2:4523468/93‑1 (MQ=255)
                                   ctctACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGa                                                      <  1:2048801/93‑1 (MQ=255)
                                    tctACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGaa                                                     <  1:3234451/93‑1 (MQ=255)
                                    tctACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGaa                                                     <  1:43828/93‑1 (MQ=255)
                                    tctACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGaa                                                     <  2:2157789/93‑1 (MQ=255)
                                     ctACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAag                                                    >  2:2524532/1‑93 (MQ=255)
                                     ctACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAag                                                    <  1:1604746/93‑1 (MQ=255)
                                      tACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAaga                                                   >  1:1435127/1‑93 (MQ=255)
                                      tACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAaga                                                   >  2:780739/1‑93 (MQ=255)
                                        ccAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAgaga                                                 >  1:3006437/1‑93 (MQ=255)
                                         cAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGaa                                                <  2:255518/93‑1 (MQ=255)
                                         cAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGaa                                                <  2:4231735/93‑1 (MQ=255)
                                          aaCTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAg                                               <  2:4004295/93‑1 (MQ=255)
                                          aaCTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAg                                               <  2:4653927/93‑1 (MQ=255)
                                           aCTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGt                                              >  1:4298729/1‑93 (MQ=255)
                                            cTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTg                                             <  2:116093/93‑1 (MQ=255)
                                            cTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTg                                             >  1:2446519/1‑93 (MQ=255)
                                             ttAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTgg                                            <  1:627925/93‑1 (MQ=255)
                                             ttAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTgg                                            >  1:2386462/1‑93 (MQ=255)
                                             ttAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTgg                                            >  1:1089657/1‑93 (MQ=255)
                                             ttAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTgg                                            >  2:198884/1‑93 (MQ=255)
                                              tAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGc                                           <  1:2824230/93‑1 (MQ=255)
                                                aTGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCgg                                         >  2:3941474/1‑93 (MQ=255)
                                                 tGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGt                                        <  2:3795390/93‑1 (MQ=255)
                                                 tGACAATATTCAGGGGGTTTACGTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGt                                        <  2:2222852/93‑1 (MQ=255)
                                                  gACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGta                                       <  1:1256729/93‑1 (MQ=255)
                                                  gACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGta                                       >  1:1741912/1‑93 (MQ=255)
                                                  gACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGta                                       >  1:2234143/1‑93 (MQ=255)
                                                   aCAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGtat                                      <  1:2708032/93‑1 (MQ=255)
                                                   aCAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGtat                                      >  1:4175417/1‑93 (MQ=255)
                                                    cAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGtata                                     <  2:3691833/93‑1 (MQ=255)
                                                    cAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGtata                                     <  2:4579352/93‑1 (MQ=255)
                                                    cAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAATGGCGGtata                                     <  1:3013007/93‑1 (MQ=255)
                                                     aaTATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGtatac                                    >  1:1419912/1‑92 (MQ=255)
                                                     aaTATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAg                                    >  1:570793/1‑93 (MQ=255)
                                                     aaTATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAg                                    >  1:4044278/1‑93 (MQ=255)
                                                      atatTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGc                                   >  2:2366845/1‑93 (MQ=255)
                                                      atatTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGc                                   <  2:2804192/93‑1 (MQ=255)
                                                      atatTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGc                                   >  1:2269025/1‑93 (MQ=255)
                                                      atatTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGc                                   <  2:3550149/93‑1 (MQ=255)
                                                        atTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTc                                 <  1:1689463/93‑1 (MQ=255)
                                                        atTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTc                                 >  1:3984648/1‑93 (MQ=255)
                                                         ttCAGGGGTTTTACCTGGGGGATATCGCAGGTCAGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGTAGAGAAGCAGCGGCATAACTCa                                >  2:4042304/1‑93 (MQ=255)
                                                            aGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAAcc                             <  2:382417/93‑1 (MQ=255)
                                                             ggggTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCg                            >  2:425756/1‑93 (MQ=255)
                                                              gggTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGc                           <  1:467991/93‑1 (MQ=255)
                                                               ggTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCt                          <  1:1245216/93‑1 (MQ=255)
                                                               ggTTTTACCTGGGGGATATCGCAGCTCGTGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCt                          <  1:67938/93‑1 (MQ=255)
                                                                  tttACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACAGCTTTg                       >  1:2386600/1‑93 (MQ=255)
                                                                    tACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCg                     >  2:1382234/1‑93 (MQ=255)
                                                                     aCCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGt                    >  2:2818766/1‑93 (MQ=255)
                                                                       cTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTaa                  >  1:3785898/1‑93 (MQ=255)
                                                                        tGGGGGATATCGCAGGTCGGGCGGGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGGATAACTCAACCGCTTTGCGTaaa                 >  2:2701800/1‑93 (MQ=255)
                                                                        tGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTaaa                 >  2:65435/1‑93 (MQ=255)
                                                                        tGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGGATAGCTCAACCGGTTTGCGTaaa                 >  2:1759549/1‑93 (MQ=255)
                                                                         gggggATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAAt                >  2:3536243/1‑93 (MQ=255)
                                                                            ggATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCa             <  1:4508193/93‑1 (MQ=255)
                                                                            ggATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCa             <  2:1085432/93‑1 (MQ=255)
                                                                              atatCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGa           <  2:3454199/93‑1 (MQ=255)
                                                                              atatCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGa           <  2:2755809/93‑1 (MQ=255)
                                                                               tatCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGAt          >  1:2279967/1‑93 (MQ=255)
                                                                                atCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGAtt         <  1:1119470/93‑1 (MQ=255)
                                                                                    cAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCa     <  2:1836464/93‑1 (MQ=255)
                                                                                      ggTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCaga   >  2:3456734/1‑93 (MQ=255)
                                                                                      ggTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCaga   >  1:363221/1‑93 (MQ=255)
                                                                                      ggTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCaga   <  2:872823/93‑1 (MQ=255)
                                                                                          |                                                                                          
CTGTCAGAAAGTTTCGCCACCGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCAGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAG  >  NC_000913/988465‑988644

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