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,5550AG100.0% 432.9 / NA 137L342P (CTG→CCG) pncBnicotinate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (73/64);  total (73/64)

Only 100 of 142 total aligned reads displayed.

CTGTCAGAAAGTTTCGCCACCGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCAGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAG  >  NC_000913/988465‑988644
                                                                                          |                                                                                         
cTGTCAGAAAGTTTCGCCACCGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCggg                                                                                         >  1:3909362/1‑93 (MQ=255)
  gTCAGAAAGTTTCGCCACCGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCg                                                                                       >  1:3491761/1‑93 (MQ=255)
   tCAGAAAGTTTCGCCACCGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGCTTTACCTGGGGGATATCGCAGGTCGGGCGa                                                                                      >  1:3060258/1‑93 (MQ=255)
     aGAAAGTTTCGCCACCGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAgc                                                                                    >  1:1832299/1‑92 (MQ=255)
        aaGTTTCGCCACCGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTccc                                                                                 >  1:4300379/1‑93 (MQ=255)
          gTTTCGCCACCGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCaa                                                                               <  1:405995/93‑1 (MQ=255)
           tttCGCCACCGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAAt                                                                              <  2:1641974/93‑1 (MQ=255)
               gCCACCGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCa                                                                          >  1:3113052/1‑93 (MQ=255)
               gCCACCGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCa                                                                          >  1:2979169/1‑93 (MQ=255)
                ccACCGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCaa                                                                         >  2:1081394/1‑93 (MQ=255)
                  aCCGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCaaaa                                                                       >  1:3485610/1‑93 (MQ=255)
                  aCCGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGGCCCAATACCaaaa                                                                       >  1:510447/1‑93 (MQ=255)
                   ccGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAAc                                                                      >  2:1135709/1‑93 (MQ=255)
                   ccGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAAc                                                                      <  1:3459795/93‑1 (MQ=255)
                    cGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACt                                                                     <  1:1252965/93‑1 (MQ=255)
                     ggTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACtt                                                                    <  1:3470947/93‑1 (MQ=255)
                      gTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTa                                                                   >  1:4477727/1‑93 (MQ=255)
                      gTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTa                                                                   >  1:1408447/1‑93 (MQ=255)
                       tttACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTaa                                                                  <  2:1742852/93‑1 (MQ=255)
                        ttACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAAt                                                                 >  2:1259940/1‑93 (MQ=255)
                        ttACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAAt                                                                 <  1:2010573/93‑1 (MQ=255)
                        ttACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAAt                                                                 <  2:4634950/93‑1 (MQ=255)
                          aCCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTg                                                               >  1:128262/1‑93 (MQ=255)
                           ccGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGc                                                              >  2:3778029/1‑93 (MQ=255)
                             gTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCAc                                                            <  1:4073054/93‑1 (MQ=255)
                             gTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCAc                                                            <  1:4368973/93‑1 (MQ=255)
                              ttACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCAcg                                                           >  1:775791/1‑93 (MQ=255)
                              ttACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCAcg                                                           >  2:2223828/1‑93 (MQ=255)
                                acacTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCAcgcg                                                         >  2:3506363/1‑93 (MQ=255)
                                 ccctctACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCgg                                                        <  2:2928230/91‑1 (MQ=255)
                                 cacTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCgg                                                        <  1:4255287/93‑1 (MQ=255)
                                 cacTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCgg                                                        <  1:4279203/93‑1 (MQ=255)
                                  acTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCggg                                                       <  1:1051090/93‑1 (MQ=255)
                                     ctACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAag                                                    <  1:1377012/93‑1 (MQ=255)
                                      tACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAaga                                                   <  1:3653254/93‑1 (MQ=255)
                                       aCCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAagag                                                  <  1:285319/93‑1 (MQ=255)
                                         cAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGaa                                                <  2:1329314/93‑1 (MQ=255)
                                         cAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGaa                                                >  2:2030214/1‑93 (MQ=255)
                                         cAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGaa                                                <  1:1685452/93‑1 (MQ=255)
                                         cAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGaa                                                <  2:139717/93‑1 (MQ=255)
                                           aCTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGt                                              <  1:2032142/93‑1 (MQ=255)
                                           aCTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACAAAAACTTAATTGCACGCGGGAAGAGAAGt                                              <  1:3036553/93‑1 (MQ=255)
                                            cTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTg                                             >  2:3189688/1‑93 (MQ=255)
                                             ttAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTgg                                            >  1:2531517/1‑93 (MQ=255)
                                             ttAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTgg                                            >  2:4261605/1‑93 (MQ=255)
                                             ttAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTgg                                            >  1:1645454/1‑93 (MQ=255)
                                              tAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGc                                           >  2:2710663/1‑93 (MQ=255)
                                              tAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGc                                           <  2:2897198/93‑1 (MQ=255)
                                              tAATGACAATATTCAGGGGTTATACCTGGGGGATATCGCAGGTCGGGCTAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGc                                           <  2:1136738/93‑1 (MQ=255)
                                               aaTGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCg                                          >  2:2295698/1‑93 (MQ=255)
                                               aaTGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCg                                          >  2:253080/1‑93 (MQ=255)
                                                aTGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCgg                                         >  2:107408/1‑93 (MQ=255)
                                                 tGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGt                                        >  2:2071959/1‑93 (MQ=255)
                                                 tGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGt                                        >  1:2147971/1‑93 (MQ=255)
                                                 tGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGt                                        <  1:13056/93‑1 (MQ=255)
                                                  gACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAACACTTAATTGCACGCGGGAAGAGAAGGGGCggga                                       >  1:1681620/1‑91 (MQ=255)
                                                   aCAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGtat                                      >  1:768470/1‑93 (MQ=255)
                                                    cAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGtata                                     <  2:4513157/93‑1 (MQ=255)
                                                     aaTATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAg                                    >  2:2954850/1‑93 (MQ=255)
                                                      atatTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGc                                   <  2:871475/93‑1 (MQ=255)
                                                        atTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAGCTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTc                                 <  1:1345537/93‑1 (MQ=255)
                                                        atTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTc                                 <  2:322932/93‑1 (MQ=255)
                                                            aGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAAcc                             >  2:562060/1‑93 (MQ=255)
                                                              gggTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGc                           <  2:3979492/93‑1 (MQ=255)
                                                                gTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCtt                         >  2:734269/1‑93 (MQ=255)
                                                                gTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCtt                         >  1:4370710/1‑93 (MQ=255)
                                                                gTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCtt                         >  1:3214876/1‑93 (MQ=255)
                                                                 ttttACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCttt                        >  1:1652008/1‑93 (MQ=255)
                                                                 ttttACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCttt                        >  2:311946/1‑93 (MQ=255)
                                                                   ttACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGc                      <  2:3891725/93‑1 (MQ=255)
                                                                     aCCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGt                    >  2:1279769/1‑93 (MQ=255)
                                                                      ccTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTa                   >  1:2275575/1‑93 (MQ=255)
                                                                       cTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTaa                  >  2:946472/1‑93 (MQ=255)
                                                                       cTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTaa                  <  1:3605913/93‑1 (MQ=255)
                                                                       cTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTaa                  >  2:2248617/1‑93 (MQ=255)
                                                                        tGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTaaa                 >  1:2378255/1‑93 (MQ=255)
                                                                        tGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTaaa                 >  2:2561798/1‑93 (MQ=255)
                                                                        tGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTaaa                 <  2:4457822/93‑1 (MQ=255)
                                                                        tGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTaaa                 >  1:357451/1‑93 (MQ=255)
                                                                        tGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTaaa                 <  2:336902/93‑1 (MQ=255)
                                                                          ggggATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATc               <  2:3973844/93‑1 (MQ=255)
                                                                          ggggATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATc               >  2:115802/1‑93 (MQ=255)
                                                                          ggggATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATc               <  2:2275462/93‑1 (MQ=255)
                                                                           gggATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATcc              <  2:1848361/93‑1 (MQ=255)
                                                                            ggATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCa             <  1:1911395/93‑1 (MQ=255)
                                                                             gATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAg            >  1:1012514/1‑93 (MQ=255)
                                                                              atatCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGa           <  2:4163845/93‑1 (MQ=255)
                                                                              atatCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGa           >  1:2294274/1‑93 (MQ=255)
                                                                              atatCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGa           <  1:3158838/93‑1 (MQ=255)
                                                                               tatCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGAt          <  2:4033018/93‑1 (MQ=255)
                                                                                atCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGAtt         <  1:3323172/93‑1 (MQ=255)
                                                                                atCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGAtt         >  2:2995031/1‑93 (MQ=255)
                                                                                 tCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTg        >  1:2555871/1‑93 (MQ=255)
                                                                                 tCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACAGCTTTGCGTAAATTCAGATTg        >  2:3195616/1‑93 (MQ=255)
                                                                                    cAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCa     <  1:802868/93‑1 (MQ=255)
                                                                                     aGGTCGGGCGAGTCCCAATTCCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCag    >  2:4577079/1‑93 (MQ=255)
                                                                                     aGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCag    <  1:4013264/93‑1 (MQ=255)
                                                                                      ggTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCaga   >  1:3150631/1‑93 (MQ=255)
                                                                                      ggTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCaga   >  2:74886/1‑93 (MQ=255)
                                                                                       gTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCagag  >  1:347306/1‑93 (MQ=255)
                                                                                          |                                                                                         
CTGTCAGAAAGTTTCGCCACCGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCAGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAG  >  NC_000913/988465‑988644

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