Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,592,875 C→T 5.3% pseudogene (1251/1461 nt) yneO ← pseudogene, AidA homolog

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,592,8750CT 5.3% 237.4 / 2.7 76pseudogene (1251/1461 nt)yneOpseudogene, AidA homolog
Reads supporting (aligned to +/- strand):  ref base C (36/35);  new base T (2/2);  total (39/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.16e-01

ACCGACATTCCCCCACAGGCTAAGGTTGTTATTTACTTTCGCCTCAACCCCGGTACGTACTTCACCGAGATTACGCGCACCTTCACGGCCTACGGTTTGACCATTCATCTTCACGGCGTAGACTTTGCTATTGTTGATCCAGTTCGCTTCAATGTAAGGCTGGAACTCACGCTGTTTACCATCGTCACGCTGGTGATGGCTGTTCAGGTAGGTTTTCACCCCAAGTCGCGTTTGCACATTTCCGTCGCCTTCCGTTTCAATGCGCGTTCCGTCTTTCCGGGTATGGTCGGAA  >  NC_000913/1592731‑1593022
                                                                                                                                                |                                                                                                                                                   
acCGACATTCCCCCACAGGCTAAGGTTGTTATTTACTTTCGCCTCAACCCCGGTACGTACTTCACCGAGATTACGCGCACCTTCACGGCCTACGGGTTGACCATTCATCTTCACGGCGTAGACTTTGCTATTGTTGATCCAGTTCGCtt                                                                                                                                                 >  1:127800/1‑149 (MQ=255)
  cGACATTCCCCCACAGGCTAAGGTTGTTATTTACTTTCGCCTCAACCCCGGTACGTACTTCACCGAGATTACGCGCACCTTCACGGCCTACGGTTTGACCATTCATCTTCACGGCGTAGACTTTGCTATTGTTGATCCAGTTCGCTTCa                                                                                                                                               <  6:311413/149‑1 (MQ=255)
    aCATTCCCCCACAGGCTAAGGTTGTTATTTACTTTCGCCTCAACCCCGGTACGTACTTCACCGAGATTACGCGCACCTTCACGGCCTACGGTTTGACCATTCATCTTCACGGCGTAGACTTTGCTATTGTTGATCCAGTTCGCTTCAAt                                                                                                                                             >  8:278823/1‑149 (MQ=255)
    aCATTCCCCCACAGGCTAAGGTTGTTATTTACTTTCGCCTCAACCCCGGTACGTACTTCACCGAGATTACGCGCACCTTCACGGCCTACGGTTTGACCATTCATCTTCACGGCGTAGACTATGCTATTGTTGATCCAGTTCGCTTCAAt                                                                                                                                             >  8:81611/1‑149 (MQ=255)
     cATTCCCCCACAGGCTAAGGTTGTTATTTACTTTCGCCTCAACCCCGGTACGTACTTCACCGAGATTACGCGCACCTTCACGGCCTACGGTTTGACCATTCATCTTCACGGCGTAGACTTTGCTATTGTTGATCCAGTTCGCTTCAATg                                                                                                                                            <  7:278823/149‑1 (MQ=255)
      aTTCCCCCACAGGCTAAGGTTGTTATTTTCTTTCGCCTCAACCCCGGTACGTACTTCACCGAGATTACGCGCACCTTCACGGCCTACGGTTTGACCATTCATCTTCACGGCGTAGACTTTGCTATTGTTGATCCAGTTCGCTTCAATg                                                                                                                                            >  2:405486/1‑148 (MQ=255)
      aTTCCCCCACAGGCTAAGGTTGTTATTTACTTTCGCCTCAACCCCGGTACGTACTTCACCGAGATTACGCGCACCTTCACGGCCTACGGTTTGACCATTCATCTTCACGGCGTAGACTTTGCTATTGTTGATCCAGTTCGCTTCAATg                                                                                                                                            >  7:293276/1‑148 (MQ=255)
       ttCCCCCACAGGCTAAGGTTGTTATTTACTTTCGCCTCAACCCCGGTACGTACTTCACCGAGATTACGCGCACCTTCACGGCCTACGGTTTGACCATTCATCTTCACGGCGTAGACTTTGCTATTGTTGATCCAGTTCGCTTCAATGTa                                                                                                                                          <  5:314051/149‑1 (MQ=255)
       ttCCCCCACAGGCTAAGGTTGTTATTTACTTTCGCCTCAACCCCGGTACGTACTTCACCGAGATTACGCGCACCTTCACGGCCTACGGTTTGACCATTCATCTTCACGGCGTAGACTTTGCTATTGTTGATCCAGTTCGCTTCAATGTa                                                                                                                                          <  7:81611/149‑1 (MQ=255)
         cccccACAGGCTAAGGTTGTTATTTACTTTCGCCTCAACCCCGGTACGTACTTCACCGAGATTACGCGCACCTTCACGGCCTACGGTTTGACCATTCATCTTCACGGCGTAGACTTTGCTATTGTTGATCCAGTTCGCTTCAATGTaa                                                                                                                                         <  2:16745/148‑1 (MQ=255)
         ccaccACAGGCTAAGGTTGTTATTTACTTTCGCCTCAACCCCGGTACGTACTTCACCGAGATTACGCGCACCTTCACGGCCTACGGTTTGACCATTCATCTTCACGGCGTAGACTTTGCTATTGTTGATCCAGTTCGCTTCAATGTAAg                                                                                                                                        >  8:151130/4‑149 (MQ=255)
           cccACAGTCTAAGGTTGTTATTTACTTTCGCCTCAACCCCGGTACGTACTTCACCTAGATTACGCGCACCTTCACGGCCTACGGATTGACCATTCATCTTCACGGCGTAGACTATGCTATTGTTGATCCAGTACGCTTCAATGTAAGGc                                                                                                                                      >  8:267991/1‑149 (MQ=255)
              acaGGCTAAGGTTGTTATTTACTTTCGCCTCAACCCCGGTACGTACTTCACCGAGATTACGCGCACCTTCACGGCCTACGGTTTGACCATTCATCTTCACGGCGTAGACTTTGCTATTGTTGATCCAGTTTGCTTCAATGTAAGGCTg                                                                                                                                    >  8:191823/1‑148 (MQ=255)
                aGGCTAAGGTTGTTATTTACTTTCGCCTCAACCCCGGTACGTACTTCACCGAGATTACGCGCACCTTCACGGCCTACGGTTTGACCATTCATCTTCACGGCGTAGACTTTGCTATTGTTGATCCAGTTCGCTTCAATGTAAGGCTGGaa                                                                                                                                 <  8:330555/149‑1 (MQ=255)
                  gCTAAGGTTGTTATTTACTTTCGCCTCAACCCCGGTACGTACTTCACCGAGATTACGCGCACCTTCACGGCCTACGGTTTGACCATTCATCTTCACGGCGTAGACTTTGCTATTGTTGATCCAGTTCGCtt                                                                                                                                                 >  4:84272/1‑131 (MQ=255)
                   cTAAGGTTGTTATTTACTTTCGCCTCAACCCCGGTACGTACTTCACCGAGATTACGCGCACCTTCACGGCCTACGGTTTGACCATTCATCTTCACGGCGTAGACTTTGCTATTGTTGATCCAGTTCGCTTCAATGTAAGGCTGGAACt                                                                                                                               >  6:36056/1‑148 (MQ=255)
                   cTAAGGTTGTTATTTACTTTCGCCTCAACCCCGGTACGTACTTCACCGAGATTACGCGCACCTTCACGGCCTACGGTTTGACCATTCATCTTCACGGCGTAGACTTTGCTATTGTTGATCCAGTTCGCTTCAATGTAAGGCTGGAACt                                                                                                                               >  8:69941/1‑148 (MQ=255)
                      aGGTTGTTATTTACTTTCGCCTCAACCCCGTTACGTACTTCACCGAGATTACCCGCACCTTCACGGCCTACGGTTTGACCATTCATCTTCACGGCGTAGACTTTGCTATTGTTGATCCAGTTCGCTTCAATGTAAGGCTGGAAATCACg                                                                                                                           >  7:115741/1‑149 (MQ=255)
                      aGGTTGTTATTTACTTTCGCCTCAACCCCGGTACGTACTTCACCGAGATTACGCGCACCTTCACGGCCTACGGTTTGACCATTCATCTTCACGGCGTAGACTTTGCTATTGTTGATCCAGTTCGCTTCAATGTAAGGCTGGAACTCAc                                                                                                                            <  1:88723/148‑1 (MQ=255)
                            ttATTTACTTTCGCCTCAACCCCGGTACGTACTTCACCGAGATTACGCGCACCTTCACGGCCTACGGTTTGACCATTCATCTTCACGGCGTAGACTTTGCTATTGTTGATCCAGTTCGCTTCAATGTAAGGCTGGAACTCACGCTGtt                                                                                                                      <  8:293276/148‑1 (MQ=255)
                              aTTTACTTTCGCCTCAACACCGGTACGTACTTCACCGAGATTACGCGCACCTTCACGGCCTACGGTTTGACCATTCATCTTCACGGCGTAGACTTTGCTATTGTTGATCCAGTTCGCTTCAATGTAAGGCTGGAACTCACGCTGTTTAc                                                                                                                   >  2:34386/1‑149 (MQ=255)
                               tttACTTTCGCCTCAACCCCGGTACGTACTTCACCGAGATTACGCGCACCTTCACGGCCTACGGTTTGACCATTCATCTTCACGGCGTAGACTTTGCTATTGTTGATCCAGTTCGCTTCAATGTAAGGCTGGAACTCACGCTGTTTAc                                                                                                                   >  4:29699/1‑148 (MQ=255)
                                ttACTTTCGCCTCAACCCCGGTACGTACTTCACCGAGATTACGCGCACCTTCACGGCCTACGGTTTGACCATTCATCTTCACGGCGTAGACTTTGCTATTGTTGATCCAGTTCGCTTCAATGTAAGGCTGGAACTCACGCTGTTTACCa                                                                                                                 >  7:105562/1‑149 (MQ=255)
                                ttACTTTCGCCTCAACCCCGGTACGTACTTCACCGAGATTACGCGCACCTTCACGGCCTACGGTTTGACCATTCATCTTCACGGCGTAGACTTTGCTATTGTTGATCCAGTTCGCTTCAATGTAAGGCTGGAACTCACGCTGTTTACCa                                                                                                                 >  8:81845/1‑149 (MQ=255)
                                 tACTTTCGCCTCAACCCCGGTACGTACTTCACCGAGATTACGCGCACCGTCACGGCCTACGGTTTGACCATTCATCTTCACGGCGTAGACTTTGCTATTGTTGATCCAGTTCGCTTCAATGTAAGGCTGGAACTCACGCTGTTTACCAt                                                                                                                >  3:100086/1‑149 (MQ=255)
                                      tCGCGTCAACACCGGTACGTACTTCAACGAGATTACTCGCACCTTCACGGCCTACGGTTTAACCATTGATCTTCACGGCGTAGACTTTGCTATTGTTGATCCAGTTCGCTTCAATGTAAGGCTGGAACTCACGCTGTTTACCATCGTc                                                                                                            <  4:275071/148‑1 (MQ=255)
                                      tCGCCTCAACCCCGGTACGTACTTCACCGAGATTACGCGCACCTTCACGGCCTACGGTTTGACCATTCATCTTCACGGCGTAGACTTTGCTATTGTTGATCCAGTTCGCTTCAATGTAAGGCTGGAACTCACGCTGTTTACCATCGt                                                                                                             <  7:151130/147‑1 (MQ=255)
                                      tCGCCTCAACCCCGGTACGTACTTCACCGAGATTACGCGCACCTTCACGGCCTACGGTTTGACCATTCATCTTCACGGCGTAGACTTTGCTATTGTTGATCCAGTTCGCTTCAATGTAAGGCTGGAACTCACGCTGTTTACCATCGt                                                                                                             <  5:36056/147‑1 (MQ=255)
                                      tCGCCTCAACCCCGGTACGTACTTCACCGAGATTACGCGCACCTTCACGGCCTACGGTTTGACCATTCATCTTCACGGCGTAGACTTTGCTATTGTTGATCCAGTTCGCTTCAATGTAAGGCTGGAACTCACGCTGTTTACCATCGt                                                                                                             <  7:69941/147‑1 (MQ=255)
                                        gCCTCAACCCCGGTACGTACTTCACCGAGATTACGCGCACCTTCACGGCCTACGGTTTGACCATTCATCTTCACGGCGTAGACTTTGCTATTGTTGATCCAGTTCGCTTCAATGTAAGGCTGGAACTCACGCTGTTTACCATCGTCAc                                                                                                          <  1:405486/148‑1 (MQ=255)
                                         ccTCAACCCCGGTACGTACTGCAACGAGATTACGCGCACCTTCACGGCCTACGGGTTGACCATTCATCTTCACGGCGTAGACTTTGCTATTGTTGATCCAGTTAGCTTCAATGTAAGGATGGAACTCACGCTGTTTAACATCGTCACGc                                                                                                        >  7:48240/1‑149 (MQ=255)
                                          cTCAACCCCGGTACGTACTTCACCGAGATTACGCGCACCTTCACGGCCTACGGTTTGACCATTCATCTTCACGGCGTAGACTTTGCTATTGTTGATCCAGTTCGCTTCAATGTAAGGCTGGAACTCACGCTGTTTACCATCGTCACGCt                                                                                                       >  5:281055/1‑149 (MQ=255)
                                           tCAACCCCGGTACGTACTTCACCGAGATTACGCGCACCTTCACGGCCTACGGTTTGACCATTCATCTTCACGGCGTAGACTTTGCTATTGTTGATCCAGTTCGCTTCAATGTAAGGCTGGAACTCACGCTGTTTACCATCGTCACGCTg                                                                                                      >  1:51648/1‑149 (MQ=255)
                                             acccccGGTACGTACTTCACCGAGATTACGCGCACCTTCACGGCCTACGGTTTGACCATTCATCTTCACGGCGTAGACTTTGCTATTGTTGATCCAGTTCGCTTCAATGTAAGGCTGGAACTCACGCTGTTTACCATCGTCACGCTGGt                                                                                                    <  1:272580/147‑1 (MQ=255)
                                             aaCCCCGGTACGTACTTCACCGAGATTACGCGCACCTTCACGGCCTACGGTTTGACCATTCATCTTCACGGCGTAGACTTTGCTATTGTTGATCCAGTTCGCTTCAATGTAAGGCTGGAACTCACGCTGTTTACCATCGTCACGCTg                                                                                                      <  8:105562/147‑1 (MQ=255)
                                             aaCCCCGGTACGTACTTCACCGAGATTACGCGCACCTTCACGGCCTACGGTTTGACCATTCATCTTCACGGCGTAGACTTTGCTATTGTTGATCCAGTTCGCTTCAATGTAAGGCTGGAACTCACGCTGTTTACCATCGTCACGCTGGt                                                                                                    <  2:51648/149‑1 (MQ=255)
                                              aCCCCGGTACGTACTTCACCGAGATTACGCGCACCTTCACGGCCTACGGTTTGACCATTCATCTTCACGGCGTAGACTTTGCTATTGTTGATCCAGTTCGCTTCAATGTAAGGCTGGAACTCACGCTGTTTACCATCGTCACGCTGGt                                                                                                    >  8:15921/1‑148 (MQ=255)
                                               ccccGGTACGTACTTCACCGAGATTACGCGCACCTTCACGGCCTACGGTTTGACCATTCATCTTCACGGCGTAGACTTTGCTATTGTTGATCCAGTTCGCTTCAATGTAAGGCTGGAACTCACGCTGTTTACCATCGTCACGCTGGTGa                                                                                                  <  5:266855/149‑1 (MQ=255)
                                                  cGGTACGTACTTCACCGAGATTACGCGCACCTTCACGGCCTACGGTGTGACCATTCATCTTCACGGCGTAGACTTTGCTATTGTTGATCCAGTTCGCTTCAATGTAAGGCTGGAACTCACGCTGTTTACCATCGTCACGCTGGTGATgg                                                                                               <  3:338064/149‑1 (MQ=255)
                                                     tacgtacTTCACCGAGATTACGCGCACCTTCACGGCCTACGGTTTGACCATTCATCTTCACGGCGTAGACTTTGCTATTGTTGATCCAGTTCGCTTCAATGTAAGGCTGGAACTCACGCTGTTTACCATCGTCACGCTGGTGATGGCTg                                                                                            >  8:341039/1‑149 (MQ=255)
                                                        gtacTTCACCGAGATTACGCGCACCTTCACGGCCTACGGTTTGACCATTCATCTTCACGGCGTAGACTTTGCTATTGTTGATCCAGTTCGCTTCAATGTAAGGCTGGAACTCACGCTGTTTACCATCGTCACGCTGGTGATGGCTGTTc                                                                                         <  8:180231/149‑1 (MQ=255)
                                                              cACCGAGATTACGCGCCCCTTCCCGGCCGACGGTTTGACCATTCCTCTTCACGGCGGAGACTTTGCTATTGTTGATCCAGTTCGCTTCAATGTAAGGCTGGAACTCACGCTGTTTACCATCGTCACGCTGGTGATGGCTGTTCaggt                                                                                     <  4:56256/147‑1 (MQ=255)
                                                               aCCGAGATTACGCGCACCTTCACGGCCTACGGTTTGACCATTCATCTTCACGGCGTAGACTTTGCTATTGTTGATCCAGTTCGCTTCAATGTAAGGCTGGAACTCACGCTGTTTACCATCGTCACGCTGGTGATGGCTGTTCaggtag                                                                                   <  7:81845/148‑1 (MQ=255)
                                                                 cGAGATTACGCGCACCTTCACGGCCTACGGTTTGACCATTCATCTTCACGGCGTAGACTTTGCTATTGTTGATCCAGTTCGCTTCAATGTAAGGCTGGAACTCACGCTGTTTACCATCGTCACGCTGGTGa                                                                                                  >  4:273695/1‑131 (MQ=255)
                                                                     aTTACGCGCTCCTTCACGGCCTACGGTTTGACCATTCATCTTCACGGCGTAGACTTTGCTATTGTTGATCCAGTTCGCTTCAATGTTAGGCTGGAACTCACGCTGTTTACCATCGTCACGCTGGTGATGGCTGTTCAGGTAGGTTTTc                                                                             <  6:236918/148‑1 (MQ=255)
                                                                     aTTACGCGCACCTTCACGGCCTACGGTTTGACCATTCATCTTCACGGCGTAGACTTTGCTATTGTTGATCCAGTTCGCTTCAATGTAAGGCTGGAACTCACGCTGTTTACCATCGTCACGCTGGTGATGGCTGTTCAGGTAGGTTTTCa                                                                            <  2:391662/149‑1 (MQ=255)
                                                                      ttACGCGCACCTTCCCGGCCTACGGTTTGCCCATTCATCTTCACGGCGTAGACGTTGCTATTGTTGATCCAGGTCGCTTCAATGTAAGGCTGGAACTCACGCTGTTTACCATCGTCACGCTGGTGATGGCTGTTCAGGTAGGTTTTCAc                                                                           >  2:151551/1‑149 (MQ=255)
                                                                      ttACGCGCACCTTCACGGCCTACGGTTTGACCATTCATCTTCACGGCGTAGACTTTGCTATTGTTGATCCAGTTCGCTTCAATGTAAGGCTGGAACTCACGCTGTTTACCATCGTCACGCTGGTGATGGCTGTTCAGGTAGGTTTTCAc                                                                           <  7:15921/149‑1 (MQ=255)
                                                                          gcgcACCTTCACGGCCTACGGTTTGACCATTCATCTTCACGGCGTAGACTTTGCTATTGTTGATCCAGTTCGCTTCAATGTAAGGCTGGAACTCACGCTGTTTACCATCGTCACGCTGGTGATGGCTGTTCAGGTAGGTTTTCAcccc                                                                        >  7:135166/1‑148 (MQ=255)
                                                                                cTTCACGGCCTACGGTTTGACCATTCATCTTCACGGCGTAGACTTTGCTATTGTTGATCCAGTTTGCTTCAATGTAAGGCTGGAACTCACGCTGTTTACCATCGTCACGCTGGTGATGGCTGTTCAGGTAGGTTTTCACCCCAAGTcgc                                                                 <  7:191823/149‑1 (MQ=255)
                                                                                    aCGGCCTACGGTTTGACCATTCATCTTCACGGCGTAGACTTTGCTATTGTTGATCCAGTTCGCTTCAATGTAAGGCTGGAACTCACGCTGTTTACCATCGTCACGCTGGTGATGGCTGTTCAGGTAGGTTTTCACCCCAAGTCGCGtt                                                              <  4:100086/148‑1 (MQ=255)
                                                                                         cTACGGTTTGACCATTCATCTTCACGGCGTAGACTTTGCTATTGTTGATCCAGTTCGCTTCAATGTAAGGCTGGAACTCACGCTGTTTACCATCGTCACGCTGGTGATGGCTGTTCAGGTAGGTTTTCACCCCAAGTCGCGTTTGcaca                                                        >  8:212156/1‑149 (MQ=255)
                                                                                          tACGGTTTGACCATTCATCTTCACGGCGTAGACTTTGCTATTGTTGATCCAGTTTGCTTCAATGTAAGGCTGGAACTCACGCTGTTTACCATCGTCACGCTGGTGATGGCTGTTCAGGTAGGTTTTCACCCCAAGTCGCGTTTGCACAt                                                       >  2:278114/1‑149 (MQ=255)
                                                                                              gTTTGACCATTCATCTTCACGGCGTAGACTTTGCTATTGTTGATCCAGTTCGCTTCAATGTAAGGCTGGAACTCACGCTGTTTACCATCGTCACGCTGGTGATGGCTGTTCAGGTAGGTTTTCACCCCAAGTCGCGTTTGCACATTTcc                                                   <  2:127800/149‑1 (MQ=255)
                                                                                               tttGACCATTCATCTTCACGGCGTAGACTTTGCTATTGTTGATCCAGTTCGCTTCAATGTAAGGCTGGAACTCACGCTGTTTACCATCGTCACGCTGGTGATGGCTGTTCAGGTAGGTTTTCACCCCAAGTCGCGTTTGCACATTTCCg                                                  >  7:7526/1‑149 (MQ=255)
                                                                                                  gACCATTCATCTTCACGGCGTAGACTTTGCTATTGTTGATCCAGTTCGCTTCAATGTAAGGCTGGAACTCACGCTGTTTACCATCGTCACGCTGGTGATGGCTGTTCAGGTAGGTTTTCACCCCAAGTCGCGTTTGCACATTTCCGTc                                                <  3:29699/148‑1 (MQ=255)
                                                                                                     cATTCATCTTCACGGCGTAGACTTTGCTATTGTTGATCCAGTTCGCTTCAATGTAAGGCTGGAACTCACGCTGTTTACCATCGTCACGCTGGTGATGGCTGTTCAGGTAGGTTTTCACCCCAAGTCGCGTTTGCACATTTCCGTCGCCt                                            >  1:242892/1‑149 (MQ=255)
                                                                                                     cATTCATCTTCACGGCGTAGACTTTGCTATTGTTGATCCAGTTCGCTTCAATGTAAGGCTGGAACTCACGCTGTTTACCATCGTCACGCTGGTGATGGCTGTTCAGGTAGGTTTTCACCCCAAGTCGCGTTTGCACATTTCCGTCGCCt                                            >  4:167530/1‑149 (MQ=255)
                                                                                                     cATTCATCTTCACGGCGTAGACTTTGCTATTGTTGATCCAGTTCGCTTCAAAGTAAGGCTGGAACTCACGCTCTTTACCATCGTCACGCTCGTGATGGCTGTTAAGGTAGGTTTTCACCCCAAGTCGCGTTTGCACATTTCCGACGCCt                                            >  2:43038/1‑149 (MQ=255)
                                                                                                          aTCTTCACGGCGTAGACTTTGCTATTGTTGATCCAGTTCGCTTCAATGTAAGGCTGGAACTCACGCTGTTTACCATCGTCACGCTGGTGATGGCTGTTCAGGTAGGTTTTCACCCCAAGTCGCGTTTGCACATTTCCGTCGCCTTCCGt                                       >  2:293861/1‑149 (MQ=255)
                                                                                                               cACGGCGTAGACTTTGCTATTGTTGATCCAGTTCGCTTCAATGTAAGGCTGGAACTCACGCTGTTTACCATCGTCACGCTGGTGATGGCTGTTCAGGTAGGTTTTCACCCCAAGTCGCGTTTGCACATTTCCGTCGCCTTCCGttt                                     <  8:135166/146‑1 (MQ=255)
                                                                                                               cACGGCGTAGACTTTGCTATTGTTGATCCAGTTCGCTTCAATGTAAGGCTGGAACTCACGCTGTTTACCATCGTCACGCTGGTGATGGCTGTTCAGGTAGGTTTTCACCCCAAGTCGCGTTTGCACATTTCCGTCGCCTTCCGttt                                     <  7:212156/146‑1 (MQ=255)
                                                                                                                 cGGCGTAGACTTTGCTATTGTTGATCCAGTTCGCTTCAATGTAAGGCTGGAACTCACGCTGTTTACCATCGTCACGCTGGTGATGGCTGTTCAGGTAGGTTTTCACCCCAAGTCGCGTTTGCACATTTCCGTCGCCTTCCGTTTCAATg                                <  7:341039/149‑1 (MQ=255)
                                                                                                                      tAGACTTTGCTATTGTTGATCCAGTTCGCTTCAATGTAAGGCTGGAACTCACGCTGTTTACCATCGTCACGCTGGTGATGGCTGTTCAGGTAGGTTTTCACCCCAAGTCGCGTTTGCACATTTCCGTCGCCTTCCGTTTCAATGCGCGt                           >  6:285754/1‑149 (MQ=255)
                                                                                                                           gttGCTATTGTTGATCCAGTTCGCTTCAATGTAAGGCTGGAACTCACGCTTTTTACCATCGTCACGCTGGTGATGGCTGTTCAGGTAGGTTTTCACCCCAAGTCGCGTTTGCACATTTCCGTCGCCTTCCGTTTCAATGCGCGTTCCGt                      >  7:286391/2‑149 (MQ=255)
                                                                                                                                 aTTGTTGATCCAGTTCGCTTCAATGTAAGGCTGGAACTCACGCTGTTTACCATCGTCACGCTGGTGATGGCTGTTCAGGTAGGTTTTCACCCCAAGTCGCGTTTGCACATTTCCGTCGCCTTCCGTTTCAATGCGCGTTCCGTCTTTcc                >  1:223819/1‑149 (MQ=255)
                                                                                                                                    gttgATCGAGTTCGCTTCAATGTAAGGCTGGAACTCTCGCTGTTTACCATCGTCACGCTGGTGATGGCTGTTCAGGTAGGTTTTCACCCCAAGTCGCGTTTGCACATTTCCGTCGCCTTCCGTTTCAATGCGCGTTCCGTCTTTCCgg              <  8:115741/148‑1 (MQ=255)
                                                                                                                                     ttgATCCAGTTTGCTTCAATGTAAGGCTGGAACTCACGCTGTTTACCATCGTCACGCTGGTGATGGCTGTTCAGGTAGGTTTTCACCCCAAGTCGCGTTTGCACATTTCCGTCGCCTTCCGTTTCAATGCGCGTTCCGTCTTTCCGGGt            <  1:278114/149‑1 (MQ=255)
                                                                                                                                       gATCCAGTTCGCTTCAATGTAAGGCTGGAACTCACGCTGTTTACCATCGTCACGCTGGTGATGGCTGTTCAGGTAGGTTTTCACCCCAAGTCGCGTTTGCACATTTCCGTCGCCTTCCGTTTCAATGCGCGTTCCGTCTTTCCGGGTa           <  5:285754/148‑1 (MQ=255)
                                                                                                                                        aTCCAGTTCGCTTCAATGTAAGGCTGGAACTCACGCTGTTTACCATCGTCACGCTGGTGATGGCTGTTCAGGTAGGTTTTCACCCCAAGTCGCGTTTGCACATTTCCGTCGCCTTCCGTTTc                                    >  7:320179/1‑122 (MQ=255)
                                                                                                                                        aTCCAGTTCGCTTCAATGTAAGGCTGGAACTCACGCTGTTTACCATCGTCACGCTGGTGATGGCTGTTCAGGTAGGTTTTCACCCCAAGTCGCGTTTGCACATTTCCGTCGCCTTCCGTTTc                                    <  8:320179/122‑1 (MQ=255)
                                                                                                                                        aTCCAGTTCGCTTCAATGTAAGGCTGGAACTCACGCTGTTTACCATCGTCACGCTGGTGATGGCTGTTCAGGTAGGTTTTCACCCCAAGTCGCGTTTGCACATTTCCGTCGCCTTCCGTTTCAATgcgc                             <  2:225714/129‑1 (MQ=255)
                                                                                                                                        aTCCAGTTCGCTTCAATGTAAGGCTGGAACTCACGCTGTTTACCATCGTCACGCTGGTGATGGCTGTTCAGGTAGGTTTTCACCCCAAGTCGCGTTTGCACATTTCCGTCGCCTTCCGTTTCAATgcgc                             >  1:225714/1‑129 (MQ=255)
                                                                                                                                         tCCAGTTCGCTTCAATGAAAGCCTGGAACTCACGCTGTTTACCATCGTCACGCTGGTGATGGCTGTTCAGGTAGGTTTTCACCCCAAGTCGCGTTTGCACATTTCCGTCGCCTTCAGTTTCAATTCGCGTTCCGTCTTTCCGGGTATgg        >  2:146197/1‑149 (MQ=255)
                                                                                                                                           cAGTTCGCTTCAATGTAAGGCTGGAACTCACGCTGTTTACCATCGTCACGCTGGTGATGGCTGTTCAGGTAGGTTTTCACCCCAAGTCGCGTTTGCACATTTCCGTCGCCTTCCGTTTCAATGCGCGTTCCGTCTTTCCGGGTATGGTc      <  6:281055/149‑1 (MQ=255)
                                                                                                                                              ttCGCTTCAATGTAAGGCTGGAACTCACGCTGTTTACCATCGTCACGCTGGTGATGGCTGTTCAGGTAGGTTTTCACCCCAAGTCGCGTTTGCACATTTCCGTCGCCTTCCGTTTCAATGCGCGTTCCGTCTTTCCGGGTATGGTCgg    >  6:367317/1‑148 (MQ=255)
                                                                                                                                                cGCTTCAATGTAAGGCTGGAACTCACGCTGTTTACCATCGTCACGCTTGTGATGGCTGTTCAGGTAGGTTTTCACCCCAAGTCGCGTTTGCACATTTCCGTCGCCTTCCGTTTCAATGCGCGTTCCGTCTTTCCGGGTATGGTCGGaa  <  1:34386/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
ACCGACATTCCCCCACAGGCTAAGGTTGTTATTTACTTTCGCCTCAACCCCGGTACGTACTTCACCGAGATTACGCGCACCTTCACGGCCTACGGTTTGACCATTCATCTTCACGGCGTAGACTTTGCTATTGTTGATCCAGTTCGCTTCAATGTAAGGCTGGAACTCACGCTGTTTACCATCGTCACGCTGGTGATGGCTGTTCAGGTAGGTTTTCACCCCAAGTCGCGTTTGCACATTTCCGTCGCCTTCCGTTTCAATGCGCGTTCCGTCTTTCCGGGTATGGTCGGAA  >  NC_000913/1592731‑1593022

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