Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,627,567 C→T Q85* (CAG→TAG)  nadR → trifunctional protein: nicotinamide mononucleotide adenylyltransferase, ribosylnicotinamide kinase, transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,627,5670CT97.3% 283.0 / ‑4.1 73Q85* (CAG→TAG) nadRtrifunctional protein: nicotinamide mononucleotide adenylyltransferase, ribosylnicotinamide kinase, transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base T (32/39);  total (32/41)
Fisher's exact test for biased strand distribution p-value = 5.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACGATCGGTGTCGTATTCGGTAAGTTCTACCCACTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGCAGCCGACCGTGCCGGATCGTCTGCGTTGGTTATTGCAAACTTTTAAATAT  >  NC_000913/4627421‑4627713
                                                                                                                                                  |                                                                                                                                                  
tGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACGATCGGTGTCGTATTCGGTAAGTTCTACCCACTGCATACCGGACATATCTACCTTATCTAgcg                                                                                                                                                >  2:563351/1‑151 (MQ=255)
 gAATGCCAAGATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACGATCGGTGTCGTATTCGGTAAGTTCTACCCACTGCATACCGGACATATCTACCTTATCTAgcgc                                                                                                                                               <  1:150125/151‑1 (MQ=255)
  aaTGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACGATCGGTGTCGTATTCGGTAAGTTCTACCCACTGCATACCGGACATATCTACCTTATCTAgcgcg                                                                                                                                              <  2:133597/151‑1 (MQ=255)
  aaTGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACGATCGGTGTCGTATTCGGTAAGTTCTACCCACTGCATACCGGACATATCTACCTTATCTAgcgcg                                                                                                                                              >  1:828279/1‑151 (MQ=255)
      ccAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACGATCGGTGTCGTATTCGGTAAGTTCTACCCACTGCATACCGGACATATCTACCTTATCTAGCGCGCCTg                                                                                                                                          <  2:242079/151‑1 (MQ=255)
       cAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACGATCGGTGTCGTATTCGGTAAGTTCTACCCACTGCATACCGGACATATCTACCTTATCTAGCGCGCCTGt                                                                                                                                         <  2:989053/151‑1 (MQ=255)
        aaaaTCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACGATCGGTGTCGTATTCGGTAAGTTCTACCCACTGCATACCGGACATATCTACCTTATCTAGCGCGCCTGt                                                                                                                                         >  1:107066/1‑150 (MQ=255)
        aaaaTCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACGATCGGTGTCGTATTCGGTAAGTTCTACCCACTGCATACCGGACATATCTACCTTATCTAGCGCGCCTGTa                                                                                                                                        <  2:93197/151‑1 (MQ=255)
        aaaaTCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACGATCGGTGTCGTATTCGGTAAGTTCTACCCACTGCATACCGGACATATCTACCTTATCTAGCGCGCCTGTa                                                                                                                                        <  2:225651/151‑1 (MQ=255)
          aaTCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACGATCGGTGTCGTATTCGGTAAGTTCTACCCACTGCATACCGGACATATCTACCTTATCTAGCGCGCCTGTAGc                                                                                                                                      >  1:1062572/1‑151 (MQ=255)
          aaTCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACGATCGGTGTCGTATTCGGTAAGTTCTACCCACTGCATACCGGACATATCTACCTTATCTAGCGCGCCTGTAGc                                                                                                                                      <  2:180828/151‑1 (MQ=255)
           aTCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACGATCGGTGTCGTATTCGGTAAGTTCTACCCACTGCATACCGGACATATCTACCTTATCTAGCGCGCCTGTAGcc                                                                                                                                     <  2:1368479/151‑1 (MQ=255)
             cAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACGATCGGTGTCGTATTCGGTAAGTTCTACCCACTGCATACCGGACATATCTACCTTATCTAGCGCGCCTGTAGcc                                                                                                                                     >  2:272188/1‑149 (MQ=255)
              aaaaGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACGATCGGTGTCGTATTCGGTAAGTTCTACCCACTGCATACCGGACATATCTACCTTATCTAGCGCGCCTGTAGCCAgg                                                                                                                                  >  2:1361886/1‑151 (MQ=255)
              aaaaGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACGATCGGTGTCGTATTCGGTAAGTTCTACCCACTGCATACCGGACATATCTACCTTATCTAGCGCGCCTGTAGCCAgg                                                                                                                                  >  2:658853/1‑151 (MQ=255)
              aaaaGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACGATCGGTGTCGTATTCGGTAAGTTCTACCCACTGCATACCGGACATATCTACCTTATCTAGCGCGCCTGTAGCCAgg                                                                                                                                  <  2:141075/151‑1 (MQ=255)
                     ccAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACGATCGGTGTCGTATTCGGTAAGTTCTACCCACTGCATACCGGACATATCTACCTTATCTAGCGCGCCTGTAGCCAGGTTGACg                                                                                                                            >  2:1235294/1‑150 (MQ=255)
                         cgcgcAAAAGCTGGGGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACGATCGGTGTCGTATTCGGTAAGTTCTACCCACTGCATACCGGACATATCTACCTTATCTAGCGCGCCTGTAGCCAGGTTGACGAGCTg                                                                                                                       <  2:43037/151‑1 (MQ=255)
                         cgcgcAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACGATCGGTGTCGTATTCGGTAAGTTCTACCCACTGCATACCGGACATATCTACCTTATCTAGCGCGCCTGTAGCCAGGTTGACGAGCTg                                                                                                                       >  1:530744/1‑151 (MQ=255)
                          gcgcAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACGATCGGTGTCGTATTCGGTAAGTTCTACCCACTGCATACCGGACATATCTACCTTATCTAGCGCGCCTGTAGCCAGGTTGACGAGCTGc                                                                                                                      >  2:473036/1‑151 (MQ=255)
                                    tGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACGATCGGTGTCGTATTCGGTAAGTTCTACCCACTGCATACCGGACATATCTACCTTATCTAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTAt                                                                                                            >  1:238134/1‑151 (MQ=255)
                                    tGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACGATCGGTGTCGTATTCGGTAAGTTCTACCCACTGCATACCGGACATATCTACCTTATCTAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTAt                                                                                                            >  1:453142/1‑151 (MQ=255)
                                     ggAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACGATCGGTGTCGTATTCGGTAAGTTCTACCCACTGCATACCGGACATATCTACCTTATCTAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTAt                                                                                                            <  1:95611/150‑1 (MQ=255)
                                     ggAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACGATCGGTGTCGTATTCGGTAAGTTCTACCCACTGCATACCGGACATATCTACCTTATCTAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTAt                                                                                                            <  1:1134700/150‑1 (MQ=255)
                                          cGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACGATCGGTGTCGTATTCGGTAAGTTCTACCCACTGCATACCGGACATATCTACCTTATCTAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGttt                                                                                                      >  2:1050826/1‑151 (MQ=255)
                                           gTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACGATCGGTGTCGTATTCGGTAAGTTCTACCCACTGCATACCGGACATATCTACCTTATCTAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGtttt                                                                                                     <  2:268919/151‑1 (MQ=255)
                                              gCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACGATCGGTGTCGTATTCGGTAAGTTCTACCCACTGCATACCGGACATATCTACCTTATCTAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGa                                                                                                   <  2:1131503/150‑1 (MQ=255)
                                                      ttttGGGGCTTGAGTTTCCCCGGCAGAAGAAAACGATCGGTGTCGTATTCGGTAAGTTCTACCCACTGCATACCGGACATATCTACCTTATCTAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCg                                                                                          >  2:856262/1‑151 (MQ=255)
                                                           gggCTTGAGTTTCCCCGGCAGAAGAAAACGATCGGTGTCGTATTCGGTAAGTTCTACCCACTGCATACCGGACATATCTACCTTATCTAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGAcc                                                                                     >  1:804642/1‑151 (MQ=255)
                                                              cTTGAGTTTCCCCGGCAGAAGAAAACGATCGGTGTCGTATTCGGTAAGTTCTACCCACTGCATACCGGACATATCTACCTTATCTAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACcgcg                                                                                  >  1:174543/1‑151 (MQ=255)
                                                                 gAGTTTCCCCGGCAGAAGAAAACGATCGGTGTCGTATTCGGTAAGTTCTACCCACTGCATACCGGACATATCTACCTTATCTAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACcgcgcg                                                                                <  1:403675/150‑1 (MQ=255)
                                                                   gTTTCCCCGGCAGAAGAAAACGATCGGTGTCGTATTCGGTAAGTTCTACCCACTGCATACCGGACATATCTACCTTATCTAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCgttg                                                                             >  1:760255/1‑151 (MQ=255)
                                                                      tCCCCGGCAGAAGAAAACGATCGGTGTCGTATTCGGTAAGTTCTACCCACTGCATACCGGACATATCTACCTTATCTAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCgttgtt                                                                           <  1:294364/150‑1 (MQ=255)
                                                                               gaagaaAACGAACGGTGTCGTATTCGGTAAGTTCTACCCACTGCATACCGGACATATCTACCTTATCTAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACa                                                                   >  1:1020971/1‑149 (MQ=255)
                                                                                aagaaAACGATCGGTGTCGTATTCGGTAAGTTCTACCCACTGCATACCGGACATATCTACCTTATCTAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTg                                                                <  2:78766/151‑1 (MQ=255)
                                                                                  gaaAACGATCGGTGTCGTATTCGGTAAGTTCTACCCACTGCATACCGGACATATCTACCTTATCTAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGcc                                                              >  1:619804/1‑151 (MQ=255)
                                                                                        gATCGGTGTCGTATTCGGTAAGTTCTACCCACTGCATACCGGACATATCTACCTTATCTAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCg                                                        <  1:906859/151‑1 (MQ=255)
                                                                                        gATCGGTGTCGTATTCGGTAAGTTCTACCCACTGCATACCGGACATATCTACCTTATCTAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCg                                                        <  2:115971/151‑1 (MQ=255)
                                                                                        gATCGGTGTCGTATTCGGTAAGTTCTACCCACTGCATACCGGACATATCTACCTTATCTAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCg                                                        <  2:378961/151‑1 (MQ=255)
                                                                                        gATCGGTGTCGTATTCGGTAAGTTCTACCCACTGCATACCGGACATATCTACCTTATCTAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCg                                                        >  1:662934/1‑151 (MQ=255)
                                                                                           cGGTGTCGTATTCGGTAAGTTCTACCCACTGCATACCGGACATATCTACCTTATCTAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCgc                                                       >  2:1259761/1‑149 (MQ=255)
                                                                                           cGGTGTCGTATTCGGTAAGTTCTACCCACTGCATACCGGACATATCTACCTTATCTAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCgc                                                       >  2:1327785/1‑149 (MQ=255)
                                                                                             gtgtCGTATTCGGTAAGTTCTACCCACTGCATACCGGACATATCTACCTTATCTAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCgcagc                                                    >  2:926072/1‑150 (MQ=255)
                                                                                               gtCGTATTCGGTAAGTTCTACCCACTGCATACCGGACATATCTACCTTATCTAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGcagcag                                                  <  1:1197334/150‑1 (MQ=255)
                                                                                               gtCGTATTCGGTAAGTTCTACCCACTGCATACCGGACATATCTACCTTATCTAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCagcagc                                                 <  1:1288224/151‑1 (MQ=255)
                                                                                                  gTATTCGGTAAGTTCTACCCACTGCATACCGGACATATCTACCTTATCTAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGCAGCCg                                               >  1:885139/1‑150 (MQ=255)
                                                                                                   tATTCGGTAAGTTCTACCCACTGCATACCGGACATATCTACCTTATCTAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGCAGCCGa                                              <  1:787320/150‑1 (MQ=255)
                                                                                                    aTTCGGTAAGTTCTACCCACTGCATACCGGACATATCTACCTTATCTAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGCAGCCGAcc                                            >  2:850464/1‑151 (MQ=255)
                                                                                                       cGGTAAGTTCTACCCACTGCATACCGGACATATCTACCTTATCTAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGCAGCCGACCGTg                                         >  1:695207/1‑151 (MQ=255)
                                                                                                       cGGTAAGTTCTACCCACTGCATACCGGACATATCTACCTTATCTAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGCAGCCGACCGTg                                         <  1:354325/151‑1 (MQ=255)
                                                                                                        ggTAAGTTCTACCCACTGCATACCGGACATATCTACCTTATCTAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGCAGCCGACCGTGc                                        <  1:1418721/151‑1 (MQ=255)
                                                                                                        ggTAAGTTCTACCCACTGCATACCGGACATATCTACCTTATCTAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGCAGCCGACCGTGc                                        >  1:363472/1‑151 (MQ=255)
                                                                                                           aaGTTCTACCCACTGCATACCGGACATATCTACCTTATCTAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGCAGCCGACCGTGCCgg                                     >  2:643586/1‑151 (MQ=255)
                                                                                                           aaGTTCTACCCACTGCATACCGGACATATCTACCTTATCTAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGCAGCCGACCGTGCCg                                      <  1:272141/150‑1 (MQ=255)
                                                                                                            aGTTCTACCCACTGCATACCGGACATATCTACCTTATCTAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGCAGCCGACCGTGCCGGa                                    <  2:107492/151‑1 (MQ=255)
                                                                                                             gTTCTACCCACTGCATACCGGACATATCTACCTTATCTAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGCAGCCGACCGTGCCGGAt                                   <  1:1223178/151‑1 (MQ=255)
                                                                                                             gTTCTACCCACTGCATACCGGACATATCTACCTTATCTAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGCAGCCGACCGTGCCGGAt                                   <  2:731273/151‑1 (MQ=255)
                                                                                                              ttCTACCCACTGCATACCGGACATATCTACCTTATCTAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGCAGCCGACCGTGCCGGATc                                  >  2:685971/1‑151 (MQ=255)
                                                                                                               tCTACCCACTGCATACCGGACATATCTACCTTATCTAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGCAGCCGACCGTGCCGGATc                                  <  1:140063/150‑1 (MQ=255)
                                                                                                                   cccACTGCATACCGGACATATCTACCTTATCTAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGCAGCCGACCGTGCCGGATCGTCTg                             <  2:1338348/151‑1 (MQ=255)
                                                                                                                    ccACTGCATACCGGACATATCTACCTTATCTAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGCAGCCGACCGTGCCGGATCGTCTGc                            <  1:213793/151‑1 (MQ=255)
                                                                                                                         gCATACCGGACATATCTACCTTATCTAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGCAGCCGACCGTGCCGGATCGTCTGCGTTgg                       <  2:642718/151‑1 (MQ=255)
                                                                                                                         gCATACCGGACATATCTACCTTATCTAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGCAGCCGACCGTGCCGGATCGTCTGCGTTg                        <  2:966316/150‑1 (MQ=255)
                                                                                                                             aCCGGACATATCTACCTTATCTAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGCAGCCGACCGTGCCGGATCGTCTGCGTTGGTTa                    <  2:680335/150‑1 (MQ=255)
                                                                                                                             aCCGGACATATCTACCTTATCTAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGCAGCCGACCGTGCCGGATCGTCTGCGTTGGTTAt                   >  1:326343/1‑151 (MQ=255)
                                                                                                                             aCCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGCAGCCGACCGTGCCGGATCGTCTGCGTTGGTTAt                   <  2:1368791/151‑1 (MQ=255)
                                                                                                                              ccGGACATATCTACCTTATCTAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGCAGCCGACCGTGCCGGATCGTCTGCGTTGGTTAtt                  <  1:1086297/151‑1 (MQ=255)
                                                                                                                                ggACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGCAGCCGACCGTGCCGGATCGTCTGCGTTGGTTATTGc                <  2:354882/151‑1 (MQ=255)
                                                                                                                                   cATATCTGCCTTATCTAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGCAGCCGACCGTGCCGGATCGTCTGCGTTGGTTATTGCaaa             <  1:642593/151‑1 (MQ=255)
                                                                                                                                     tatCTACCTTATCTAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGCAGCCGACCGTGCCGGATCGTCTGCGTTGGTTATTGCAAACt           >  2:750586/1‑151 (MQ=255)
                                                                                                                                          aCCTTATCTAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGCAGCCGACCGTGCCGGATCGTCTGCGTTGGTTATTGCAAACTTTTa       <  2:67725/150‑1 (MQ=255)
                                                                                                                                           ccTTATCTAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGGCAGTGCCATGTCGCAGCAGCCGACCGTGCCGGATCGTCTGCGTTGGTTATTGCAAACTTTTaa      <  1:733789/150‑1 (MQ=255)
                                                                                                                                              tATCTAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGCAGCCGACCGTGCCGGATCGTCTGCGTTGGTTATTGCAAACTTTTAAatat  >  1:244939/1‑151 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACGATCGGTGTCGTATTCGGTAAGTTCTACCCACTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGCAGCCGACCGTGCCGGATCGTCTGCGTTGGTTATTGCAAACTTTTAAATAT  >  NC_000913/4627421‑4627713

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