Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,632,949 C→T A20A (GCC→GCT ynfT → Qin prophage; protein YnfT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,632,9490CT99.1% 247.9 / NA 108A20A (GCC→GCTynfTQin prophage; protein YnfT
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (64/43);  total (65/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 108 total aligned reads displayed.

AAAAATTGTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCCCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCAGACATACAATATCGCAATGAAAATCAATA  >  NC_000913/1632841‑1633051
                                                                                                            |                                                                                                      
aaaaaTTGTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCCGCATTTTTAGTCATGCCCTAATAAAATTAGGTATTGGCATATCCAATa                                                                         >  2:1097627/1‑140 (MQ=255)
       gTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGa                                                                  <  2:71012/140‑1 (MQ=14)
       gTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGa                                                                  >  1:508839/1‑140 (MQ=17)
               tatCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCg                                                          <  1:31299/140‑1 (MQ=14)
                  cATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATg                                                       <  1:1354661/140‑1 (MQ=17)
                  cATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATg                                                       <  1:1729758/140‑1 (MQ=14)
                    tatataACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGta                                                     <  1:2031858/140‑1 (MQ=14)
                      tataACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGtata                                                   >  1:1744741/1‑140 (MQ=18)
                      tataACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGtata                                                   >  1:1606920/1‑140 (MQ=18)
                      tataACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATGGACATATCCAATTACCCCGACAAAACAGATGtata                                                   >  2:1175244/1‑140 (MQ=11)
                       ataACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATaa                                                  >  1:81719/1‑140 (MQ=18)
                        taACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAg                                                 <  1:2349190/140‑1 (MQ=14)
                        taACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAg                                                 <  2:187825/140‑1 (MQ=14)
                         aaCACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGt                                                >  1:320018/1‑140 (MQ=18)
                         aaCACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGt                                                <  1:1405518/140‑1 (MQ=14)
                           cacaTTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGt                                                >  1:1637987/1‑138 (MQ=17)
                           cacaTTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGt                                                <  2:1637987/138‑1 (MQ=14)
                           cacaTTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCa                                              >  1:2059325/1‑140 (MQ=18)
                           cacaTTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCa                                              >  1:2253949/1‑140 (MQ=18)
                           cacaTTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCa                                              >  2:1347816/1‑140 (MQ=18)
                             caTTACTAATTCGAGGCTATATCAACAGCATACTGCTAAACACATCTCGCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAAc                                            <  2:884880/140‑1 (MQ=11)
                              aTTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATCTCCAATAGCCCTGAGAAAACCAATGTATAAGTCAACa                                           >  1:697569/1‑140 (MQ=11)
                              aTTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATTCAAAAACCCAGACAAAAACGATGTATAAGTCAACa                                           >  2:1029031/1‑140 (MQ=255)
                              aTTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAg                                                 <  2:1336573/134‑1 (MQ=17)
                              aTTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACa                                           >  1:35034/1‑140 (MQ=17)
                              aTTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACa                                           >  2:2251609/1‑140 (MQ=18)
                              aTTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAAACAAAACCGATGTATAAg                                                 >  1:1336573/1‑134 (MQ=17)
                              aTTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAAGGATTGGCATATACAAAAACCCAGAAAAAATCGATGTATCAGTCAACa                                           >  2:1175903/1‑140 (MQ=11)
                               ttACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAAGCCAGACGAAACCGGTGTATAAGTCAACat                                          >  1:884880/1‑140 (MQ=255)
                                tACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGTATTGGCATCTCCAATAACCCAGACAAAACCAATGTATAAGTCAACata                                         >  2:906891/1‑140 (MQ=255)
                                tACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATTCAATAACACAGTCACAACCGAGGTATAAGTCAACata                                         >  1:758552/1‑140 (MQ=11)
                                tACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGATAAAACCGATGta                                                     <  1:1825682/128‑1 (MQ=9)
                                tACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGATAAAACCGATGta                                                     >  2:1825682/1‑128 (MQ=21)
                                tACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGta                                                     >  2:1313023/1‑128 (MQ=21)
                                tACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGta                                                     <  1:1313023/128‑1 (MQ=14)
                                tACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACata                                         >  2:143875/1‑140 (MQ=18)
                                tACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACata                                         >  1:1360039/1‑140 (MQ=18)
                                tACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACata                                         >  2:1904998/1‑140 (MQ=18)
                                tACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACata                                         >  1:2424649/1‑140 (MQ=17)
                                tACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACata                                         >  1:1546861/1‑140 (MQ=18)
                                tACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGAAGTATAAGTCAACata                                         >  1:888385/1‑140 (MQ=17)
                                tACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACACAGTCAAAACCGATGTATAAGTCAACata                                         >  1:1093377/1‑140 (MQ=255)
                                tACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCATTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACata                                         >  2:2423072/1‑140 (MQ=16)
                                  cTAATTCGAGGCTATATGAACAGCATACTGATAATCATATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATc                                       <  1:1175903/140‑1 (MQ=17)
                                  cTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATc                                       <  1:143875/140‑1 (MQ=17)
                                     aTTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGTCATATCCAATAAACCAGACAACACCAATGTGTAAGTCAACATATCCTg                                    >  2:730599/1‑140 (MQ=11)
                                        cGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACa                                           >  2:2061759/1‑130 (MQ=17)
                                        cGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACa                                           <  1:2061759/130‑1 (MQ=17)
                                        cGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAAt                                 <  1:2251609/140‑1 (MQ=17)
                                        cGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAAt                                 <  2:1360039/140‑1 (MQ=17)
                                        cGAGGCTAGATGAACATCATACTGATAATCACACCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAAt                                 <  2:1209878/140‑1 (MQ=255)
                                         gAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATc                                <  2:1410805/140‑1 (MQ=17)
                                         gAGGCTATATGAACAGCATACTCATCATCGCATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATc                                <  1:1175244/140‑1 (MQ=255)
                                          aGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAAt                                 <  1:101666/138‑1 (MQ=17)
                                          aGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAAt                                 >  2:101666/1‑138 (MQ=18)
                                          aGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCa                               >  1:332506/1‑140 (MQ=18)
                                          aGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCa                               >  2:2439577/1‑140 (MQ=18)
                                          aGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCa                               >  1:1082991/1‑140 (MQ=18)
                                          aGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATAACCAATAACTCAGACAAAACCGATGTATAAGGCAACATATACTGAATCa                               >  2:1192517/1‑140 (MQ=11)
                                           ggCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCAg                              <  2:2059325/140‑1 (MQ=17)
                                           ggCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCAg                              <  2:1324559/140‑1 (MQ=17)
                                           ggCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCAg                              <  1:1347816/140‑1 (MQ=17)
                                            gCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGAATGAGTCAAGATATCCTGAATCAGa                             >  1:788827/1‑140 (MQ=255)
                                            gCTATATGAACAGCAGACTGATAATCACAACTCTCCTGATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCAGa                             <  1:1029031/140‑1 (MQ=255)
                                              tataTGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCAGACa                           <  1:1220131/140‑1 (MQ=17)
                                                tatGAACAGCATACTGATAATCTCATCTCTCCTTATCATTTTCAGCATGTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCAGACATa                         <  2:1093377/140‑1 (MQ=255)
                                                tatGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAATACCGATGGATAAGTGAACATATCCTGAATCAGACATa                         >  1:1209878/1‑140 (MQ=255)
                                                tatGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCAGACATa                         >  2:207477/1‑140 (MQ=18)
                                                tatGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCAGACATa                         >  1:2412144/1‑140 (MQ=18)
                                                      cAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCAGACATACAatat                   >  1:65372/1‑140 (MQ=18)
                                                       aGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAAt                                 <  1:2153244/125‑1 (MQ=17)
                                                       aGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAAt                                 >  2:2153244/1‑125 (MQ=18)
                                                        gCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCAGACATACAATATCg                 >  1:2409340/1‑140 (MQ=17)
                                                         cATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTCTAAGTCAACATATCATGAATCAGACATACAATATCGc                >  1:871054/1‑140 (MQ=9)
                                                         cATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCAGACATACAATATCGc                >  2:1674645/1‑140 (MQ=18)
                                                            aCTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCAGACATACAATATCGCAAt             >  2:414569/1‑140 (MQ=17)
                                                             cTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCAGACATACAATATCGCAATg            <  2:65372/140‑1 (MQ=17)
                                                              tGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACGGATGTATAAGTCCACATATCCTGAAGCAGACTTACAATATCGCAATGa           >  2:1160318/1‑140 (MQ=11)
                                                                aTAATCACATCTCTCGTTATCAGATGCAGTATTTTTAATCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCAGACATACAATATCGCAATGaaa         <  1:1160318/140‑1 (MQ=11)
                                                                aTAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCAGACa                           >  2:2147472/1‑122 (MQ=17)
                                                                aTAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCAGACa                           <  1:2147472/122‑1 (MQ=17)
                                                                aTAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCAGACATACAATATCGCAATGaaa         >  1:1376141/1‑140 (MQ=18)
                                                                aTAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCAGACATACAATATCGCAATGaaa         <  2:2424649/140‑1 (MQ=17)
                                                                aTAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGGCAAGATAGCCTGAATCAGACAAACAATATCGCAATGaaa         >  2:783830/1‑140 (MQ=11)
                                                                aTAATCACATCTCTCCTTATCACATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCAGACATACAATATCGCAATGaaa         <  1:1084414/140‑1 (MQ=17)
                                                                  aaTCACATCTCTCCTTATCATATTCCGCATTTTTAGTAATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCAGACATACAATATCGCAATGAAAAt       <  1:1192517/140‑1 (MQ=12)
                                                                  aaTCACATCTCTCCTTATCATATTCAGCATTCTTAATCAGGCTCTAATCAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCAGACATACAATATCGCAATGAAAAt       <  2:888385/140‑1 (MQ=11)
                                                                   aTCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGTGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCAGACATACAATATCGCAATGAAAATc      >  2:189935/1‑140 (MQ=17)
                                                                   aTCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCAGACATACAATATCGCAATGaaa         <  2:1502782/137‑1 (MQ=17)
                                                                   aTCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCAGACATACAATATCGCAATGaaa         >  1:1502782/1‑137 (MQ=18)
                                                                   aTCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCAGACATACAATATCGCAATGAAATTc      >  2:1487603/1‑140 (MQ=18)
                                                                   aTCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCAGACATACAATATCGCAATGAAAATc      >  2:2201774/1‑140 (MQ=18)
                                                                   aTCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCAGACATACAATATCGCAATGAAAATc      >  2:2367237/1‑140 (MQ=17)
                                                                   aTCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCAGACATACAATATCGCAATGAAAATc      >  2:1933485/1‑140 (MQ=18)
                                                                   aTCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCAGACATACAATATCGCAATGAAAATc      >  2:555975/1‑140 (MQ=18)
                                                                   aTCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCAGACATACAATATCGCAATGAAAATc      >  1:1602278/1‑140 (MQ=18)
                                                                      acaTCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCAGACATACAATATCGCAATg            >  1:1241117/1‑131 (MQ=17)
                                                                      acaTCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCAGACATACAATATCGCAATg            <  2:1241117/131‑1 (MQ=17)
                                                                       caTCTCTCCTTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCAGACATACAATATCGCAATGAAAATCaata  <  1:414569/140‑1 (MQ=14)
                                                                       caTCTCTCATTATCATATTCAGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCAGACATACAATATCGCAATGAAAATCaata  <  1:2423072/140‑1 (MQ=14)
                                                                                                            |                                                                                                      
AAAAATTGTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCCCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCAGACATACAATATCGCAATGAAAATCAATA  >  NC_000913/1632841‑1633051

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