Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,181,603 C→T 6.7% intergenic (+22/+16) ygiC → / ← ygiD ATP‑Grasp family ATPase/4,5‑ DOPA‑extradiol‑dioxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,181,6030CT 6.7% 425.4 / 11.6 135intergenic (+22/+16)ygiC/ygiDATP‑Grasp family ATPase/4,5‑ DOPA‑extradiol‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base C (71/54);  new base T (3/6);  total (74/61)
Fisher's exact test for biased strand distribution p-value = 2.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 138 total aligned reads displayed.

GCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGTTCG  >  NC_000913/3181454‑3181752
                                                                                                                                                     |                                                                                                                                                     
gCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTc                                                                                                                                                       >  1:1755694/1‑150 (MQ=255)
 cGACAGCTATATGCTGATTGGTAGCTGGCGGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCg                                                                                                                                                      <  2:1013889/150‑1 (MQ=255)
      gCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATc                                                                                                                                                 >  1:1582830/1‑150 (MQ=255)
       cTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATcc                                                                                                                                                >  1:1647399/1‑150 (MQ=255)
         atatTCTGATGGGTAGCTGGCTGGTGAACGCTCAACCCGCCGGCATTGGCATTCGGGAAGACCGTGCAGTGATCACCCAGGATAGGGCGAGGGTTTATCCACATATTTTTGGTGAATACGCCACGAAACCGGATGGCGAGTGCCATCCgg                                                                                                                                              <  1:1644442/150‑1 (MQ=255)
         atatGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGTAATTGGCAGTCGTGTAGATCGTGCATTGATCACCCAGGATATGTCTCGGGTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCGAGTGCCATCCgg                                                                                                                                              <  1:1059755/150‑1 (MQ=255)
         atatGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCgg                                                                                                                                              <  1:261356/150‑1 (MQ=255)
          tatGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGt                                                                                                                                             >  1:1266812/1‑150 (MQ=255)
          tatGCTGATTGGTAGCTGGCTGGAGACCGATCAAACCGCCGGAATTGGCATTCGGGAAGACCGTTCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCGAGTGCCATCCGGt                                                                                                                                             <  1:886412/150‑1 (MQ=255)
            tGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTaa                                                                                                                                           >  2:1597730/1‑150 (MQ=255)
            tGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTaa                                                                                                                                           >  2:1314093/1‑150 (MQ=255)
            tGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTaa                                                                                                                                           >  2:1805/1‑150 (MQ=255)
             gCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAAt                                                                                                                                          >  1:1103591/1‑150 (MQ=255)
              cTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAAtt                                                                                                                                         >  2:868593/1‑150 (MQ=255)
               ggATTGGTAGCTGGCTGGTGAACGATCACCCCGCCGGAATGGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTg                                                                                                                                        <  2:788228/149‑1 (MQ=255)
                  ttGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGCTCACCCAGGCTATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCCCGATACCGGATGGCACTCGCCATCCGGTAATTGTTa                                                                                                                                     <  1:393300/150‑1 (MQ=255)
                   tGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAg                                                                                                                                    <  1:743621/150‑1 (MQ=255)
                   tGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAg                                                                                                                                    <  2:1081842/150‑1 (MQ=255)
                        gctggctgGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTAt                                                                                                                               >  1:557745/1‑150 (MQ=255)
                        gctggctgGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATAACCCAGGATATGTCTCGGTATTATCCACATATTTTTGTTGAATAAGCAACGATACCGGATGGCTCTCTACATACGGTAATTGTTAGACTAt                                                                                                                               >  1:240310/1‑150 (MQ=255)
                         ctggctggTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATc                                                                                                                              <  2:1236664/150‑1 (MQ=255)
                         ctggctggTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATc                                                                                                                              <  1:1265649/150‑1 (MQ=255)
                           ggctggTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCAACCGGTAATTGTTATCCTATCTg                                                                                                                            >  1:1623803/1‑150 (MQ=255)
                           ggctggTGAACGAGCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTg                                                                                                                            >  2:674440/1‑150 (MQ=255)
                            gctggTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGc                                                                                                                           >  1:622459/1‑150 (MQ=255)
                               ggTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCATATCTGCAcc                                                                                                                        >  1:373225/1‑150 (MQ=255)
                                  gAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGAc                                                                                                                     >  2:1022084/1‑150 (MQ=255)
                                  gAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGAc                                                                                                                     >  2:1015534/1‑150 (MQ=255)
                                    aCGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAg                                                                                                                   >  2:1590303/1‑150 (MQ=255)
                                    aCGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAg                                                                                                                   >  1:533360/1‑150 (MQ=255)
                                      gATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACTGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCa                                                                                                                 <  2:424662/150‑1 (MQ=255)
                                      gATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCa                                                                                                                 <  1:1863609/150‑1 (MQ=255)
                                         cAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTAGTTGAATAAGCCACGATACCGGATGGCACTCGACATCCGGTAATTGTAAGCCTATATGCACCGACAGCATaa                                                                                                              >  2:414200/1‑149 (MQ=255)
                                            cccGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCa                                                                                                           <  1:978436/150‑1 (MQ=255)
                                             ccgccgGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAg                                                                                                          <  1:1326166/150‑1 (MQ=255)
                                             ccgccgGAATTGGCATTCGTGAAGACAGTGCATTGATCACCCAGGATATGTCTCGGGTTTATCCACATATTTTTGTTGAATATGCCACGATACCGGAAGACACTCTCCATCCGGTAAATGTTAGCCTATCTAAACCGACATCATACTCAg                                                                                                          >  1:321047/1‑150 (MQ=255)
                                               gccgGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGc                                                                                                        <  1:793894/150‑1 (MQ=255)
                                                     tttGGCATTCGTGAAGAGCGTGCATTGATCACCCAGGAGATGTATCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGAGGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCa                                                                                                  <  1:792674/149‑1 (MQ=255)
                                                      ttGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCAt                                                                                                 >  1:868507/1‑150 (MQ=255)
                                                         gCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTc                                                                                              <  2:753623/150‑1 (MQ=255)
                                                         gCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTc                                                                                              <  2:1546949/150‑1 (MQ=255)
                                                          cATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCt                                                                                             >  1:1282814/1‑150 (MQ=255)
                                                          cATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCt                                                                                             <  2:1045366/150‑1 (MQ=255)
                                                          cATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCt                                                                                             <  2:1721145/150‑1 (MQ=255)
                                                           aTTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCta                                                                                            >  1:1139399/1‑150 (MQ=255)
                                                             tCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCtata                                                                                          <  1:1232463/150‑1 (MQ=255)
                                                              cGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATAc                                                                                         >  2:1108149/1‑150 (MQ=255)
                                                              cGTGAAGACCGGGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATAc                                                                                         <  2:1398039/150‑1 (MQ=255)
                                                                tGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATAccc                                                                                       >  2:1423701/1‑150 (MQ=255)
                                                                      ccGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAAcc                                                                                 >  2:832075/1‑150 (MQ=255)
                                                                       cGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCg                                                                                <  2:1488257/150‑1 (MQ=255)
                                                                          gCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGaa                                                                             >  2:1824043/1‑150 (MQ=255)
                                                                                  cACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAAtt                                                                     <  1:84942/150‑1 (MQ=255)
                                                                                   aCCCAGGATATGTCTCGGTTTTATCCACATTTTTTTGTTGAATAAGCCACGATTCCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTg                                                                    >  1:1121416/1‑150 (MQ=255)
                                                                                   aCCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTg                                                                    >  2:1178951/1‑150 (MQ=255)
                                                                                   aCCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCAACGATACCTGATGGCACTCGACATCCGGTAATTGTTAGCCTATCTGCACAGACAGCATACTCAGGCTGCCCATTTCTATACCCTAAACCGTAATGGTAATTg                                                                    >  1:1223500/1‑150 (MQ=255)
                                                                                    cccAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCCCTCGCAATCCGGTAATTGTTAGCCTATCTGCACCGACAGAATACTCAGGCTGCCCAATTCTATACCCTCAACCTGAATGGTAATTgg                                                                   >  1:1342382/1‑150 (MQ=255)
                                                                                     ccAGGTTATGGCTCGGTTTTATCAACATATTTTTGTTGAATAATCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGc                                                                  <  1:495414/150‑1 (MQ=255)
                                                                                     ccAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGc                                                                  <  1:464726/150‑1 (MQ=255)
                                                                                           tatGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCAACCGGTAATTGTTAGCATATCAGCACAGACAGCAGACTCAGGCTGCCCATTTCTATCCCCTAAACCGGAATGCTAATTGGATCCTGc                                                            >  2:1611591/1‑150 (MQ=255)
                                                                                           tatGTCTCGGTTGTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCCACCGGAATGGTAATTGGCTCCTGc                                                            <  1:1436786/150‑1 (MQ=255)
                                                                                              gTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTCCCCa                                                         >  1:301186/1‑150 (MQ=255)
                                                                                                ctcGGTTTTTTCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAAGCGGATTGGTAATTGGCGGCTGCCCTTc                                                       >  2:1527057/1‑150 (MQ=255)
                                                                                                ctcGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATc                                                       >  1:1084620/1‑150 (MQ=255)
                                                                                                   ggTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCa                                                    <  2:1647399/150‑1 (MQ=255)
                                                                                                    gTTTTATCCACATTTTTTTGTTTAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGTCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCAc                                                   <  2:275495/150‑1 (MQ=255)
                                                                                                    gTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCAc                                                   <  1:1373724/150‑1 (MQ=255)
                                                                                                      tttATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGc                                                 >  2:1865239/1‑150 (MQ=255)
                                                                                                         aTCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCCTCCCACGCAcc                                              >  1:610295/1‑150 (MQ=255)
                                                                                                         aTCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCAATACCACGCAcc                                              >  2:1177213/1‑150 (MQ=255)
                                                                                                           ccACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTa                                            <  2:1274535/150‑1 (MQ=255)
                                                                                                             acaTATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAac                                          <  1:511768/150‑1 (MQ=255)
                                                                                                               atatTTTTGTTTTATAAGCCTCTATACCGGATGGCTCTCGCCATCCGGAAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAacac                                        <  2:285839/150‑1 (MQ=255)
                                                                                                                tatTTTTGTTGAATTAGCCACGATACCGGATGTCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCTGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAAcaca                                       <  2:1295442/150‑1 (MQ=255)
                                                                                                                tatTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCCCCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGACCATCCCACGCACCTAAcaca                                       >  1:1428058/1‑150 (MQ=255)
                                                                                                                   ttttGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATac                                    <  2:1755694/150‑1 (MQ=255)
                                                                                                                         tgAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACaacaac                              >  2:87706/1‑150 (MQ=255)
                                                                                                                         tgAATAAGCCACGATACCGGATGCCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACaacaac                              <  2:206986/150‑1 (MQ=255)
                                                                                                                          gAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCACCCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACg                             >  2:1571432/1‑150 (MQ=255)
                                                                                                                          gAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACg                             >  2:1578985/1‑150 (MQ=255)
                                                                                                                          gAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACg                             >  2:853154/1‑150 (MQ=255)
                                                                                                                          gAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACg                             >  2:516013/1‑150 (MQ=255)
                                                                                                                            aTAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCAGCACAACCCACGCACCTAACACATACAACAAAGGc                           >  2:1281415/1‑150 (MQ=255)
                                                                                                                              aaGCCACGATACCGGATGGCACTCGCCATCCGGTAAGTGTTAGCCTCTCTGCTCCGGCAGCATACTGAGGCTGCCCTTTTCTATAGCCTCAGCCGGAATGGTGATTGGCTGCTGCCCAGCGCACGCGCGTAGCACAGACAACAACGGCaa                         >  2:1614839/1‑150 (MQ=255)
                                                                                                                               aGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCaaa                        >  2:1319724/1‑150 (MQ=255)
                                                                                                                               aGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCaaa                        >  1:868658/1‑150 (MQ=255)
                                                                                                                                gCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCCCGCACCTAACACATACAACAACGGCAaat                       >  1:1648907/1‑150 (MQ=255)
                                                                                                                                 ccACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACACCGGCAaata                      >  1:1014033/1‑150 (MQ=255)
                                                                                                                                     gATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGc                  <  2:903171/150‑1 (MQ=255)
                                                                                                                                       tACCGGATGGCTCTTGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCACATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGctc                <  2:1869848/150‑1 (MQ=255)
                                                                                                                                        aCCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATTCCCTCTACCGGAATGGTAATTGGCTCCAGCCAAACCCACGCACCGAACACATACAACAACGGCAAAGAATGctct               >  2:1719412/1‑150 (MQ=255)
                                                                                                                                        aCCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGctct               >  2:1028513/1‑150 (MQ=255)
                                                                                                                                         ccGGATGGCGAGTGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACCCATACATCCACGGCAAATAATGCTCTg              >  1:687253/1‑150 (MQ=255)
                                                                                                                                         ccGGATGGCGAGTGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACATCAACGGCAAATAATGCTCTg              >  1:1690646/1‑150 (MQ=255)
                                                                                                                                          cGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTgg             <  1:1590303/150‑1 (MQ=255)
                                                                                                                                            gATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCg           <  2:1569914/150‑1 (MQ=255)
                                                                                                                                               ggCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTg        <  2:625873/150‑1 (MQ=255)
                                                                                                                                                gCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGGTgg       >  1:851982/1‑150 (MQ=255)
                                                                                                                                                     cGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGTTCg  <  2:557745/150‑1 (MQ=255)
                                                                                                                                                     cGCCATCCGGTAATTGTTAGACTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGTTCg  <  1:414200/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGTTCG  >  NC_000913/3181454‑3181752

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: