Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,792,240 T→A 8.3% intergenic (‑220/+27) ydiV ← / ← nlpC anti‑FlhD4C2 factor, inactive EAL family phosphodiesterase/putative C40 clan peptidase lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,792,2400TA 8.3% 213.0 / 8.8 72intergenic (‑220/+27)ydiV/nlpCanti‑FlhD4C2 factor, inactive EAL family phosphodiesterase/putative C40 clan peptidase lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (31/35);  new base A (3/3);  total (34/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.37e-01

ACAGAAGGAGTAACTTTCATTTGTTCCATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGGTTGTTGGTATCGTAAATAC  >  NC_000913/1792093‑1792382
                                                                                                                                                   |                                                                                                                                              
aCAGAAGGTGTAACTTTCATTTGTTCCATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAATATACTACTGTATTCAGGATTAAATGGACGGAGATGt                                                                                                                                                <  2:193823/148‑1 (MQ=255)
   gAAGGAGTAACTTTCATTTGTTCCATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTTAGGAGTAAATGGACGGCGATGTatca                                                                                                                                            <  1:48398/149‑1 (MQ=255)
       gAGTAACTTTCATTTGTTCCATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTatcatc                                                                                                                                          >  6:340011/1‑147 (MQ=255)
       gAGTAACTTTCATTTGTTCCATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTatcatc                                                                                                                                          <  5:340011/147‑1 (MQ=255)
          tAACTTTCATTTGTTCCATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGt                                                                                                                                     >  5:64552/1‑149 (MQ=255)
           aaCTTTCATTTGTTCCATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTc                                                                                                                                    >  2:110057/1‑149 (MQ=255)
              tttCATTTGTTCCATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCAt                                                                                                                                 >  5:197632/1‑149 (MQ=255)
                  aTTTGTTCCATGTTAACCACTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTc                                                                                                                            <  4:39498/149‑1 (MQ=255)
                   tttGTGCCATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTc                                                                                                                            <  5:7027/149‑1 (MQ=255)
                      gTTCCATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGa                                                                                                                          <  3:293159/148‑1 (MQ=255)
                         ccATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGATACATCGCCGt                                                                                                                                     >  6:97565/1‑134 (MQ=255)
                         ccATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGATACATCGCCGt                                                                                                                                     <  5:97565/134‑1 (MQ=255)
                          cATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATAc                                                                                                                     >  2:283308/1‑149 (MQ=255)
                                 aCCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCAttttt                                                                                                                             >  8:49244/1‑134 (MQ=255)
                                 aCCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCAttttt                                                                                                                             <  7:49244/134‑1 (MQ=255)
                                 aCCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTc                                                                                                              <  2:154589/149‑1 (MQ=255)
                                 aCCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTc                                                                                                              <  7:142058/149‑1 (MQ=255)
                                  ccACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCg                                                                                                             >  7:15316/1‑149 (MQ=255)
                                   cACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGt                                                                                                            <  4:96880/149‑1 (MQ=255)
                                      tttttCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGtctc                                                                                                         >  7:273676/1‑149 (MQ=255)
                                           cAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCt                                                                                                     <  1:110057/148‑1 (MQ=255)
                                            aGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGc                                                                                                   >  2:180363/1‑149 (MQ=255)
                                              gggTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCa                                                                                                 <  3:301110/149‑1 (MQ=255)
                                                gTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAg                                                                                                >  1:93863/1‑148 (MQ=255)
                                                 tCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATAGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGa                                                                                               >  1:331866/1‑148 (MQ=255)
                                                   cTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGaaa                                                                                             <  7:296933/148‑1 (MQ=255)
                                                   cTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAAt                                                                                            >  2:350825/1‑149 (MQ=255)
                                                   cTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAAt                                                                                            <  1:203690/149‑1 (MQ=255)
                                                    tttttttAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAAtt                                                                                           >  3:243982/1‑149 (MQ=255)
                                                      tttttAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAAtttt                                                                                         >  4:190895/1‑149 (MQ=255)
                                                       ttttAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAAtttt                                                                                         <  4:66311/148‑1 (MQ=255)
                                                            gATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGc                                                                                    <  7:315462/148‑1 (MQ=255)
                                                            aaTTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGATACATCGCCg                                                                                                                                      >  8:278512/2‑98 (MQ=39)
                                                            aaTTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGATACATCGCCg                                                                                                                                      <  7:278512/97‑1 (MQ=39)
                                                              ttATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCaa                                                                                 <  2:93863/149‑1 (MQ=255)
                                                                aTCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCaa                                                                                 >  7:294156/1‑147 (MQ=255)
                                                                 tCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTTCTGTATTCAGGAGTAAATGGACGGCGATGTAACATCGCCGTCCATTTTTCGACATACTAGATTCGTCGCGCCTGCCAGAAATTTTTCAGCCaa                                                                                 >  6:325974/1‑146 (MQ=255)
                                                                      aTTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACAtt                                                                         <  6:64552/149‑1 (MQ=255)
                                                                         ataAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTAt                                                                       <  8:347024/148‑1 (MQ=255)
                                                                               gAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGt                                                                  >  1:78298/1‑147 (MQ=255)
                                                                               gAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGa                                                                <  2:78298/149‑1 (MQ=255)
                                                                                aataatCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAg                                                               >  3:211929/1‑149 (MQ=255)
                                                                                    atCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAAc                                                           <  8:338317/149‑1 (MQ=255)
                                                                                     tCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACg                                                          >  7:16789/1‑149 (MQ=255)
                                                                                     tCCATTTACGGAATTTTTGTCTGCAAAATACTACGGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATAAAGTGAGGAACg                                                          >  1:64123/1‑149 (MQ=255)
                                                                                      ccATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACg                                                          >  6:10411/1‑148 (MQ=255)
                                                                                        aTTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGtctc                                                                                                         >  8:123075/1‑99 (MQ=255)
                                                                                        aTTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGtctc                                                                                                         <  7:123075/99‑1 (MQ=255)
                                                                                            aCGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGc                                                         <  1:19407/143‑1 (MQ=255)
                                                                                            aCGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGc                                                         >  2:19407/1‑143 (MQ=255)
                                                                                              ggAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCTTCTCGCCTGACAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAAAGCATACCTcc                                                 >  5:41440/1‑149 (MQ=255)
                                                                                              ggAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTcc                                                 <  6:197632/149‑1 (MQ=255)
                                                                                              ggAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTcc                                                 <  5:10411/149‑1 (MQ=255)
                                                                                               gAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTccc                                                <  5:183366/149‑1 (MQ=255)
                                                                                                   ttttGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACAt                                                                          >  1:207942/1‑119 (MQ=255)
                                                                                                   ttttGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACAt                                                                          <  2:207942/119‑1 (MQ=255)
                                                                                                      tGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTg                                          <  8:15316/148‑1 (MQ=255)
                                                                                                            cAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAgg                                   >  1:350903/1‑149 (MQ=255)
                                                                                                               aaTACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGt                                <  1:180363/149‑1 (MQ=255)
                                                                                                                   ctactGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGAt                            <  3:190895/149‑1 (MQ=255)
                                                                                                                           ttCAGGAGTAAATGGAGGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGAAGAGGCGTGGATAAATTgg                     <  4:243982/148‑1 (MQ=255)
                                                                                                                              aGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTgg                     >  2:20376/1‑145 (MQ=255)
                                                                                                                              aGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTgg                     <  1:20376/145‑1 (MQ=255)
                                                                                                                              aGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGgttgt                 >  7:216988/1‑149 (MQ=255)
                                                                                                                              aGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCGTGCTGGTAGAGGCGTGGATAAATTGgttgt                 >  3:205603/1‑149 (MQ=255)
                                                                                                                               ggAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGgttgtt                <  7:68068/149‑1 (MQ=255)
                                                                                                                                gAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGGttgttg               <  8:294156/149‑1 (MQ=255)
                                                                                                                                     aaTGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGGTTGTTGGTATc          <  2:64123/149‑1 (MQ=255)
                                                                                                                                     aaTGGACGGCGATGTATCATCGCCGGCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGGTTGTTGGTATc          <  2:350903/149‑1 (MQ=255)
                                                                                                                                         gACGGCGATGAATCATCGCCGTCCATTTTTCGCCATACTAGCTTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAATTGAGGAACGCCGAACTCCCTTGCTGGTAGAGGCGTGGATAAATTGGTTGTTTGTATCgga       >  2:266874/1‑146 (MQ=255)
                                                                                                                                          aCGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCaa                                                                    >  8:203136/1‑86 (MQ=255)
                                                                                                                                           cGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGGTTGTTGGTATCGTaaa     >  3:16209/1‑148 (MQ=255)
                                                                                                                                             gagATGAATCATCGCCGTCCCTTTTTCGACATACTAGATTCGTATCGGCTTCCAGAAATTTTTCTGCCAATAGACATTATGAAGTTAGGAACGCATCACTCCCTTGCTGGTCGCGGCGTGGATAAATTGGTTGTTGGGATCGTAAATAc  <  1:266874/147‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
ACAGAAGGAGTAACTTTCATTTGTTCCATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGGTTGTTGGTATCGTAAATAC  >  NC_000913/1792093‑1792382

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