Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,597,846 (T)5→4 intergenic (‑286/‑138) livK ← / → panM leucine transporter subunit/PanD autocleavage accelerator, panothenate synthesis

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,597,8420T.99.3% 649.9 / ‑2.7 145intergenic (‑282/‑142)livK/panMleucine transporter subunit/PanD autocleavage accelerator, panothenate synthesis
Reads supporting (aligned to +/- strand):  ref base T (0/1);  major base . (81/62);  minor base A (0/1);  total (81/64)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 156 total aligned reads displayed.

TTATTGATAGATAAACAGAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATAAAAAAGAGAAAGCCTCCGATTAAAT  >  NC_000913/3597744‑3597937
                                                                                                  |                                                                                               
 tATTGATAGATAAACAGAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACAttt                                                                                               >  1:1163407/1‑100 (MQ=255)
  aTTGATAGATAAACAGAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACAtttt                                                                                              <  2:3158557/100‑1 (MQ=255)
  aTTGATAGATAAACAGAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACAtttt                                                                                              >  2:465314/1‑100 (MQ=255)
  aTTGATAGATAAACAGAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACAtttt                                                                                              <  1:1751709/100‑1 (MQ=255)
  aTTGATAGATAAACAGAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACAtttt                                                                                              >  2:2530220/1‑100 (MQ=255)
    tGATAGATAAACAGAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAAAAATTCACTATCATTCAGGGGATTATGCTGGACAttttca                                                                                            >  1:2442392/1‑98 (MQ=255)
     gatagataAACAGAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACAttttcat                                                                                           >  2:1194862/1‑97 (MQ=255)
     gatagataAACAGAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACAttttcat                                                                                           >  1:270459/1‑97 (MQ=255)
       tagataAACAGAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATtc                                                                                        <  2:1780651/100‑1 (MQ=255)
       tagataAACAGAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATtc                                                                                        <  2:2542931/100‑1 (MQ=255)
        agataAACAGAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCCGAATAATTCGCTATCGTTCAGGGGATTATGCTGGACAttttccttct                                                                                        >  1:2341555/1‑94 (MQ=255)
         gataAACAGAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATtctc                                                                                      >  2:856788/1‑100 (MQ=255)
         gataAACAGAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATtctc                                                                                      >  2:3262640/1‑100 (MQ=255)
         gataAACAGAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATtctc                                                                                      >  1:6249/1‑100 (MQ=255)
         gataAACAGAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATtctc                                                                                      >  2:2182869/1‑100 (MQ=255)
         gataAACAGAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATtcac                                                                                      >  1:2113652/1‑98 (MQ=255)
           taacccgcaaaaaGTGCCTTTGTCAGCATAAAATAAAGGAACAAAGGGCGGAATAATTCACTATCATTAAAAGGAATATACTGGACATTTTCATTCTCTa                                                                                    >  1:1575416/9‑100 (MQ=255)
           taAACAGAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTa                                                                                    >  1:2586456/1‑100 (MQ=255)
           taAACAGAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTa                                                                                    <  2:2260194/100‑1 (MQ=255)
             aaCAGAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAAt                                                                                  <  2:1757714/100‑1 (MQ=255)
               cAGAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGt                                                                                >  2:2849758/1‑100 (MQ=255)
                aGAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGtt                                                                               >  2:109882/1‑100 (MQ=255)
                aGAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGtt                                                                               >  2:1970316/1‑100 (MQ=255)
                 gAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGttt                                                                              >  1:1630686/1‑100 (MQ=255)
                 gAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGttt                                                                              >  2:2592829/1‑100 (MQ=255)
                   aaaaaGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTa                                                                            >  2:1512225/1‑100 (MQ=255)
                      aaGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAAtt                                                                         >  1:1888447/1‑100 (MQ=255)
                        gTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAAtttt                                                                       >  2:1850021/1‑100 (MQ=255)
                         tGCCTTTGTCGGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTg                                                                      >  2:1007549/1‑100 (MQ=255)
                         tGCCTTTGTCATCATAAAATAACGGCACAAAGGGCGGAATTATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTg                                                                      <  1:2727714/100‑1 (MQ=255)
                         tGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTg                                                                      <  2:1148753/100‑1 (MQ=255)
                         tGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTg                                                                      <  1:776520/100‑1 (MQ=255)
                         tGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTg                                                                      >  1:379634/1‑100 (MQ=255)
                         tGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTg                                                                      <  2:2403884/100‑1 (MQ=255)
                         tGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCAAATGTTTTAATTTTg                                                                      <  1:2417733/100‑1 (MQ=255)
                          gCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGt                                                                     <  1:3193852/100‑1 (MQ=255)
                          gCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGt                                                                     >  2:2214948/1‑100 (MQ=255)
                          gCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGt                                                                     >  1:1177060/1‑100 (MQ=255)
                               tGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAatta                                                                >  1:3040082/1‑100 (MQ=255)
                                 tCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAattatt                                                              <  2:852595/100‑1 (MQ=255)
                                  cAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTg                                                             >  1:3217381/1‑100 (MQ=255)
                                  cAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTg                                                             >  2:1816208/1‑100 (MQ=255)
                                   aGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGc                                                            <  2:540680/100‑1 (MQ=255)
                                   aGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGc                                                            <  1:2881272/100‑1 (MQ=255)
                                   aGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGc                                                            <  2:1879953/100‑1 (MQ=255)
                                    gCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCt                                                           <  1:1291528/100‑1 (MQ=255)
                                    gCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCt                                                           >  2:3169100/1‑100 (MQ=255)
                                     cATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTg                                                          <  2:2885378/100‑1 (MQ=255)
                                       tAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGtt                                                        >  2:3433851/1‑100 (MQ=255)
                                          aaTAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTaaa                                                     <  2:2530385/100‑1 (MQ=255)
                                            tAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTaaaaa                                                   <  2:297770/100‑1 (MQ=255)
                                             aaCGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTaaaaaa                                                  <  1:604551/100‑1 (MQ=255)
                                              aCGGCACAAAGGGCGGAATAATTCACTATCATTCAGTGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAAt                                                 <  1:2056802/100‑1 (MQ=255)
                                              aCGGCACAAAGGGCGGAATAATTCACTATCATTCAGGTGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAAt                                                 <  1:1551487/100‑1 (MQ=255)
                                              aCGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAAt                                                 >  2:3400217/1‑100 (MQ=255)
                                              aCGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGAACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAAt                                                 <  1:1633366/100‑1 (MQ=255)
                                                ggCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAAtt                                                >  1:2863795/1‑99 (MQ=255)
                                                 gCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTaa                                              >  1:860151/1‑100 (MQ=255)
                                                 gCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTa                                               <  1:74504/99‑1 (MQ=255)
                                                  cacaAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAAt                                             >  1:1136363/1‑100 (MQ=255)
                                                  cacaAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAAt                                             >  1:668165/1‑100 (MQ=255)
                                                   acaAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATc                                            <  2:972517/100‑1 (MQ=255)
                                                    caAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCa                                           >  2:345685/1‑100 (MQ=255)
                                                      aaGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCAcc                                         <  1:1194847/100‑1 (MQ=255)
                                                      aaGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCAcc                                         <  2:1040470/100‑1 (MQ=255)
                                                      aaGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCAcc                                         >  1:300571/1‑100 (MQ=255)
                                                      aaGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCAcc                                         <  1:2592794/100‑1 (MQ=255)
                                                           cGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCa                                    >  2:3225277/1‑100 (MQ=255)
                                                             gAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCaaa                                  >  2:597529/1‑100 (MQ=255)
                                                              ttccgatctcTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCaaaa                                 <  1:3252133/91‑1 (MQ=255)
                                                               ataatTCAATATCAGGCAGGAGATTATGATGGAAAGTTTCAGATGCTCATGATTTAATTCTGTAATTATTGCTGTAAAAAAATTAATAACATGACaaaa                                 <  1:2950219/99‑1 (MQ=255)
                                                                     cACTATCATTCAGGGGATTATGCTGGACATTTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATa                           <  2:715046/100‑1 (MQ=255)
                                                                     cACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATaa                          >  1:1704501/1‑100 (MQ=255)
                                                                     cACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATaa                          <  1:1740211/100‑1 (MQ=255)
                                                                      aCTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATaaa                         >  1:439896/1‑100 (MQ=255)
                                                                       cTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATaaaa                        <  2:1428559/100‑1 (MQ=255)
                                                                        tATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATaaaaa                       <  1:2172265/100‑1 (MQ=255)
                                                                        tATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATaaaaa                       >  2:627014/1‑100 (MQ=255)
                                                                         ataaTTCAGGAGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATATTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATaaaaaa                      <  1:2714718/97‑1 (MQ=255)
                                                                         aTCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAATAAATTAATCACTTGCCAAAATAAATaaaa                        >  1:2254798/1‑98 (MQ=255)
                                                                         aTCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATaaaaaa                      >  2:2668642/1‑100 (MQ=255)
                                                                          tCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATAAAAAag                     <  1:744875/100‑1 (MQ=255)
                                                                          tCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATAAAAAag                     <  2:3362148/100‑1 (MQ=255)
                                                                            aTTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATAAAAAagag                   >  2:2831359/1‑100 (MQ=255)
                                                                            aTTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATAAAAAagag                   <  1:1107246/100‑1 (MQ=255)
                                                                            aTTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATAAAAAagag                   >  2:1176727/1‑100 (MQ=255)
                                                                            aTTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATAAAAAagag                   >  1:1486695/1‑100 (MQ=255)
                                                                             ttCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATAAAAAAgaga                  >  2:1676524/1‑100 (MQ=255)
                                                                              tCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATAAAAAAGAGaa                 <  2:1971844/100‑1 (MQ=255)
                                                                               cAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATAAAAAAGAGaaa                >  1:3351889/1‑100 (MQ=255)
                                                                                aGGGGATTATGCTGGACATTTTCATTCTCTAGTGTTTTAATTATGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATAAAAAAGAGAAAg               >  2:2349039/1‑100 (MQ=255)
                                                                                     aTTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATAAAAAAGAGAAAGCCTct          >  2:739488/1‑99 (MQ=255)
                                                                                     aTTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATAAAAAAGAGAAAGCCTcc          <  1:2685352/100‑1 (MQ=255)
                                                                                        aTGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATAAAAAAGAGAAAGCCTCCGAt       >  1:1240204/1‑100 (MQ=255)
                                                                                         tGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATAAAAAAGAGAAAGCCTCCGAtt      >  1:1396657/1‑100 (MQ=255)
                                                                                         tGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATAAAAAAGAGAAAGCCTCCGAtt      <  2:66626/100‑1 (MQ=255)
                                                                                          gCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATAAAAAAGAGAAAGCCTCCGATTa     >  1:782427/1‑100 (MQ=255)
                                                                                          gCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATAAAAAAGAGAAAGCCTCCGATTa     >  1:2481099/1‑100 (MQ=255)
                                                                                            tGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCATAAGAAATAAAAAAGAGAAAGCCTCCGATTaaa   <  2:3300600/100‑1 (MQ=255)
                                                                                             ggACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATAAAAAAGAGAAAGCCTCCGATTAAat  <  2:1879222/100‑1 (MQ=255)
                                                                                                  |                                                                                               
TTATTGATAGATAAACAGAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATAAAAAAGAGAAAGCCTCCGATTAAAT  >  NC_000913/3597744‑3597937

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