Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,712,011 T→C 6.0% L531P (CTT→CCT)  nadB → quinolinate synthase, L‑aspartate oxidase (B protein) subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,712,0110TC 6.0% 205.1 / 3.2 66L531P (CTT→CCT) nadBquinolinate synthase, L‑aspartate oxidase (B protein) subunit
Reads supporting (aligned to +/- strand):  ref base T (30/32);  new base C (2/2);  total (32/34)
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.09e-01

GCGTCTCAAATAATTTGCTGGAGCTGCGTAATCTGGTACAGGTTGCCGAGTTGATTGTTCGCTGTGCAATGATGCGTAAAGAGAGTCGGGGGTTGCATTTCACGCTGGATTATCCGGAACTGCTCACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGA  >  NC_000913/2711864‑2712156
                                                                                                                                                   |                                                                                                                                                 
gcgTCTCAAATAATTTGCTGGAGCTGCGTAATCTGGTACAGGTTGCCGAGTTGATTGTTCGCTGTGCAATGATGCGTAAAGAGAGTCGGGGGTTGCATTTCACGCTGGATTATCCGGAACTGCTCACCCATTCCGGTCCGTCGATCCtt                                                                                                                                                  >  4:225862/1‑149 (MQ=255)
gcgTCTCAAATAATTTGCTGGAGCTGCGTAATCTGGTACAGGTTGCCGAGTTGAGTGTTCGCTGTGCAATGATGCGTAAAGAGAGTCGGGGGTTGCATTTCACGCTGGATTATCCGGAACTGCTCACCCATTCCGGTCCGTCGATCCtt                                                                                                                                                  >  1:71031/1‑149 (MQ=255)
 cgTCTCAAATAATTTGCTGGAGCTGCGTAATCTGGTACAGGTTGCCGAGTTGATTGTTCGCTGTGCAATGATGCGTAAAGAGAGTCGGGGGTTGCATTTCACGCTGGATTATCCGGAACTGCTCACCCATTCCGGTCCGTCGATCCttt                                                                                                                                                 <  3:225862/149‑1 (MQ=255)
                gCTGGAGCTGCGTAATCTGGTACAGGTTGCCGAGTTGATTGTTCGCTGTGCAATGATGCGTAAAGAGAGTCGGGGGTTGCATTTCACGCTGGATTATCCGGAACTGCTCACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCAt                                                                                                                                   <  3:110527/148‑1 (MQ=255)
                 cTGGAGCTGCGTAATCTGGTACAGGTTGCCGAGTTGATTGTTCGCTGTGCAATGATGCGTAAAGAGAGTCGGGGGTTGCATTTCACGCTGGATTATCCGGAACTGCTCACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCAtt                                                                                                                                  >  7:306438/1‑148 (MQ=255)
                          cGTAATCTGGTACAGGTTGCCGAGTTGATTGTTCGCTGTGCAATGATGCGTAAAGAGAGTCGGGGGTTGCATTTCACGCTGGATTATCCGGAACTGCTCACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACa                                                                                                                         <  8:277860/148‑1 (MQ=255)
                          cGTAATCTGGTACAGGTTGCCGAGTTGATTGTGCGCTGTGCAATGATGCGTAAAGAGAGTCGGGGGTTGCATTTCACGCTGGATTATCCGGAACTGCTCACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAAc                                                                                                                          <  6:157445/147‑1 (MQ=255)
                            tAATCTGGTACAGGTTGCCGAGTTGATTGTTCGCTGTGCAATGATGCGTAAAGAGAGTCGGGGGTTGCATTTCACGCTGGATTATCCGGAACTGCTCACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGAt                                                                                                                      <  8:2775/149‑1 (MQ=255)
                                       aGGTTGCCGAGTTGATTGTTCGCTGTGCAATGATGCGTAAAGAGAGTCGGGGGTTGCATTTCACGCTGGATTATCCGGAACTGCTCACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTgg                                                                                                           <  4:250711/149‑1 (MQ=255)
                                              cGAGTTGATTGTTCGCTGTGCAATGATGCGTAAAGAGAGTCGGGGGTTGCATTTCACGCTGGATTATCCGGAACTGCTCACCCATTCCGGTCCGTCGATCCTTTCCCCCg                                                                                                                                           <  1:216578/110‑1 (MQ=255)
                                              cGAGTTGATTGTTCGCTGTGCAATGATGCGTAAAGAGAGTCGGGGGTTGCATTTCACGCTGGATTATCCGGAACTGCTCACCCATTCCGGTCCGTCGATCCTTTCCCCCg                                                                                                                                           >  2:216578/1‑110 (MQ=255)
                                                aGTTGATTGTTCGCTGTGCAATGATGCGTAAAGAGAGTCGGGGGTTGCATTTCACGCTGGATTATCCGGAACTGCTCACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGcc                                                                                                  <  8:261840/149‑1 (MQ=255)
                                                 gTTGATTGTTCGCTGTGCAATGATGCGTAAAGAGAGTCGGGGGTTGCATTTCACGCTGGATTATCCGGAACTGCTCACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCg                                                                                                 <  3:265296/149‑1 (MQ=255)
                                                  ttGATTGTTCGCTGTGCAATGATGCGTAAAGAGAGTCGGGGGTTGCATTTCACGCTGGATTATCCGGAACTGCTCACCCATTCCGGTCCGTCGATCCCTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCg                                                                                                 >  2:3928/1‑148 (MQ=255)
                                                    gATTGTTCGCTGTGCAATGATGCGTAAAGAGAGTCGGGGGTTGCATTTCACGCTGGATTATCCGGAACTGCTCACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGcc                                                                                                              >  7:34165/1‑133 (MQ=255)
                                                    gATTGTTCGCTGTGCAATGATGCGTAAAGAGAGTCGGGGGTTGCATTTCACGCTGGATTATCCGGAACTGCTCACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGcc                                                                                                              <  8:34165/133‑1 (MQ=255)
                                                           cGCTGTGCAATGATGCGTAAAGAGAGTCGGGGGTTGCATTTCACGCTGGATTATCCGGAACTGCTCACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCtt                                                                                        <  1:316362/148‑1 (MQ=255)
                                                            gCTGTGCAATGATGCGTAAAGAGAGTCGGGGGTTGCATTTCACGCTGGATTATCCGGAACTTCTCACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTAc                                                                                                                                <  3:75220/107‑1 (MQ=255)
                                                            gCTGTGCAATGATGCGTAAAGAGAGTCGGGGGTTGCATTTCACGCTGGATTATCCGGAACTACTCACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTAc                                                                                                                                >  4:75220/1‑107 (MQ=255)
                                                             cTGTGCAATGATGCGTAAAGAGAGTCGGGGGTTGCATTTCACGCTGGATTATCCGGAACTGCTCACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAGACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCgg                                                                                     <  2:214292/149‑1 (MQ=255)
                                                                 gCAATGATGCGTAAAGAGAGTCGGGGGTTGCATTTCACGCTGGATTATCCGGAACTGCTCACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCAtt                                                                                                                                  <  6:342152/100‑1 (MQ=255)
                                                                 gCAATGATGCGTAAAGAGAGTCGGGGGTTGCATTTCACGCTGGATTATCCGGAACTGCTCACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCAtt                                                                                                                                  >  5:342152/1‑100 (MQ=255)
                                                                   aaTGATGCGTAAAGAGAGTCGGGGGTTGCATTTCACGCTGGATTATCCGGAACTGCTCACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGt                                                                               >  4:162730/1‑149 (MQ=255)
                                                                       atgCGTAAAGAGAGTCGGGGGTTGCATTTCACGCTGGATTATCCGGAACTGCTCACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGtata                                                                                             <  4:326785/131‑1 (MQ=255)
                                                                       atgCGTAAAGAGAGTCGGGGGTTGCATTTCACGCTGGATTATCCGGAACTGCTCACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGtata                                                                                             >  3:326785/1‑131 (MQ=255)
                                                                       atgCGTAAAGAGAGTCGGGGGTTGCATTTCACGCTGGATTATCCGGAACTGCTCACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTtct                                                                            <  3:311971/148‑1 (MQ=255)
                                                                               agagagTCGGGGGTTGCATTTCACGCTGGATTATCCGGAACTGCTCACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGcc                                                                   >  3:181214/1‑149 (MQ=255)
                                                                                gagagTCGGGGGTTGCATTTCACGCTGGATTATCCGGAACTGCTCACCCATTCCGGTCCGGCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGTATAGTTCTGGTCTGGccc                                                                  <  4:4449/149‑1 (MQ=255)
                                                                                        ggggTTGCATTTCACGCTGGATTATCCGGAACTGCTCACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGa                                                          >  8:300842/1‑149 (MQ=255)
                                                                                              gCATTTCACGCTGGATTATCCGGAACTGCTCACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGGCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAg                                                     <  8:67741/148‑1 (MQ=255)
                                                                                                aTTTCACGCTGGATTATCCGGAACTGCTCACCCATTCCGGTCCGTCGATCCTTTCCCCCgg                                                                                                                                          <  3:303441/61‑1 (MQ=255)
                                                                                                aTTTCACGCTGGATTATCCGGAACTGCTCACCCATTCCGGTCCGTCGATCCTTTCCCCCgg                                                                                                                                          >  4:303441/1‑61 (MQ=255)
                                                                                                 tttCACGCTGGATTATCCGGAACTGCTCCCCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGAt                                                 >  1:15843/1‑149 (MQ=255)
                                                                                                 tttCACGCTGGATTATCCGGAACTGCTCACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTaa                                                      >  2:343444/1‑144 (MQ=255)
                                                                                                 tttCACGCTGGATTATCCGGAACTGCTCACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGAt                                                 >  4:110741/1‑149 (MQ=255)
                                                                                                 tttCACGCTGGATTATCCGGAACTGCTCACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGAt                                                 >  7:355509/1‑149 (MQ=255)
                                                                                                 attCACGCTGGATTCTCCGGAACTGCTCACCCATTCCGGTCCGTCGCTCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCGGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGAt                                                 <  7:283429/148‑1 (MQ=255)
                                                                                                  ttCACGCTGGATTATCCGGAACTGCTCACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAAt                                                            >  5:308640/1‑137 (MQ=255)
                                                                                                  ttCACGCTGGATTATCCGGAACTGCTCACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAAt                                                            <  6:308640/137‑1 (MQ=255)
                                                                                                    cACGCTGGATTATCCGGGACTGCTCACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGc                                              >  6:249337/1‑149 (MQ=255)
                                                                                                     aCGCTGGATTATCCGGAACTGCTCACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAACGCCTGGGTCAGCGCCGTATACGCTTCGGAAAAGTTCTGGTCTGGCCCACGACTGACTAAGCGATCGCt                                             <  2:71031/149‑1 (MQ=255)
                                                                                                      cGCTGGATTATCCGGAACTGCTCACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTa                                            <  1:343444/149‑1 (MQ=255)
                                                                                                        cTGGATTATCCGGAACTGCTCACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTaaa                                          <  7:258688/149‑1 (MQ=255)
                                                                                                         tGGATTATCCGGAACTGCTCACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAg                                         >  1:77065/1‑149 (MQ=255)
                                                                                                                tCCGGAACTGCTCACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATtctc                                  <  3:110741/149‑1 (MQ=255)
                                                                                                                tCCGGAACTGCTCACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATAAGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATtctc                                  >  8:96022/1‑149 (MQ=255)
                                                                                                                  cGGAACTGCTCACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTccc                                <  4:181214/149‑1 (MQ=255)
                                                                                                                   ggAACTGCTCACCCATTCCGGTCCGTCGATCGTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCg                               >  2:107940/1‑149 (MQ=255)
                                                                                                                    gAACTGCTCAACCATTCCGGTCCGGCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGc                              <  2:302847/149‑1 (MQ=255)
                                                                                                                        tGCTCACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGc                          >  4:92350/1‑149 (MQ=255)
                                                                                                                           tCACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGACGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCggg                       >  5:238069/1‑149 (MQ=255)
                                                                                                                            cACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCggg                       >  3:131165/1‑148 (MQ=255)
                                                                                                                             aCCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCggg                       >  4:331208/1‑147 (MQ=255)
                                                                                                                             aCCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGga                     <  8:355509/149‑1 (MQ=255)
                                                                                                                             aCCCATTCCGGTCCGTCGATCCCTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGga                     <  1:3928/149‑1 (MQ=255)
                                                                                                                               ccATTCCGGTCCGTCGATCCCTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGccc                                                                  >  7:361438/1‑102 (MQ=255)
                                                                                                                               ccATTCCGGTCCGTCGATCCCTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGccc                                                                  <  8:361438/102‑1 (MQ=255)
                                                                                                                                      ggTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCcag            <  2:15843/149‑1 (MQ=255)
                                                                                                                                       gTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCcagc           <  4:131165/149‑1 (MQ=255)
                                                                                                                                       gTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCCGCAGGCGGGAGAATGcttta           <  8:306438/149‑5 (MQ=255)
                                                                                                                                            tCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATtc                                    >  8:180839/1‑119 (MQ=255)
                                                                                                                                            tCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATtc                                    <  7:180839/119‑1 (MQ=255)
                                                                                                                                            tCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCAc      >  5:277303/1‑149 (MQ=255)
                                                                                                                                              gATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCAccc    >  6:278451/1‑149 (MQ=255)
                                                                                                                                               aTCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATa                                                                                 <  5:9335/71‑1 (MQ=255)
                                                                                                                                               aTCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATa                                                                                 >  6:9335/1‑71 (MQ=255)
                                                                                                                                               aTCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCAccc    >  8:174635/1‑148 (MQ=255)
                                                                                                                                               aTCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCAccc    >  4:335055/1‑148 (MQ=255)
                                                                                                                                                tCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGa  >  3:367015/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GCGTCTCAAATAATTTGCTGGAGCTGCGTAATCTGGTACAGGTTGCCGAGTTGATTGTTCGCTGTGCAATGATGCGTAAAGAGAGTCGGGGGTTGCATTTCACGCTGGATTATCCGGAACTGCTCACCCATTCCGGTCCGTCGATCCTTTCCCCCGGCAATCATTACATAAACAGATAAAAAGCCTGGGTCAGCGCCGTATACGCTTCGGAATAGTTCTGGTCTGGCCCACGAATGACTAAGCGATCGCTAAAGCATTCTCCCGCCTGCGGGGAGAATGCCAGCAGCACCCGA  >  NC_000913/2711864‑2712156

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