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.6% 742.4 / ‑6.0 237L342P (CTG→CCG) pncBnicotinate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (117/119);  minor base T (0/1);  total (117/120)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 248 total aligned reads displayed.

CTGTCAGAAAGTTTCGCCACCGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCAGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAG  >  NC_000913/988465‑988644
                                                                                          |                                                                                         
cTGTCAGAAAGTTTCGCCACCGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCggg                                                                                         >  2:6479809/1‑93 (MQ=255)
cTGTCAGAAAGTTTCGCCACCGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCggg                                                                                         <  1:4579139/93‑1 (MQ=255)
 tGTCAGAAAGTTTCGCCACCGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGc                                                                                        <  2:4977619/93‑1 (MQ=255)
 tGTCAGAAAGTTTCGCCACCGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGc                                                                                        >  2:2905409/1‑93 (MQ=255)
 tGTCAGAAAGTTTCGCCACCGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGc                                                                                        >  1:774586/1‑93 (MQ=255)
   tCAGAAAGTTTCGCCACCGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGa                                                                                      >  2:5213169/1‑93 (MQ=255)
   tCAGAAAGTTTCGCCACCGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGa                                                                                      <  1:5208801/93‑1 (MQ=255)
     aGAAAGTTTCGCCACCGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAgc                                                                                    >  2:1338988/1‑92 (MQ=255)
     aGAAAGTTTCGCCACCGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGt                                                                                    <  2:7675759/93‑1 (MQ=255)
      gAAAGTTTCGCCACCGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTc                                                                                   >  1:1297297/1‑93 (MQ=255)
      gAAAGTTTCGCCACCGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTc                                                                                   <  1:6505359/93‑1 (MQ=255)
          gTTTCGCCACCGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCaa                                                                               <  2:2430023/93‑1 (MQ=255)
           tttCGCCACCGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAAt                                                                              <  1:1795342/93‑1 (MQ=255)
           tttCGCCACCGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAAt                                                                              <  2:189452/93‑1 (MQ=255)
            ttCGCCACCGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATa                                                                             <  2:453169/93‑1 (MQ=255)
             tCGCCACCGGTTTACCGTTACACTCTACCACCTTAATGCCAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATAc                                                                            <  2:6970552/93‑1 (MQ=255)
               gCCACCGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCa                                                                          >  1:5425043/1‑93 (MQ=255)
                 cACCGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCaaa                                                                        >  1:278326/1‑93 (MQ=255)
                 cACCGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCaaa                                                                        >  2:1149406/1‑93 (MQ=255)
                  aCCGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCaaaa                                                                       >  1:681163/1‑93 (MQ=255)
                  aCCGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCaaaa                                                                       >  1:2520891/1‑93 (MQ=255)
                   ccGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAAc                                                                      <  2:425995/93‑1 (MQ=255)
                    cGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACt                                                                     <  2:6905136/93‑1 (MQ=255)
                     ggTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACtt                                                                    >  2:7534121/1‑93 (MQ=255)
                      gTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTa                                                                   >  1:6096231/1‑93 (MQ=255)
                      gTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTa                                                                   >  1:3479599/1‑93 (MQ=255)
                       tttACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTaa                                                                  >  2:2764560/1‑93 (MQ=255)
                       tttACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTaa                                                                  >  1:1573338/1‑93 (MQ=255)
                       tttACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTaa                                                                  <  1:1965615/93‑1 (MQ=255)
                        ttACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAAt                                                                 <  2:5750437/93‑1 (MQ=255)
                        ttACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAAt                                                                 >  1:5724970/1‑93 (MQ=255)
                            cGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCa                                                             <  2:1988703/93‑1 (MQ=255)
                             gTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCAc                                                            <  2:4855329/93‑1 (MQ=255)
                             gTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCAc                                                            <  2:1240193/93‑1 (MQ=255)
                               tACACTCTACCAACTTAATGACAATCTTCAGGGGTTTTACCTGGGGGATATCGAAGGTAGGGCGAGTCACAATACCAATACTTAACTGCAcgc                                                          >  2:7419247/1‑93 (MQ=255)
                               tACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCAcgc                                                          <  1:1417778/93‑1 (MQ=255)
                                    tctACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGaa                                                     <  2:5462602/93‑1 (MQ=255)
                                    tctACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGaa                                                     <  2:6404865/93‑1 (MQ=255)
                                     ctACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAag                                                    <  1:7141348/93‑1 (MQ=255)
                                      tACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAaga                                                   >  2:1528570/1‑93 (MQ=255)
                                      tACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAaga                                                   <  2:4046482/93‑1 (MQ=255)
                                      tACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAaga                                                   <  2:4312098/93‑1 (MQ=255)
                                       aCCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAagag                                                  >  2:4201444/1‑93 (MQ=255)
                                        ccAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAgaga                                                 <  2:3841753/93‑1 (MQ=255)
                                        ccAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAgaga                                                 <  1:3904250/93‑1 (MQ=255)
                                         cAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGaa                                                <  1:4710884/93‑1 (MQ=255)
                                         cAACTTAATGACAATATTCAGGGGTTTTACCTGGGAGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGaa                                                <  2:4604350/93‑1 (MQ=255)
                                          aaCTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAg                                               >  1:552840/1‑93 (MQ=255)
                                          aaCTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAg                                               <  1:6263350/93‑1 (MQ=255)
                                           aCTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGt                                              <  1:7666592/93‑1 (MQ=255)
                                           aCTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGt                                              <  2:3920992/93‑1 (MQ=255)
                                              tAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGc                                           <  2:7278000/93‑1 (MQ=255)
                                              tAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGc                                           >  1:6330300/1‑93 (MQ=255)
                                              tAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGc                                           >  2:1890049/1‑93 (MQ=255)
                                              tAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGc                                           <  1:3627183/93‑1 (MQ=255)
                                               aaTGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCACAACTTAATTGCACGCGGGAAGAGAAGGGGCg                                          >  1:4416299/1‑93 (MQ=255)
                                               aaTGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCg                                          >  1:4416297/1‑93 (MQ=255)
                                               aaTGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCg                                          <  1:5493887/93‑1 (MQ=255)
                                               aaTGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCg                                          >  1:5056866/1‑93 (MQ=255)
                                                   aCAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGtaa                                      >  1:6572619/1‑92 (MQ=255)
                                                    cAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGtata                                     <  1:1055926/93‑1 (MQ=255)
                                                        atTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGcgc                                 <  1:293692/93‑3 (MQ=255)
                                                         ttCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCa                                >  1:6320285/1‑93 (MQ=255)
                                                          tCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCaa                               >  2:7154302/1‑93 (MQ=255)
                                                            aGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAAcc                             <  2:5686745/93‑1 (MQ=255)
                                                            aGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAAcc                             >  1:4608029/1‑93 (MQ=255)
                                                             ggggTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCg                            >  1:906621/1‑93 (MQ=255)
                                                             ggggTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCg                            <  2:7908810/93‑1 (MQ=255)
                                                             ggggTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCg                            >  1:624115/1‑93 (MQ=255)
                                                              gggTTTTACCTGGGGGATATCGCAGGTCTGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGc                           <  2:3225632/93‑1 (MQ=255)
                                                                 ttttACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCttt                        >  1:177877/1‑93 (MQ=255)
                                                                 ttttACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCttt                        <  2:2827667/93‑1 (MQ=255)
                                                                  tttACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTg                       <  1:6390392/93‑1 (MQ=255)
                                                                  tttACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTg                       <  2:4394230/93‑1 (MQ=255)
                                                                     aCCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGt                    >  1:580526/1‑93 (MQ=255)
                                                                      ccTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTa                   >  1:3201165/1‑93 (MQ=255)
                                                                         gggggATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTGAAt                >  2:5512293/1‑93 (MQ=255)
                                                                         gggggATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAAt                >  2:1159157/1‑93 (MQ=255)
                                                                         gggggATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAAACGCTTTGCGTAAAt                >  2:4471502/1‑93 (MQ=255)
                                                                         gggggATATCGCAGGGCGGGCGAGTCCCAATACCAAAAGTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAAt                <  1:4164667/93‑1 (MQ=255)
                                                                          ggggATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATc               >  1:6000363/1‑93 (MQ=255)
                                                                          ggggATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATc               <  1:6662225/93‑1 (MQ=255)
                                                                           gggATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATcc              <  2:1635199/93‑1 (MQ=255)
                                                                            ggATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCa             >  1:2063672/1‑93 (MQ=255)
                                                                             gATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAg            >  2:6648017/1‑93 (MQ=255)
                                                                             gATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAg            >  2:2708607/1‑93 (MQ=255)
                                                                             gATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAg            >  2:3118118/1‑93 (MQ=255)
                                                                              atatCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGa           <  1:2380149/93‑1 (MQ=255)
                                                                               tatCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGAt          <  1:3074193/93‑1 (MQ=255)
                                                                               tatCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGAt          >  1:7429393/1‑93 (MQ=255)
                                                                                atCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGAtt         >  2:7796048/1‑93 (MQ=255)
                                                                                atCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGAtt         <  1:2729803/93‑1 (MQ=255)
                                                                                 tCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTg        <  1:1501347/93‑1 (MQ=255)
                                                                                  cGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGt       <  1:4872565/93‑1 (MQ=255)
                                                                                   gCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTc      >  2:3681210/1‑93 (MQ=255)
                                                                                   gCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTc      >  1:6418797/1‑93 (MQ=255)
                                                                                    cAGGTCGGACGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCa     <  1:7804410/93‑1 (MQ=255)
                                                                                     aGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCag    <  1:3730904/93‑1 (MQ=255)
                                                                                      ggTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCaga   >  1:7219868/1‑93 (MQ=255)
                                                                                      ggTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCaga   >  2:5332406/1‑93 (MQ=255)
                                                                                          |                                                                                         
CTGTCAGAAAGTTTCGCCACCGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCAGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAG  >  NC_000913/988465‑988644

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