Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,406,129 T→C intergenic (+144/‑61) purA → / → nsrR adenylosuccinate synthetase/nitric oxide‑sensitive repressor for NO regulon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,406,1290TC99.0% 288.3 / NA 97intergenic (+144/‑61)purA/nsrRadenylosuccinate synthetase/nitric oxide‑sensitive repressor for NO regulon
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (46/50);  minor base G (1/0);  total (47/50)
Fisher's exact test for biased strand distribution p-value = 4.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAATTCTGGTACGCCTGGCAGATATTTTGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCTTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTACATGGCGTCATTGCCAGAAGGGCGGATGACCAGTATTT  >  NC_000913/4405983‑4406274
                                                                                                                                                  |                                                                                                                                                 
tAATTCAGGTACGCCTGGCAGATAATTTGCCTGCCGGGCGAACAGTGGGATACATTGCTGTGTCGGGTAAACCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCtg                                                                                                                                                 <  2:834174/149‑1 (MQ=255)
     cTGGTACGCCTGGCAGATATTTTGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCt                                                                                                                                            >  1:688395/1‑149 (MQ=255)
     cTGGTACGCCTGGCAGAGTTTTTGCCTGCCGGGCGAACAGTGTGATAGATTGCTGTGTCGGGTAAGCCATTACGCTATCGGACACAGTGTTTAATCCTCGCTTTTTTCCTTCCCCGAAGTGAAATAAATTAGCGACACAGCCTGTGGCt                                                                                                                                            <  1:1138947/149‑1 (MQ=255)
            gCCTGGCAGATATTTTGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTat                                                                                                                                     <  1:245564/149‑1 (MQ=255)
            gCCTGGCAGATATTTTGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTat                                                                                                                                     <  1:141510/149‑1 (MQ=255)
            gCCTGGAAGATATTTTGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTat                                                                                                                                     <  1:1383932/149‑1 (MQ=255)
            gCCGGGCAGATATTTTGCCTGCCGGGAGAACAGTGAGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTat                                                                                                                                     <  1:1389341/149‑1 (MQ=255)
             ccTGGCAGAAATTTTGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTatc                                                                                                                                    <  1:293292/149‑1 (MQ=255)
               tGGCAGATATTTTGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTatcat                                                                                                                                  >  1:1069478/1‑149 (MQ=255)
               tGGCAGATATTTTGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTatcat                                                                                                                                  >  1:735286/1‑149 (MQ=255)
                ggCAGATATTTTGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTatcata                                                                                                                                 >  2:498470/1‑148 (MQ=255)
                 gCAGATATTTTGCCTGGCGTGCGAACAGTTTGAGACATTGCTGTGTCGGGTAATACATTACCGTATCCGACACAGTGTTAAATCCTCTCTTTTTTCCTTCCCCGAACTGAAATAAATTATCTACACAGCCTGTGGCTGGTTTAtcatca                                                                                                                                <  2:989518/149‑1 (MQ=255)
                   aGATATTTTGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACCTCCTGTGGATGGTTTATCATCAat                                                                                                                              >  2:885898/1‑149 (MQ=255)
                    gATATTTTGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCCTTACGCTAGGCGAGAGAGTGTTATATCGTCGCTTTTTTCATTCCCGGAACTGAAATACTTTAGAGCCTGCGCGTGGGGCTGGTTTatctgctata                                                                                                                             >  2:1400130/1‑142 (MQ=255)
                        ttttGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATaaa                                                                                                                         >  1:946612/1‑149 (MQ=255)
                         tttGCCTGCCGGGCGAACAGTGTGATAGATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAAt                                                                                                                        <  2:122310/149‑1 (MQ=255)
                         tttGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAAt                                                                                                                        <  1:412137/149‑1 (MQ=255)
                         tttGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAAt                                                                                                                        <  2:435602/149‑1 (MQ=255)
                            gcctgccGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTa                                                                                                                     <  2:655158/149‑1 (MQ=255)
                               tgccGGGCGAACAGTGTGATACATTGCTCTGTCGTGTAAGCCATTACGCAATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTAtat                                                                                                                  <  2:1233661/149‑3 (MQ=255)
                                gccGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTAtttt                                                                                                                 <  2:543456/149‑1 (MQ=255)
                                       gAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGa                                                                                                          <  2:1374275/149‑1 (MQ=255)
                                           aGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTc                                                                                                      >  1:139256/1‑149 (MQ=255)
                                           aGTGTGATACATTGCTGAGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTc                                                                                                      <  2:1069478/149‑1 (MQ=255)
                                            gtgtGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTcc                                                                                                     <  2:557342/149‑1 (MQ=255)
                                              gtgATACATTGCTGTGTCGGGTAAGCCATTACGCTATGCGACAAAGTGTTAAATCCTCGCGTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACATCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCt                                                                                                   <  1:885644/149‑1 (MQ=255)
                                               tgATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCtt                                                                                                  >  2:664038/1‑149 (MQ=255)
                                                gATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCttt                                                                                                 >  1:1290840/1‑149 (MQ=255)
                                                 aTACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCtttt                                                                                                <  2:431175/149‑1 (MQ=255)
                                                 aTACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCtttt                                                                                                <  2:1443296/149‑1 (MQ=255)
                                                  tACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACCGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTGTTg                                                                                               >  1:1115531/1‑149 (MQ=255)
                                                    cATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAg                                                                                             >  2:629734/1‑149 (MQ=255)
                                                     aTTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGa                                                                                              <  2:1163431/147‑1 (MQ=255)
                                                     aTTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGa                                                                                              >  1:1163431/1‑147 (MQ=255)
                                                            tgtCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTGCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAATTGTATTTTTTCCCGATTTCCCTTTTGTGGTTGAtgt                                                                                     >  1:1156182/1‑149 (MQ=255)
                                                             gtCGGGTAAGCCATTACGCTATCAGACACAGAGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGAGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCAGATTTCCCTTTTGAGGTTGAtgtg                                                                                    >  2:744313/1‑149 (MQ=255)
                                                               cGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCa                                                                                  <  2:1146860/149‑1 (MQ=255)
                                                                  ttAAGACATTACGCTATCAGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATATATTAGCGACACAGCCTGTGGCTGGTTTATCATCATTATAAATGTATTTTTTCCCGTTTTCCCTTTTGAGGTTGATGTGCAGtt                                                                               <  1:758630/148‑1 (MQ=255)
                                                                  gTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGtt                                                                               >  1:1349144/1‑149 (MQ=255)
                                                                  gTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATAAGAGACACAGCCTGTGGCTGGTTTATCAGCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGtt                                                                               <  2:393970/149‑1 (MQ=255)
                                                                       ccATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGa                                                                          >  2:721757/1‑149 (MQ=255)
                                                                       ccATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGa                                                                          >  2:281107/1‑149 (MQ=255)
                                                                        cATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAg                                                                         >  1:1331235/1‑149 (MQ=255)
                                                                         aTTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGGTGATGTGCAGTTAACGAGt                                                                        <  2:1410154/149‑1 (MQ=255)
                                                                          ttACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGtt                                                                       >  2:303495/1‑149 (MQ=255)
                                                                             cGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCa                                                                    >  1:858998/1‑149 (MQ=255)
                                                                              gCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCAc                                                                   <  1:721757/149‑1 (MQ=255)
                                                                              gCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCAc                                                                   >  2:67914/1‑149 (MQ=255)
                                                                               caaTCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGCCACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTAGAGGTTGATGTGCAGTTAACGAGTTTCACt                                                                  <  1:337019/147‑1 (MQ=255)
                                                                               cTATCCGACACAGTGTTAAATCCTCGGTTTTTTCCATCCCCGAACTGAAATAAATTAGCGACACATCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGATAACGAGTTTCACt                                                                  <  1:331508/149‑1 (MQ=255)
                                                                               cTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACt                                                                  <  2:308237/149‑1 (MQ=255)
                                                                               atATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACAAATCCTTTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACt                                                                  <  1:498470/148‑1 (MQ=255)
                                                                                tATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTg                                                                 >  2:1049109/1‑149 (MQ=255)
                                                                                tATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTg                                                                 >  1:455695/1‑149 (MQ=255)
                                                                                     gACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTgattac                                                            <  1:1149791/149‑1 (MQ=255)
                                                                                       cacaGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTgattacgg                                                          >  2:577569/1‑149 (MQ=255)
                                                                                       cacaGTGTTAAATCCTCGCTTTTTTCCATCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTgattacgg                                                          <  2:89428/149‑1 (MQ=255)
                                                                                         caGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTgattacggat                                                        <  1:1221908/149‑1 (MQ=255)
                                                                                           gtgaTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGATTATCATCAATATAAATCTAATTATACCCTATTTCCATTTTCAGGTTGATGTGCAGTTAACGAGTTTCAATgattacggataa                                                      >  2:633066/1‑147 (MQ=255)
                                                                                             gtTAAATCCTCGC‑TTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGt                                                   >  1:265218/1‑149 (MQ=255)
                                                                                              ttAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGt                                                   >  2:1129126/1‑149 (MQ=255)
                                                                                               tAAATCCTTGCTTTTTTCCTTCCCCGTACTTAAATAAATTATCGTCACAGCCGGTGGCTGGTTTATAATCAATATAAATTTATTTTTTCCCGATTTCCCTTTTGAGGTTTATGTGCATTTAAAGATTTTCACTGATTACGGATTACGTg                                                  <  2:688395/149‑1 (MQ=255)
                                                                                                aaaTCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTgc                                                 >  1:1325882/1‑149 (MQ=255)
                                                                                                aaaTCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTgc                                                 >  2:907831/1‑149 (MQ=255)
                                                                                                 aaTCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTgcg                                                >  1:1240620/1‑149 (MQ=255)
                                                                                                    ccTCGC‑TTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGa                                            <  2:265218/149‑1 (MQ=255)
                                                                                                     cTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGa                                            >  2:1408019/1‑149 (MQ=255)
                                                                                                      tCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACTTTCGCTGAt                                           >  1:1361503/1‑149 (MQ=255)
                                                                                                      tCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGAt                                           >  1:1334615/1‑149 (MQ=255)
                                                                                                      tCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGAt                                           >  1:481194/1‑149 (MQ=255)
                                                                                                      tCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGAt                                           <  1:409188/149‑1 (MQ=255)
                                                                                                       cGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATc                                          <  1:1331834/149‑1 (MQ=255)
                                                                                                          ttttttCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTAc                                       <  1:629734/149‑1 (MQ=255)
                                                                                                              ttccttccCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTACATgg                                   >  1:499176/1‑149 (MQ=255)
                                                                                                              ttccttccCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTACATgg                                   >  1:1288057/1‑149 (MQ=255)
                                                                                                                ccttccGCGCACTGAAATAAAATAACTACACAGCCTTTGTCTGGTGTATCATCAATATAAATGTATTTTTTCCCGATTTACATTTTGAGGATGATGTGCCGTTAACGAGTTTCCCTGATTTCTGATTCCGTGCTCTGATCTACATGGCg                                 <  2:62705/149‑1 (MQ=255)
                                                                                                                  ttccCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTACATGGCGTc                               >  1:1183317/1‑149 (MQ=255)
                                                                                                                  ttccCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATAATCAATATAAATGTATTTTTTCCCGATTTCCCTTTAGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTACATGGCGTc                               <  1:885898/149‑1 (MQ=255)
                                                                                                                     cccGAACTGAAATAAATTAGCGACACAGCCTTTGGCTGGTTTATCATCAAAATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTACATGGCGTCAtt                            <  2:1420208/149‑1 (MQ=255)
                                                                                                                     cccGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTACATGGCGTCAtt                            >  1:192557/1‑149 (MQ=255)
                                                                                                                      ccGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTACATGGCGTCATTg                           <  2:768595/149‑1 (MQ=255)
                                                                                                                       ataaCTTAAATAAATTAGCGACAAAGCCTGTGTCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGATTTTCACTGATTACGGATTACTTGCGCTGATCTACATGGCGTCATTGc                          <  2:864591/147‑1 (MQ=255)
                                                                                                                            tGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTACATGGCGTCATTGCCAGaa                     <  2:858998/149‑1 (MQ=255)
                                                                                                                              aaataaatTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTACATGGCGTCATTGCCAGAAgg                   >  2:618095/1‑149 (MQ=255)
                                                                                                                                ataaatTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTACATGGCGTCATTGCCAGAAGGGc                 >  1:856465/1‑149 (MQ=255)
                                                                                                                                  aaatTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTACATGGCGTCATTGCCAGAAGGGCgg               >  1:20239/1‑149 (MQ=255)
                                                                                                                                    atTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTACATGGCGTCATTGCCAGAAGGGCGGAt             <  1:595910/149‑1 (MQ=255)
                                                                                                                                       aGCGACACAGGCTGTGGCTGGTTTATAATCAATATTAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTACATGGCGTCATTGCCAGAAGGGCGGATGAc          <  2:1349144/149‑1 (MQ=255)
                                                                                                                                       aGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTACATGGCGTCATTGCCAGAAGGGCGGATGAc          <  1:281107/149‑1 (MQ=255)
                                                                                                                                          gACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTACATGGCGTCATTGCCAGAAGGGCGGATGACCAg       >  2:1196155/1‑149 (MQ=255)
                                                                                                                                          gACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTACATGGCGTCATTGCCAGAAGGGCGGATGACAAg       >  2:408321/1‑149 (MQ=255)
                                                                                                                                           acacAGCCTGTGGCTGGTTTATCATCAATATAATTGTATTTTTTCCAGATTTCGGTTTTGTGGTTGATGTGGTGTTAACGAGTTTCACTGATTTCGGTTTACTTGAGGTGATCTGGATGGAGTGATTGCCAGAATGTCGGAGGGCCatt      >  2:1430249/1‑147 (MQ=255)
                                                                                                                                               aGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTACATGGCGTCATTGCCAGAAGGGCGGATGACCAGTAttt  >  2:1021781/1‑149 (MQ=255)
                                                                                                                                               aGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTACATGGCGTCATTGCCAGAAGGGCGGATGACCAGTAttt  >  1:669032/1‑149 (MQ=255)
                                                                                                                                               aGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTACATGGCGTCATTGCCAGAAGGGCGGATGACCAGTAttt  <  1:67914/149‑1 (MQ=255)
                                                                                                                                               aGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTACATGGCGTCATTGCCAGAAGGGCGGATGACCAGTAttt  <  1:744313/149‑1 (MQ=255)
                                                                                                                                               aGCCTGTGGCGGGTTTATCATCAATAAAAATGTATTTTTTCCCGATTTCCCTATTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTACATGGCGTCATTGCCAGAAGGGCGGATGACCAGTAttt  <  1:1340305/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TAATTCTGGTACGCCTGGCAGATATTTTGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCTTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTACATGGCGTCATTGCCAGAAGGGCGGATGACCAGTATTT  >  NC_000913/4405983‑4406274

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