Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,429,024 A→C 6.1% intergenic (+40/‑25) lomR → / → stfR pseudogene, Rac prophage lom homolog;Phage or Prophage Related; interrupted by IS5 and N‑ter deletion/Rac prophage; putative tail fiber protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,429,0240AC 6.1% 257.9 / 7.8 82intergenic (+40/‑25)lomR/stfRpseudogene, Rac prophage lom homolog;Phage or Prophage Related; interrupted by IS5 and N‑ter deletion/Rac prophage; putative tail fiber protein
Reads supporting (aligned to +/- strand):  ref base A (38/39);  new base C (2/3);  total (40/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.04e-01

CAGTTTAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGA  >  NC_000913/1428877‑1429170
                                                                                                                                                   |                                                                                                                                                  
cAGTTTAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAg                                                                                                                                                   >  2:283636/1‑149 (MQ=255)
cAGTTTAACCCGACCGAATCCGTGGCCATTGATAATGCTTATGAAGGCGCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAg                                                                                                                                                   <  5:245018/149‑1 (MQ=255)
 aGTTTAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAgg                                                                                                                                                  >  2:247123/1‑149 (MQ=255)
  gTTTAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGc                                                                                                                                                 >  5:86819/1‑149 (MQ=255)
   tttAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCg                                                                                                                                                >  7:110127/1‑149 (MQ=255)
   tttAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCg                                                                                                                                                <  5:200689/149‑1 (MQ=255)
     tAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCgg                                                                                                                                               <  5:211498/148‑1 (MQ=255)
      aaCCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGc                                                                                                                                             >  7:245489/1‑149 (MQ=255)
       aCCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCTGAACCGGCGGGCt                                                                                                                                            <  1:118889/149‑1 (MQ=255)
       aCCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCt                                                                                                                                            >  4:243051/1‑149 (MQ=255)
         ccgaccgaATCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCGGGTTCAGGCGGGCttt                                                                                                                                          >  4:120870/1‑149 (MQ=255)
            accgaATCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCtttttt                                                                                                                                       <  6:253879/149‑1 (MQ=255)
              cgaATCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTtgt                                                                                                                                     >  2:112540/1‑149 (MQ=255)
               gaATCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTtgtg                                                                                                                                    <  1:129579/149‑1 (MQ=255)
                aaTCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTgg                                                                                                                                   <  3:243051/149‑1 (MQ=255)
                aaTCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGGCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTgg                                                                                                                                   <  5:335819/149‑1 (MQ=255)
                 aTCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTggg                                                                                                                                  <  6:234197/149‑1 (MQ=255)
                  tCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCGGGTTCAGGCGGGCTTTTTTGTgggg                                                                                                                                 <  3:120870/149‑1 (MQ=255)
                     gTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTg                                                                                                                               >  5:167681/1‑148 (MQ=255)
                        gCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAata                                                                                                                           <  2:76944/149‑1 (MQ=255)
                        gCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGGCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGGTCAGGCGGGCTTTTTTGTGGGGTGAata                                                                                                                           <  1:112540/149‑1 (MQ=255)
                         ccATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAatat                                                                                                                          >  1:294428/1‑149 (MQ=255)
                          cATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTAGAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATg                                                                                                                         >  6:314621/1‑149 (MQ=255)
                           aTTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATgg                                                                                                                        >  6:215948/1‑149 (MQ=255)
                           aTTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATg                                                                                                                         >  2:154880/1‑148 (MQ=255)
                           aTTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATg                                                                                                                         <  1:154880/148‑1 (MQ=255)
                              gATATTGCTTATGAAGGCCCCGGCTGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCa                                                                                                                      <  2:318023/148‑1 (MQ=255)
                                  ttGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTCGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTaaa                                                                                                                 <  6:65927/149‑1 (MQ=255)
                                    gCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTa                                                                                                                   >  2:146719/1‑145 (MQ=255)
                                     cTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGAt                                                                                                              <  3:276117/149‑1 (MQ=255)
                                     cTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGAt                                                                                                              <  5:215948/149‑1 (MQ=255)
                                     cTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGAt                                                                                                              >  4:276117/1‑149 (MQ=255)
                                        aTGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGAttt                                                                                                            <  4:66629/148‑1 (MQ=255)
                                         tGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTc                                                                                                           <  7:50727/148‑1 (MQ=255)
                                         tGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTc                                                                                                           <  3:70804/148‑1 (MQ=255)
                                         tGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTc                                                                                                           <  8:110127/148‑1 (MQ=255)
                                           aaGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAgg                                                                                                        <  8:245489/149‑1 (MQ=255)
                                                tccGGCAGTGGCGACTGGCGTACTGACGGATTCATCGTTGGGGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGgagtcc                                                                                                   >  8:287313/2‑147 (MQ=255)
                                                cccGGCAGTGGCGCCTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTAc                                                                                                   >  1:229114/1‑149 (MQ=255)
                                                cccGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTAc                                                                                                   >  3:243774/1‑149 (MQ=255)
                                                cccGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGGGAATATGGCAGTaaa                                                                                                                 >  3:238684/1‑135 (MQ=255)
                                                cccGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGGGAATATGGCAGTaaa                                                                                                                 <  4:238684/135‑1 (MQ=255)
                                                          gCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACt                                                                                                  <  1:52551/140‑1 (MQ=255)
                                                          gCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACt                                                                                                  >  2:52551/1‑140 (MQ=255)
                                                              cTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGcaca                                                                                     <  8:105948/149‑1 (MQ=255)
                                                                 gcgcACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAg                                                                                    <  5:314621/147‑1 (MQ=255)
                                                                    cACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGaaaa                                                                               >  2:48727/1‑149 (MQ=255)
                                                                     aCTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAAc                                                                              <  6:86819/149‑1 (MQ=255)
                                                                                 aTCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACt                                                                  >  1:66577/1‑149 (MQ=255)
                                                                                  tCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATTGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTg                                                                 >  8:337483/1‑149 (MQ=255)
                                                                                     tGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGca                                                               >  4:208657/1‑148 (MQ=255)
                                                                                         tgtCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATc                                                          >  6:30559/1‑149 (MQ=255)
                                                                                           tCGGTTATAAGTTCGGATTACCCAGGGAACACAGTTTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATcc                                                         <  7:132957/148‑1 (MQ=255)
                                                                                             ggTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACCGCACCATTCAg                                                       <  7:287313/148‑1 (MQ=255)
                                                                                                  taAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGGTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGaaa                                                 <  3:221908/149‑1 (MQ=255)
                                                                                                   aaGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAg                                                >  3:67380/1‑149 (MQ=255)
                                                                                                      ttCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCaa                                             <  2:294428/149‑1 (MQ=255)
                                                                                                          gATTAGCCAGGTAACACAGTGTTATGAAAGCCCGCCTGAACCGGCGGGCtgt                                                                                                                                          <  6:265461/52‑3 (MQ=17)
                                                                                                          gATTAGCCAGGTAACACAGTGTTATGAAAGCCCGCCTGAACCGGCGGGCtgt                                                                                                                                          >  5:265461/1‑50 (MQ=17)
                                                                                                           aTTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACg                                         >  8:188110/1‑148 (MQ=255)
                                                                                                              aGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTaa                                      <  1:48727/148‑1 (MQ=255)
                                                                                                               gCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAAc                                     <  7:188110/148‑1 (MQ=255)
                                                                                                               gCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACa                                    <  3:208657/149‑1 (MQ=255)
                                                                                                                ccAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAg                                   >  6:117259/1‑149 (MQ=255)
                                                                                                                ccAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAg                                   >  8:271146/1‑149 (MQ=255)
                                                                                                                  aGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTAc                                                                        <  2:202453/110‑1 (MQ=255)
                                                                                                                  aGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTAc                                                                        >  1:202453/1‑110 (MQ=255)
                                                                                                                         acaGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACggt                          >  4:171218/1‑149 (MQ=255)
                                                                                                                         acaGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACggt                          >  6:265329/1‑149 (MQ=255)
                                                                                                                           aGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACggtgg                        <  2:66577/149‑1 (MQ=255)
                                                                                                                               ttATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACggtggtggt                    >  1:137919/1‑149 (MQ=255)
                                                                                                                               ttATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACggtggtggt                    >  6:115882/1‑149 (MQ=255)
                                                                                                                               gtATGACAGCCCGCCGGTTCAGGCGGTCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACATGAAATCCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACggtggtggt                    <  5:131777/148‑1 (MQ=255)
                                                                                                                                 aTGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGa                  >  5:86998/1‑149 (MQ=255)
                                                                                                                                  tGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGaa                 >  2:29404/1‑149 (MQ=255)
                                                                                                                                  tGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGaa                 >  2:43250/1‑149 (MQ=255)
                                                                                                                                   gACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGCAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGaa                 >  6:322504/1‑148 (MQ=255)
                                                                                                                                    aCAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTAGAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAaca               <  1:283636/149‑1 (MQ=255)
                                                                                                                                    aCAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAac                <  3:248283/148‑1 (MQ=255)
                                                                                                                                          cgccgGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAAc                                          <  8:342868/116‑1 (MQ=255)
                                                                                                                                                ttCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCt      >  7:34074/1‑146 (MQ=255)
                                                                                                                                                  cAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGa  <  1:146719/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CAGTTTAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGA  >  NC_000913/1428877‑1429170

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