Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,635,293 G→A 7.5% T185T (ACC→ACT stfQ ← Qin prophage; putative side tail fibre assembly protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,635,2930GA 7.5% 277.6 / 9.5 80T185T (ACC→ACTstfQQin prophage; putative side tail fibre assembly protein
Reads supporting (aligned to +/- strand):  ref base G (37/37);  new base A (3/3);  total (40/40)
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.30e-01

CCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTGGATGCGCTGGCGCTGTGGGTGTGCGATTTAATGCCGTCCTGTTCCTGAGATAATACGGCCCGACCACTGGCGGGCTTGCCCTTAA  >  NC_000913/1635168‑1635426
                                                                                                                             |                                                                                                                                     
ccATTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCaaa                                                                                                                         <  2:104469/140‑1 (MQ=255)
   ttAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAacga                                                                                                                      >  2:51308/1‑140 (MQ=255)
     aGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAacgact                                                                                                                    >  2:944140/1‑139 (MQ=255)
     aGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAacgacg                                                                                                                    >  1:616877/1‑140 (MQ=255)
     aGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAacgacg                                                                                                                    >  2:1671438/1‑140 (MQ=255)
     aGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAacgacg                                                                                                                    >  2:1671844/1‑140 (MQ=255)
     aGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAacgacg                                                                                                                    >  2:1685392/1‑140 (MQ=255)
           aaTGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTgg                                                                                                              >  2:775256/1‑140 (MQ=255)
                          gCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAAt                                                                                               <  2:931743/140‑1 (MQ=255)
                              ccaAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGt                                                                                           >  2:1610667/1‑140 (MQ=255)
                               caAATGCACCGGACCTCTTGTGTTGGTGCGCCCCGGCACTATTGGCGGACCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTACTCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTa                                                                                          <  2:1738454/140‑1 (MQ=37)
                                    gCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTgta                                                                                     <  2:548112/140‑3 (MQ=255)
                                    gCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTGGa                                                                                     <  2:545417/140‑1 (MQ=255)
                                        cGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTGGATgcg                                                                                 >  1:15485/1‑140 (MQ=255)
                                        cGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTGGATgcg                                                                                 >  1:15568/1‑140 (MQ=255)
                                        cGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTGGATgcg                                                                                 >  1:12706/1‑140 (MQ=255)
                                        cGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGAATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCAGTACTGGATgcg                                                                                 >  1:992105/1‑140 (MQ=255)
                                          gAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTGGATgcgct                                                                               <  2:434763/140‑1 (MQ=255)
                                          gAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTGGATgcgct                                                                               <  2:432099/140‑1 (MQ=255)
                                          gAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTGGATgcgct                                                                               <  1:1670076/140‑1 (MQ=255)
                                          gAACTCTTGTGTTGGTGCGCACCGGCACTATTGGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTGGATgcgct                                                                               <  1:1524421/140‑1 (MQ=255)
                                                    gtTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTGGATGCGCTGGCGCTGTgg                                                                     >  2:389095/1‑140 (MQ=255)
                                                     ttGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTGGATGCGCTGGCGCTGTggg                                                                    >  1:364651/1‑140 (MQ=255)
                                                         tGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTGGATGCGCTGGCGCTGTGGGtgtg                                                                <  1:1249591/140‑1 (MQ=255)
                                                         tGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTGGATGCGCTGGCGCTGTGGGtgtg                                                                <  1:1847574/140‑1 (MQ=255)
                                                              aCCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTGGATGCGCTGGCGCTGTGGGTGTGCGAtt                                                           <  1:1610667/140‑1 (MQ=255)
                                                                   cACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTGGATGCGCTGGCGCTGTGGGTGTGCGATTTAATg                                                      <  1:35823/140‑1 (MQ=255)
                                                                   cACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTGGATGCGCTGGCGCTGTGGGTGTGCGATTTAATg                                                      <  1:233662/140‑1 (MQ=255)
                                                                       aTTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTGGATGCGCTGGCGCTGTGGGTGTGCGATTTAATGCCGt                                                  >  1:1678023/1‑140 (MQ=255)
                                                                        tttGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTGGATGCGCTGGCGCTGTGGGTGTGCGATTTAATGCCGTc                                                 <  2:1755094/140‑1 (MQ=255)
                                                                        tttGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTGGATGCGCTGGCGCTGTGGGTGTGCGATTTAATGCCGTc                                                 >  2:821073/1‑140 (MQ=255)
                                                                        tttGCGGTCCCGCTAATACTATGGGAATGCCCCCCGGTGTAAATCGTGGATTTGGTTCCGTAAGCAAACGACGAAGAGGATTCCGTCCCCAAAGACCTCCTGGAGGCGCTGGCGGGCGCTGTGTGCGACTTAATgaagac                                                 >  2:1907573/1‑135 (MQ=255)
                                                                          tGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTAGTGCCGTAATCAAACGACGATGTGGTTTTCGTCCCCAAATCTGTACTGGATGCGCTGGCGCTGTGGGTATGCGATTTAATGCCGTCCt                                               >  2:785802/1‑140 (MQ=11)
                                                                             ggTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTGGATGCGCTGGCGCTGTGGGTGTGCGATTTAATGCCGTCCTGtt                                            >  2:1815208/1‑140 (MQ=255)
                                                                             ggTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTAGTGCCGTAATCAAACGACGATGTGGTTTTCGTCCCCAAATCTGTACTGGATGCGCTGGCGCTGTGGGTATGCGATTTAATGCCGTCCTGtt                                            <  1:1564964/140‑1 (MQ=11)
                                                                             ggTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTAGTGCCGTAATAAAACGACGATGTGGTTTTCGTCCCCAAATCTGTACTGGATGCGCTGGCGCTGTGGGTATGCGATTTAATGCCGTCCTGtt                                            <  1:1565056/140‑1 (MQ=11)
                                                                              gTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTGGATGCGCTGGCGCTGTGGGTGTGCGATTTAATGCCGTCCTGTTc                                           <  2:489645/140‑1 (MQ=255)
                                                                               tCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTGGATGCGCTGGCGCTGTGGGTGTGCGATTTAATGCCGTCCTGTTcc                                          <  1:296472/140‑1 (MQ=255)
                                                                               tCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTGGATGCGCTGGCGCTGTGGGTGTGCGATTTAATGCCGTCCTGTTcc                                          <  1:1663447/140‑1 (MQ=255)
                                                                               tCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTGGATGCGCTGGCGCTGTGGGTGTGCGATTTAATGCCGTCCTGTTcc                                          <  1:1663669/140‑1 (MQ=255)
                                                                               tCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTGGATGCGCTGGCGCTGTGGGTGTGCGATTTAATGCCGTCCTGTTcc                                          <  1:1663427/140‑1 (MQ=255)
                                                                               tCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTGGATGCGCTGGCGCTGTGGGTGTGCGATTTAATGCCGTCCTGTTcc                                          <  1:1663037/140‑1 (MQ=255)
                                                                               tCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTGGATGCGCTGGCGCTGTGGGTGTGCGATTTAATGCCGTCCTGTTcc                                          <  1:1865844/140‑1 (MQ=255)
                                                                               tCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTGGATGCGCTGGCGCTGTGGGTGTGCGATTTAATGCCGTCCTGTTcc                                          <  1:1049918/140‑1 (MQ=255)
                                                                               tCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTGGATGCGCTGGCGCTGTGGGTGTGCGATTTAATGCCGTCCTGTTcc                                          <  1:409898/140‑1 (MQ=255)
                                                                               tCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGGTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTGGATGCGCTGGCGCTGTGGGTGTGCGATTTAATGCCGTCCTGTTcc                                          <  2:617159/140‑1 (MQ=255)
                                                                                cccGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTGGATGCGCTGGCGCTGTGGGTGTGCGATTTAATGCCGTCCTGTTCCt                                         >  1:80818/1‑140 (MQ=255)
                                                                                cccGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTGGATGCGCTGGCGCTGTGGGTGTGCGATTTAATGCCGTCCTGTTCCt                                         >  1:80742/1‑140 (MQ=255)
                                                                                cccGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTGGATGCGCTGGCGCTGTGGGTGTGCGATTTAATGCCGTCCTGTTCCt                                         >  2:1390346/1‑140 (MQ=255)
                                                                                    ccaaTACTCGGGGTATGCGCCCCGGAGTTATTCGTGGATTTGGTTCCTACATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTGGATGCGCTGGCGCTGTGGGTGTGCGAATTAATGCCGTCCTGTTCCTgata                                     <  2:418415/138‑3 (MQ=17)
                                                                                      aaTACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTAGTGCCGTAATCAAACGACGATGTGGTTTTCGTccc                                                                                                    <  1:1835942/75‑1 (MQ=11)
                                                                                      aaTACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTAGTGCCGTAATCAAACGACGATGTGGTTTTCGTccc                                                                                                    >  2:1835942/1‑75 (MQ=11)
                                                                                         aCTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTGGATGCGCTGGCGCTGTGGGTGTGCGATTTAATGCCGTCCTGTTCCTGAGataata                                >  2:539214/1‑140 (MQ=255)
                                                                                             tGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTGGATGCGCTGGCGCTGTGGGTGTGCGATTTAATGCCGTCCTGTTCCTGAGATAATACGGc                            >  2:784172/1‑140 (MQ=255)
                                                                                               ggTATGCGCCCCGGTGTTATTCGTGCCTTTGGTTCCGTACTCAACCGACGATGTGGTTGCCGTCCCCAAATCCGTACTGGATGCGCTGGCGCTGTGGGTGTGCGATTAAATGCCGTCCGGGTCCTGAGATAATACGGccc                          <  2:187754/140‑1 (MQ=17)
                                                                                                   tGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTGGATGCGCTGGCGCTGTGGGTGTGCGATTTAATGCCGTCCTGTTCCTGAGATCATACGGCCCGAcc                      >  1:1104997/1‑140 (MQ=255)
                                                                                                   tGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTGGATGCGCTGGCGCTGTGGGTGTGCGATTTAATGCCGTCCTGTTCCTGAGATAATACGGCCCGAcc                      >  1:1920989/1‑140 (MQ=255)
                                                                                                   tGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTGGATGCGCTGGCGCTGTGGGTGTGCGATTTAATGCCGTCCTGTTCCTGAGATAATACGGCCCGAcc                      >  1:172519/1‑140 (MQ=255)
                                                                                                   tGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTGGATGCGCTGGCGCTGTGGGTGTGCGATTTAATGCCGTCCTGTTCCTGAGATAATACGGCCCGAcc                      >  1:1799724/1‑140 (MQ=255)
                                                                                                   tGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTGGATGCGCTGGCGCTGTGGGTGTGCGATTTAATGCCGTCCTGTTCCTGAGATAATACGGCCCGAcc                      >  1:6127/1‑140 (MQ=255)
                                                                                                   tGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTGGATGCGCTGGCGCTGTGGGTGTGCGATTTAATGCCGTCCTGTTCCTGAGATAATACGGCCCGAcc                      >  1:552560/1‑140 (MQ=255)
                                                                                                     tgcCCCGGTGTTATTCGTGGATTTAGTGCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTGGATGCGCTGGCGCTGTGGGTGTGCGATTTAATGCCGTCCTGTTCCTGAGATAATACGGCCCGACCAc                    >  2:804429/2‑140 (MQ=16)
                                                                                                     cgcCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTGGATGCGCTGGCGCTGTGGGTGTGCGATTTAATGCCGTCCTGTTCCTGAGATAATACGGCCCGACCAc                    >  1:1902367/1‑140 (MQ=255)
                                                                                                        cccGGTGTTATTGGTGGTTTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATACGTACTGGATGCGCTGGCGCTGTGGGTGTGCGATTTAATGCCGTCCTGTTCCTGAGATAATACGGCCCGACCACTgg                 <  1:1908928/140‑1 (MQ=255)
                                                                                                        cccGGTGTTATTCGTGTATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCAGTACTGGATGCGCTGGCGCTGTGGGTGTGCGATTTAATGCCGTCCTGTTCCTGAGATAATACGGCCCGACCACTgg                 <  1:138786/140‑1 (MQ=255)
                                                                                                        cccGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATACGTACTGGATGCGCTGGCGCTGTGGGTGTGCGATTTAATGCCGTCCTGTTCCTGAGATAATACGGCCCGACCACTgg                 <  1:1908967/140‑1 (MQ=255)
                                                                                                         ccGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTGGATGCGCTGGCGCTGTGGGTGTGCGATTTAATGCCGTCCTGTTCCTGAGATAATACGGCCCGACCACTGGc                >  2:1913473/1‑140 (MQ=255)
                                                                                                          cGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTGGATGCGCTGGCGCTGTGGGTGTGCGATTTAATGCCGTCCTGTTCCTGAGATAATACGGCCCGACCACTGGCg               >  2:392812/1‑140 (MQ=255)
                                                                                                          cGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTGGATGCGCTGGCGCTGTGGGTGTGCGATTTAATGCCGTCCTGTTCCTGAGATAATACGGCCCGACCACTGGCg               >  1:465464/1‑140 (MQ=255)
                                                                                                             tgtTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTGGATGCGCTGGCGCTGTGGGTGTGCGATTTAATGCCGTCCTGTTCCTGAGATAATACGGCCCGACCACTGGCGGGc            <  1:1671438/140‑1 (MQ=255)
                                                                                                             tgtTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTGGATGCGCTGGCGCTGTGGGTGTGCGATTTAATGCCGTCCTGTTCCTGAGATAATACGGCCCGACCACTGGCGGGc            <  1:1671844/140‑1 (MQ=255)
                                                                                                               ttATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTGGATGCGCTGGCGCTGTGGGTGTGCGATTTAATGCCGTCCTGTTCCTGAGATAATACGGCCCGACCACTGGCGGGCtt          <  2:900015/140‑1 (MQ=255)
                                                                                                               ttATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTGGATGCGCTGGCGCTGTGGGTGTGCGATTTAATGCCGTCCTGTTCCTGAGATAATACGGCCCGACCACTGGCGGGCtt          <  2:910748/140‑1 (MQ=255)
                                                                                                               ttATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTGGATGCGCTGGCGCTGTGGGTGTGCGATTTAATGCCGTCCTGTTCCTGAGATAATACGGCCCGACCACTGGCGGGCtt          <  2:900415/140‑1 (MQ=255)
                                                                                                                  ttCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTGGATGCGCTGGCGCTGTGGGTGTGCGATTTAATGCCGTCCTGTTCCTGAGATAATACGGCCCGACCACTGGCGGGCTTGcc       >  1:1078990/1‑140 (MQ=255)
                                                                                                                      tGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTGGATGCGCTGGCGCTGTGGGTGTGCGATTTAATGCCGTCCTGTTCCTGAGATAATACGGCCCGACCACTGGCGGGCTTGCCCTTa   >  1:593783/1‑140 (MQ=255)
                                                                                                                       ggATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTGGATGCGCTGGCGCTGTGGGTGTGCGATTTAATGCCGTCCTGTTCCTGAGATAATACGGCCCGACCACTGGCGGGCTTGCCCTTac  <  1:788701/140‑2 (MQ=255)
                                                                                                                       ggATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTGGATGCGCTGGCGCTGTGGGTGTGCGATTTAATGCCGTCCTGTTCCTGAGATAATACGGCCCGACCACTGGCGGGCTTGCCCTTaa  <  2:364651/140‑1 (MQ=255)
                                                                                                                       ggATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTGGATGCGCTGGCGCTGTGGGTGTGCGATTTAATGCCGTCCTGTTCCTGAGATAATACGGCCCGACCACTGGCGGGCTTGCCCTTaa  >  2:1371964/1‑140 (MQ=255)
                                                                                                                       ggATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTGGATGCGCTGGCGCTGTGGGTGTGCGATTTAATGCCGTCCTGTTCCTGAGATAATACGGCCCGACCACTGGCGGGCTTGCCCTTaa  >  1:1469657/1‑140 (MQ=255)
                                                                                                                             |                                                                                                                                     
CCATTAGGGAAAATGCTCGTGTTCGTGCCACCAAATGCACCGGAACTCTTGTGTTGGTGCGCACCGGCACTATTTGCGGTCCCGCTAATACTATGGGTATGCGCCCCGGTGTTATTCGTGGATTTGGTTCCGTAATCAAACGACGATGTGGTTTCCGTCCCCAAATCCGTACTGGATGCGCTGGCGCTGTGGGTGTGCGATTTAATGCCGTCCTGTTCCTGAGATAATACGGCCCGACCACTGGCGGGCTTGCCCTTAA  >  NC_000913/1635168‑1635426

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