Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,823,338 G→C 7.0% intergenic (+53/‑177) nadE → / → cho NAD synthetase, NH3/glutamine‑dependent/endonuclease of nucleotide excision repair

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,823,3380GC 7.0% 248.4 / 9.4 71intergenic (+53/‑177)nadE/choNAD synthetase, NH3/glutamine‑dependent/endonuclease of nucleotide excision repair
Reads supporting (aligned to +/- strand):  ref base G (33/33);  new base C (3/2);  total (36/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.42e-01

GCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACT  >  NC_000913/1823199‑1823485
                                                                                                                                           |                                                                                                                                                   
gCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGCATCCGGCGTGaa                                                                                                                                            <  3:95233/149‑1 (MQ=255)
    gAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGc                                                                                                                                         >  1:683416/1‑148 (MQ=255)
       cAATAGAGAACTGGTATCTGAAATACGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCtt                                                                                                                                      <  1:384551/148‑1 (MQ=255)
        aaTAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTa                                                                                                                                     <  2:831807/148‑1 (MQ=255)
         aTAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTAt                                                                                                                                    >  2:540306/1‑148 (MQ=255)
          tAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTc                                                                                                                                  <  3:416381/149‑1 (MQ=255)
           agagAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGa                                                                                                                                             <  1:97833/137‑1 (MQ=255)
           agagAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGa                                                                                                                                             >  2:97833/1‑137 (MQ=255)
             agaACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGa                                                                                                                                <  2:417558/148‑1 (MQ=255)
              gaACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGAcc                                                                                                                              >  3:292591/1‑149 (MQ=255)
              gaACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGAcc                                                                                                                              >  3:461119/1‑149 (MQ=255)
              gaACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGAcc                                                                                                                              >  3:284545/1‑149 (MQ=255)
                   ggTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCtata                                                                                                                          >  1:590721/1‑148 (MQ=255)
                         tGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTAt                                                                                                                   >  3:665569/1‑149 (MQ=255)
                          gAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTAt                                                                                                                   >  2:153898/1‑148 (MQ=255)
                          gAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATg                                                                                                                  <  3:187780/149‑1 (MQ=255)
                                 gAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTc                                                                                                           >  3:143116/1‑149 (MQ=255)
                                 gAACATAAACGCCGTCCGCCAATTACCGTTCTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGCATCCGGCGTGa                                                                                                                                             <  3:590140/115‑1 (MQ=255)
                                 gAACATAAACGCCGTCCGCCAATTACCGTTCTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGCATCCGGCGTGa                                                                                                                                             >  4:590140/1‑115 (MQ=255)
                                  aaCATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCa                                                                                                          >  3:484308/1‑149 (MQ=255)
                                   aCATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCaa                                                                                                         <  4:292591/149‑1 (MQ=255)
                                   aCATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCa                                                                                                          <  1:540306/148‑1 (MQ=255)
                                       aaaCGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATAc                                                                                                      >  2:58006/1‑148 (MQ=255)
                                       aaaCGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACa                                                                                                     >  3:124041/1‑149 (MQ=255)
                                          cGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACAtt                                                                                                   >  2:350186/1‑148 (MQ=255)
                                          cGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTg                                                                                                  <  4:124041/149‑1 (MQ=255)
                                            ccGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATgcggcg                                                                                                                                                >  2:551112/1‑101 (MQ=255)
                                            ccGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATgcggcg                                                                                                                                                <  1:551112/101‑1 (MQ=255)
                                               tCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAgg                                                                                              <  1:898566/148‑1 (MQ=255)
                                                 cGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGt                                                                                           <  3:37724/149‑1 (MQ=255)
                                                  gCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGt                                                                                           >  1:206008/1‑148 (MQ=255)
                                                      aTTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGaa                                                                                       >  2:719281/1‑148 (MQ=255)
                                                         aCCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTAt                                                                                                                                    >  4:839175/1‑100 (MQ=255)
                                                         aCCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTAt                                                                                                                                    <  3:839175/100‑1 (MQ=255)
                                                                 cGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCAAATAGGCCTGACa                                                                            <  2:634231/148‑1 (MQ=255)
                                                                  gatgatTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAg                                                                           <  2:683416/148‑1 (MQ=255)
                                                                   atgatTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAgg                                                                          >  2:117732/1‑148 (MQ=255)
                                                                      atTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGt                                                                       <  1:153898/148‑1 (MQ=255)
                                                                              aaaaaGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCa                                                              >  4:215298/1‑149 (MQ=255)
                                                                                aaaGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAg                                                             <  2:430384/148‑1 (MQ=255)
                                                                                aaaGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTCACAGGCGTAGCACGTCAGa                                                            <  3:531373/149‑1 (MQ=255)
                                                                                   gTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACgg                                                         >  3:895672/1‑149 (MQ=255)
                                                                                             gCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCttt                                               >  4:466869/1‑149 (MQ=255)
                                                                                             gCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCtt                                                <  1:109566/148‑1 (MQ=255)
                                                                                             gCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCtt                                                <  1:704786/148‑1 (MQ=255)
                                                                                              cacaTATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTg                                              <  4:484308/149‑1 (MQ=255)
                                                                                                 atatTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTc                                            >  1:737355/1‑148 (MQ=255)
                                                                                                 atatTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCa                                           >  4:15045/1‑149 (MQ=255)
                                                                                                        aTTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTa                                                                      <  4:660176/115‑1 (MQ=255)
                                                                                                        aTTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTa                                                                      >  3:660176/1‑115 (MQ=255)
                                                                                                         ttGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGAccc                                    >  1:876274/1‑148 (MQ=255)
                                                                                                                aaaaaGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCtt                             >  1:139704/1‑148 (MQ=255)
                                                                                                                 aaaaGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCtttt                           <  3:827577/149‑1 (MQ=255)
                                                                                                                  aaaGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCtttt                           <  2:139704/148‑1 (MQ=255)
                                                                                                                   aaGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAg                                                                                               <  3:631149/79‑1 (MQ=255)
                                                                                                                   aaGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAg                                                                                               >  4:631149/1‑79 (MQ=255)
                                                                                                                     gagTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTaa                       <  3:215298/149‑1 (MQ=255)
                                                                                                                      agTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAAt                      <  4:895672/149‑1 (MQ=255)
                                                                                                                        tAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAAt                                                                                                        >  2:79028/1‑65 (MQ=255)
                                                                                                                        tAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAAt                                                                                                        <  1:79028/65‑1 (MQ=255)
                                                                                                                         aaTTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCg                    <  2:206008/148‑1 (MQ=255)
                                                                                                                          aTTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCt                  <  4:143116/149‑1 (MQ=255)
                                                                                                                          aTTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCt                  >  3:401033/1‑149 (MQ=255)
                                                                                                                           tttGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCtt                 >  3:614240/1‑149 (MQ=255)
                                                                                                                           tttGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCt                  >  2:877140/1‑148 (MQ=255)
                                                                                                                               ttCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACAtcg                                                                                                  <  4:619242/64‑3 (MQ=255)
                                                                                                                               ttCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTg                                                                                                  >  3:619242/1‑64 (MQ=255)
                                                                                                                               ttCACGCCGCATCCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTccc              >  1:438418/1‑148 (MQ=255)
                                                                                                                                tCACGCCGCATCCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCg             >  2:158692/1‑148 (MQ=255)
                                                                                                                                 cACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGcc           <  4:401033/149‑1 (MQ=255)
                                                                                                                                 cACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTAATCGCTTCCCGcc           <  2:737355/148‑1 (MQ=255)
                                                                                                                                   cGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACg                                                          <  2:162846/100‑1 (MQ=255)
                                                                                                                                   cGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACg                                                          >  1:162846/1‑100 (MQ=255)
                                                                                                                                           gcggcgTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACt  <  1:719281/148‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                   
GCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACT  >  NC_000913/1823199‑1823485

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: