Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,274,629 A→T 5.1% P282P (CCA→CCT)  ychO → putative invasin YchO

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,274,6290AT 5.1% 260.2 / 2.8 79P282P (CCA→CCT) ychOputative invasin YchO
Reads supporting (aligned to +/- strand):  ref base A (41/34);  new base T (2/2);  total (43/36)
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.95e-01

TCGGATGCGCATGCCGTTCTATCAACACCTCAATACCAGTGTCAGCCTAGAACAGTATTTTGGTGATCGTGTTGATTTGTTTAACTCTGGTACGGGTTATCACAATCCCGTCGCGTTGAGTCTGGGATTAAATTACACCCCTGTGCCATTAGTCACTGTGACGGCCCAGCATAAACAGGGTGAAAGTGGCGAAAATCAAAATAACCTCGGGCTGAATCTTAATTACCGCTTTGGTGTACCGCTCAAAAAACAACTTTCTGCGGGCGAGGTTGCCGAAAGTCAGTCGTTACGTG  >  NC_000913/1274482‑1274774
                                                                                                                                                   |                                                                                                                                                 
tCGGATGCGCATGCCGTTCTATCAACACCTCAATACCAGTGTCAGCCTAGAACAGTATTTTGGTGATCGTGTTGATTTGTTTAACTCTGGTACGGGTTATCACAATCCCGTCGCGTTGAGTCTGGGATTAAATTACACCCCTGTGCCAtt                                                                                                                                                 <  1:1663117/150‑1 (MQ=255)
    aTGCGCATGCCGTTCTATCAACACCTCAATACCAGTGTCAGCCTAGAACAGTATTTTGGTGATCGTGTTGATTTGTTTAACTCTGGTACGGGTTATCACAATCCCGTCGCGTTGAGTCTGGGATTAAATTACACCCCTGTGCCATTAGTc                                                                                                                                             <  2:1258963/150‑1 (MQ=255)
     tGCGCATGCCGTTCTATCAACACCTCAATACCAGTGTCAGCCTAGAACAGTATTTTGGTGATCGTGTTGATTTGTTTAACTCTGGTACGGGTTATCACAATCCCGTCGCGTTGAGTCTGGGATTAAATTACACCCCTGTGCCATTAGTCa                                                                                                                                            >  2:1872590/1‑150 (MQ=255)
         cATGCCGTTCTATCAACACCTCAATACCAGTGTCAGCCTAGAACAGTATTTTGGTGATCGTGTTGATTTGTTTAACTCTGGTACGGGTTATCACAATCCCGTCGCGTTGAGTCTGGGATTAAATTACACCCCTGTGCCATTAGTCACtgt                                                                                                                                        <  1:1164561/150‑1 (MQ=255)
           tGCCGTTCTATCAACACCTCAATACCAGTGTCAGCCTAGAACAGTATTTTGGTGATCGTGTTGATTTGTTTAACTCTGGTACGGGTTATCACAATCCCGTCGCGTTGAGTCTGGGATTAAATTACACCCCTGTGCCATTAGTCACTGTGa                                                                                                                                      >  1:1468814/1‑150 (MQ=255)
             ccGTTCTATCAACACCTCAATACCAGTGTCAGCCTAGAACAGTATTTTGGTGATCGTGTTGATTTGTTTAACTCTGGTACGGGTTATCACAATCCCGTCGCGTTGAGTCTGGGATTAAATTACACCCCTGTGCCATTAGTCACTGTGACg                                                                                                                                    >  1:667965/1‑150 (MQ=255)
              cGTTCTATCAACACCTCAATACCAGTGTCAGCCTAGAACAGTATTTTGGTGATCGTGTTGATTTGTTTAACTCTGGTACGGGTTATCACAATCCCGTCGCGTTGAGTCTGGGATTAAATTACACCCCTGTGCCATTAGTc                                                                                                                                             <  2:1320644/140‑1 (MQ=255)
              cGTTCTATCAACACCTCAATACCAGTGTCAGCCTAGAACAGTATTTTGGTGATCGTGTTGATTTGTTTAACTCTGGTACGGGTTATCACAATCCCGTCGCGTTGAGTCTGGGATTAAATTACACCCCTGTGCCATTAGTc                                                                                                                                             >  1:1320644/1‑140 (MQ=255)
                  cTATCAACACCTCAATACCAGTGTCAGCCTAGAACAGTATTTTGGTGATCGTGTTGATTTGTTTAACTCTGGTACGGGTTATCACAATCCCGTCGCGTTGAGTCTGGGATTAAATTACACCCCTGTGCCATTAGTCACTGTGACGGCCCa                                                                                                                               >  2:1987771/1‑150 (MQ=255)
                   tATCAACACCTCAATACCAGTGTCAGCCTAGAACAGTATTTTGGTGATCGTGTTGATTTGTTTAACTCTGGTACGGGTTATCACAATCCCGTCGCGTTGAGTCTGGGATTAAATTACACCCCTGTGCCTTTAGTCACTGTGACGGCCCAg                                                                                                                              >  2:1552694/1‑150 (MQ=255)
                   tATCAACACCTCAATACCAGTGTCAGCCTAGAACAGTATTTTGGTGATCGTGTTGATTTGTTTAACTCTGGTACGGGTTATCACAATCCCGTCGCGTTGAGTCTGGGATTAAATTACACCCCTGTGCCTTTAGTCACTGTGACGGCCCAg                                                                                                                              >  2:1560095/1‑150 (MQ=255)
                    aTCAACACCTCAATACCAGTGTCAGCCTAGAACAGTATTTTGGTGATCGTGTTGATTTGTTTAACTCTGGTACGGGTTATCACAATCCCGTCGCGTTGAGTCTGGGATTAAATTACACCCCTGTGCCATTAGTCACTGTGACGGCCCAGc                                                                                                                             <  1:934444/150‑1 (MQ=255)
                           ccTCAATACCAGTGTCAGCCTAGAACAGTATTTTGGTGATCGTGTTGATTTGTTTAACTCTGGTACGGGTTATCACAATCCCGTCGCGTTGAGTCTGGGATTAAATTACACCCCTGTGCCATTAGTCACTGTGACGGCCCAGCATAAACa                                                                                                                      >  1:458257/1‑150 (MQ=255)
                             tCAATACCAGTGTCAGCCTAGAACAGTATTTTGGTGATCGTGTTGATTTGTTTAACTCTGGTACGGGTTATCACAATCCCGTCGCGTTGAGTCTGGGATTAAATTACACCCCTGTGCCATTAGTCACTGTGACGGCCCAGCATAAACAgg                                                                                                                    >  1:1641066/1‑150 (MQ=255)
                             tCAATACCAGTGTCAGCCTAGAACAGTATTTTGGTGATCGTGTTGATTTGTTTAACTCTGGTACGGGTTATCACAATCCCGTCGCGTTGAGTCTGGGATTAAATTACACCCCTGTGCCATTAGTCACTGTGACGGCCCAGCATAAACAgg                                                                                                                    >  2:1500449/1‑150 (MQ=255)
                             tCAATACCAGTGTCAGCCTAGAACAGTATTTTGGTGATCGTGTTGATTTGTTTAACTCTGGTACGGGTTATCACAATCCCGTCGCGTTGAGTCTGGGATTAAATTACACCCCTGTGCCATTAGTCACTGTGACGGCCCAGCATAAACAgg                                                                                                                    >  1:1412490/1‑150 (MQ=255)
                                 tACCAGTGTCAGCCTAGAACAGTATTTTGGTGATCGTGTTGATTTGTTTAACTCTGGTACGGGTTATCACAATCCCGTCGCGTTGAGTCTGGGATTAAATTACACCCCTGTGCCATTAGTCACTGTGACGGCCCAGCATAAACAGGGTGa                                                                                                                >  1:1752122/1‑150 (MQ=255)
                                     aGTGTCAGCCTAGAACAGTATTTTGGTGATCGTGTTGATTTGTTTAACTCTGGTACGGGTTATCACAATCCGGTCGCGTTGAGTCTGGGATTAAATTAAACCCCTGTGCCATTAGTCCATTTGACGGCCAAGCATAAACAGGTTGAAAGt                                                                                                            >  1:1529737/1‑150 (MQ=255)
                                      gtgtCAGCCTAGAACAGTATTTTGGTGATCGTGTTGATTTGTTTAACTCTGGTACGGGTTATCACAATCCCGTCGCGTTGAGTCTGGGATTAAATTACACCCCTGTGCCATTAGTCACTGTGACGGCCCAGCATAAACAGGGTGAAAGTg                                                                                                           >  2:925807/1‑150 (MQ=255)
                                        gtCAGCCTAGAACAGTATTTTGGTGATCGTGTTGATTTGTTTAACTCTGGTACGGGTTATCACAATCCCGTCGCGTTGAGTCTGGGATTAAATTACACCCCTGTGCCATTAGTCACTGTGACGGCCCAGCATAAACAGGGTGAAAGTGGc                                                                                                         <  1:1872590/150‑1 (MQ=255)
                                        gtCAGCCTAGAACAGTATTTTGGTGATCGTGTTGATTTGTTTAACTCTGGTACGGGTTATCACAATCCCATCGCGTTGAGTCTGGGATTAAATTACACCCCTGTGCCATTAGTCACTGTGACGGCCCAGCATAAACAGGGTGAAAGTGGc                                                                                                         <  2:1047115/150‑1 (MQ=255)
                                        gtCAGCCTAGAACAGTATTTTGGTGATCGTGTTGATTAGTTTAACTCTGGTACGGGTTATCACAATCCCGTCGCGTTGAGTCTGGGATTAAATTACACCCCTGTGCCATTAGTCACTGTGACGGCCCAGCATAAACAGGGTGAAAGTGGc                                                                                                         >  1:1012005/1‑150 (MQ=255)
                                        ggcAGCCTAGCCCAGTATTTTGGTGATCGTGTTGATTTGTTTAACTCTGGTACGGGTTATCACAATCCCGTCGCGTTGAGGCTGGGATTAAATTACACCCCTGTGCCATTAGTCACTGTGACGGCCCAGCATAAACAGGGTGAAAGTGGc                                                                                                         <  2:311647/148‑1 (MQ=255)
                                            gCCTAGAACAGTATTTTGGTGATCGTGTTGATTTGTTTAACTCTGGTACGGGTTATCACAATCCCGTCGCGTTGAGTCTGGGATTAAATTACACCCCTGTGCCATTAGTCACTGTGACGGCCCAGCATAAACAGGGTGAAAGTGGCGaaa                                                                                                     >  2:1857795/1‑150 (MQ=255)
                                               tAGAACAGTATTTTGGTGATCGTGTTGATTTGTTTAACTCTGGTACGGGTTATCACAATCCCGTCGCGTTGAGTCTGGGATTAAATTACACCCCTGTGCCATTAGTCACTGTGACGGCCCAGCATAAACAGGGTGAAAGTGGCGAAAATc                                                                                                  >  1:626993/1‑150 (MQ=255)
                                                 gAACAGTATTTTGGTGATCGTGTTGATTTGTTTAACTCTGGTACGGGTTATCACAATCCCGTCGCGTTGAGTCTGGGATTAAATTACACCCCTGTGCCATTAGTCACTGTGACGGCCCAGCATAAACAGGGTGAAAGTGGCGAAAATCaa                                                                                                >  1:1657158/1‑150 (MQ=255)
                                                 gAACAGTATTTTGGTGATCGTGTTGATTTGTTTAACTCTGGTACGGGTTATCACAATCCCGTCGCGTTGAGTCTGGGATTAAATTACACCCCTGTGCCATTAGTCACTGTGACGGCCCAGCATAAACAGGGTGAAAGTGGCGAAAATCaa                                                                                                >  2:1428656/1‑150 (MQ=255)
                                                 gAACAGTATTTTGGTGATCGTGTTGATTTGTTTAACTCTGGTACGGGTTATCACAATCCCGTCGCGTTGAGTCTGGGATTAAATTACACCCCTGTGCCATTAGTCACTGTGACGGCCCAGCATAAACAGGGTGAAAGTGGCGAAAATCaa                                                                                                >  1:909904/1‑150 (MQ=255)
                                                   aCAGTATTTTGGTGATCGTGTTGATTTGTTTAACTCTGGTACGGGTTATCACAATCCCGTCGCGTTGAGTCTGGGATTAAATTACACCCCTGTGCCATTAGTCACTGTGACGGCCCAGCATAAACAGGGTGAAAGTGGCGAAAATCaaaa                                                                                              <  1:98094/150‑1 (MQ=255)
                                                              gTGATCGTGTTGATTTGTTTAACTCTGGTACGGGTTATCACAATCCCGTCGCGTTGAGTCTGGGATTAAATTACACCCCTGTGCCATTAGTCACTGTGACGGCCCAGCATAAACAGGGTGAAAGTGGCTAAAATCAAAATAACCTCGGGc                                                                                   <  1:1697615/150‑1 (MQ=255)
                                                               tGATCGTGTTGATTTGTTTAACTCTGGTACGGGTTATCACAATCCCGTCGCGTTGAGTCTGGGATTAAATTACACCCCTGTGCCATTAGTCACTGTGACGGCCCAGCATAAACAGGGTGAAAGTGGCGAAAATCAAAATAACCTCGGGCt                                                                                  >  1:629972/1‑150 (MQ=255)
                                                               tGATCGTGTTGATTTGTTTAACTCTGGTACGGGTTATCACAATCCCGTCGCGTTGAGTCTGGGATTAAATTACACCCCTGTGCCATTAGTCACTGTGACGGCCCAGCATAAACAGGGTGAAAGTGGCGAAAATCAAAATAACCTCGGGCt                                                                                  >  1:1153834/1‑150 (MQ=255)
                                                               tGATCGTGTTGATTTGTTTAACTCTGGTACGGGTTATCACAATCCCGTCGCGTTGAGTCTGGGATTAAATTAAACCCCTGTGCCATTAGTCACTGTGACGGCCCAGCATAAACAGGGTGAAAGTGGCGAAAATCAAAATAACCTCGGGCt                                                                                  >  1:403371/1‑150 (MQ=255)
                                                                       ttGATTTGTTTAACTCTGGTACGGGTTATCACAATCCCGTCGCGTTGAGTCTGGGATTAAATTACACCCCTGTGCCATTAGTCACTGTGACGGCCCAGCATAAACAGGGTGAAAGTGGCGAAAATCAAAATAACCTCGGGCTGAATCTTa                                                                          <  1:1444644/150‑1 (MQ=255)
                                                                        tGATTTGTTTAACTCTGGTACGGGTTATCACAATCCCGTCGCGTTGAGTCTGGGATTAAATTACACCCCTGTGCCATTAGTCACTGTGACGGCCCAGCATAAACAGGGTGAAAGTGGCGAAAATCAAAATAACCTCGGGCTGAATCTTaa                                                                         <  2:1472360/150‑1 (MQ=255)
                                                                        tGATTTGTTTAACTCTGGTACGGGTTATCACAATCCCGTCGCGTTGAGTCTGGGATTAAATTACACCCCTGTGCCATTAGTCACTGTGACGGCCCAGCATAAACAGGGTGAAAGTGGCGAAAATCAAAATAACCTCGGGCTGAATCTTaa                                                                         <  2:1729588/150‑1 (MQ=255)
                                                                              gTTTAACTCTGGTACGGGTTATCATAATCCCGTCGCGTTGAGTCTGGGATTAAATTACACCCCTGTGCCATTAGTCACTGTGACGGCCCAGCATAAACAGGGTGAAAGTGGCGAAAATCAAAATAACCTCGGGCTGAATCTTAATTACCg                                                                   >  1:1739881/1‑150 (MQ=255)
                                                                              gTTTAACTCTGGTACGGGTTATCACAATCCCGTCGCGTTGAGTCTGGGATTAAATTACACCCCTGTGCCATTAGTCACTGTGACGGCCCAGCATAAACAGGGTGAAAGTGGCGAAAATCAAAATAACCTCGGGCTGAATCTTAATTACCg                                                                   <  1:982500/150‑1 (MQ=255)
                                                                               tttAACTCTGGTACGGGTTATCACAATCCCGTCGCGTTGAGTCTGGGATTAAATTACACCCCTGTGCCATTAGTCACTGTGACGGCCCAGCATAAACAGGGTGAAAGTGGCGAAAATCAAAATAACCTCGGGCTGAATCTTAATTACCGc                                                                  <  2:1493290/150‑1 (MQ=255)
                                                                                 tAACTCTGGTACGGGTTATCACAATCCCGTCGCGTTGAGTCTGGGATTAAATTACACCCCTGTGCCATTAGTCACTGTGACGGCCCAGCATAAACAGGGTGAAAGTGGCGAAAATCAAAATAACCTCGGGCTGAATCTTAATTACCGCtt                                                                <  2:1153834/150‑1 (MQ=255)
                                                                                 tAACTCTGGTACGGGTTATCACAATCCCGTCGCGTTGAGTCTGGGATTAAATTACACCCCTGTGCCATTAGTCACTGTGACGGCCCAGCATAAACAGGGTGAAAGTGGCGAAAATCAAAATAACCTCGGGCTGAATCTTAATTACCGCtt                                                                >  2:577245/1‑150 (MQ=255)
                                                                                 tAACTCTGGTACGGGTTATCACAATCCCGTCGCGTTGAGTCTGGGATTAAATTACACCCCTGTGCCATTAGTCACTGTGACGGCCCAGCATAAACAGGGTGAAAGTGGCGAAAATCAAAATAACCTCGGGCTGAATCTTAATTACCGCtt                                                                >  2:587314/1‑150 (MQ=255)
                                                                                 tAACTCTGGTACGGGTTATCACAATCCCGTCGCGTTGAGTCTGGGATTAAATTACACCCCTGTGCCATTAGTCACTGTGACGGCCCAGCATAAACAGGGTGAAAGTGGCGAAAATCAAAATAACCTCGGGCTGAATCTTAATTACCGCtt                                                                >  2:952314/1‑150 (MQ=255)
                                                                                   aCTCTGGTACGGGTTATCACAATCCCGTCGCGTTGAGTCTGGGATTAAATTACACCCCTGTGCCATTAGTCACTGTGACGGCCCAGCATAAACAGGGTGAAAGTGGCGAAAATCAAAATAACCTCGGGCTGAATCTTAATTACCGCTTTg                                                              <  2:629972/150‑1 (MQ=255)
                                                                                    ctctGGTACGGGTTATCACAATCCCGTCGCGTTGAGTCTGGGATTAAATTACACCCCTGTGCCATTAGTCACTGTGACGGCCCAGCATAAACAGGGTGAAAGTGGCGAAAATCAAAATAACCTCGGGCTGAATCTTAATTACCGCTTTgg                                                             <  1:1987771/150‑1 (MQ=255)
                                                                                     tctGGTACGGGTTATCACAATCCCGTCGCGTTGAGTCTGGGATTAAATTACACCCCTGTGCCATTAGTCACTGTGACGGCCCAGCATAAACAGGGTGAAAGTGGCGAAAATCAAAATAACCTCGGGCTGAATCTTAATTACCGCTTTGgt                                                            <  1:1277406/150‑1 (MQ=255)
                                                                                      ctGGTACGGGTTATCACAATCCCGTCGCGTTGAGTCTGGGATTAAATTACACCCCTGTGCCATTAGTCACTGTGACGGCCCAGCATAAACAGGGTGAAAGTGGCGAAAATCAAAATAACCTCGGGCTGAATCTTAATTACCGCTTTGgtg                                                           <  2:909904/150‑1 (MQ=255)
                                                                                        ggTACGGGTTATCACAATCCCGTCGCGTTGAGTCTGGGATTAAATTACACCCCTGTGCCATTAGTCACTGTGACGGCCCAGCATAAACAGGGTGAAAGTGGCGAAAATCAAAATAACCTCGGGCTGAATCTTAATTAACGCTTTGGTGTa                                                         >  2:499194/1‑150 (MQ=255)
                                                                                         gTACGGGTTATCACAATCCCGTCGCGTTGAGTCTGGGATTAAATTACACCCCTGTGCCATTAGTCACTGTGACGGCCCAGCATAAACAGGGTGAAAGTGGCGAAAATCAAAATAACCTCGGGCTGAATCTTAATTACCGCTTTGGTGTAc                                                        >  1:1909094/1‑150 (MQ=255)
                                                                                          tACGGGTTATCACAATCCCGTCGCGTTGAGTCTGGGATTAAATTACACCCCTGTGCCATTAGTCACTGTGACGGCCCAGCATAAACAGGGTGAAAGTGGCGAAAATCAAAATAACCTCGGGCTGAATCTTAATTACCGCTTTGGTGTAcc                                                       <  1:147460/150‑1 (MQ=255)
                                                                                           aCGGGTTATCACAATCCCGTCGCGTTGAGTCTGGGATTAAATTACACCCCTGTGCCATTAGTCACTGTGACGGCCCAGCATAAACAGGGTGAAAGTGGCGAAAATCAAAATAACCTCGGGCTGAATCTTAATTACCGCTTTGGTGTACCg                                                      >  2:708166/1‑150 (MQ=255)
                                                                                            cGGGTTATCACAATCCCGTCGCGTTGAGTCTGGGATTAAATTACACCCCTGTGCCATTAGTCACTGTGACGGCCCAGCATAAACAGGGTGAAAGTGGCGAAAATCAAAATAACCTCGGGCTGAATCTTAATTACCGCTTTGGTGTACCGc                                                     >  1:232398/1‑150 (MQ=255)
                                                                                              ggTTATCACAATCCCGTCGCGTTGAGTCTGGGATTAAATTACACCCCTGTGCCTTTAGTCACTGTGACGGCCCAGCATAAACAGGGTGAAAGTGGCGAAAATCAAAATAACCTCGGGCTGAATCTTAATTACCGCTTTGGTGTACCGCTc                                                   <  1:1552694/150‑1 (MQ=255)
                                                                                              ggTTATCACAATCCCGTCGCGTTGAGTCTGGGATTAAATTACACCCCTGTGCCTTTAGTCACTGTGACGGCCCAGCATAAACAGGGTGAAAGTGGCGAAAATCAAAATAACCTCGGGCTGAATCTTAATTACCGCTTTGGTGTACCGCTc                                                   <  1:1560095/150‑1 (MQ=255)
                                                                                              ggTTATCACAATCCCGTCGCGTTGAGTCTGGGATTAAATTACACCCCTGTGCCATTAGTCACTGTGACGGCCCAGCATAAACAGGGTGAAAGTGGCGAAAATCAAAATAACCTCGGGCTGAATCTTAATTACCGCTTTGGTGTACCGCTc                                                   >  1:375236/1‑150 (MQ=255)
                                                                                               gTTATCACAATCCCGTCGCGTTGAGTCTGGGATTAAATTACACCCCTGTGCCATTAGTCACTGTGACGGCCCAGCATAAACAGGGTGAAAGTGGCGAAAATCAAAATAACCTCGGGCTGAATCTTAATTACCGCTTTGGTGTACCGCTCa                                                  <  1:1428656/150‑1 (MQ=255)
                                                                                                 tATCACAATCCCGTCGCGTTGAGTCTGGGATTAAATTACACCCCTGTGCCATTAGTCACTGTGACGGCCCAGCATAAACAGGGTGAAAGTGGCGAAAATCAAAATAACCTCGGGCTGAATCTTAATTACCGCTTTGGTGTACCGCTCaaa                                                <  1:587314/150‑1 (MQ=255)
                                                                                                    cacaATCCCGTCGCGTTGAGTCTGGGATTAAATTACACCCCTGTGCCATTAGTCACTGTGACGGCCCAGCATAAACAGGGTGAAAGTGGCGAAAATCAAAATAACCTCGGGCTGAATCTTAATTACCGCTTTGGTGTACCGCTCaaaaaa                                             >  2:1308916/1‑150 (MQ=255)
                                                                                                        aTCCCGTCGCGTTGAGTCTGGGATTAAATTACACCCCTGTGCCATTAGTCACTGTGACGGCCCAGCATAAACAGGGTGAAAGTGGCGAAAATCAAAATAACCTCGGGCTGAATCTTAATTACCGCTTTGGTGTACCGCTCAAAAaacaac                                         <  1:764066/150‑1 (MQ=255)
                                                                                                         tCCCGTCGCGTTGAGTCTGGGATTAAATTACACCCCTGTGCCATTAGTCACTGTGACGGCCCAGCATAAACAGGGTGAAAGTGGCGAAAATCAAAATAACCTCGGGCTGAATCTTAATTACCGCTTTGGTGTACCGCTCAAAAAACAACt                                        <  2:458257/150‑1 (MQ=255)
                                                                                                         tCCCGTCGCGTTGAGTCTGGGATTAAATTACACCCCTGTGCCATTAGTCACTGTGACGGCCCAGCATAAACAGGGTGAAAGTGGCGAAAATCAAAATAACCTCGGGCTGAATCTTAATTACCGCTTTGGTGTACCGCTCAAAAAACAACt                                        >  1:79437/1‑150 (MQ=255)
                                                                                                                    tGAGTCTGGGATTAAATTACACCCCTGTGCCATTAGTCACTGTGACGGCCCAGCATAAACAGGGTGAAAGTGGCGAAAATCAAAATAACCTCGGGCTGAATCTTAATTACCGCTTTGGTGTACCGCTCAAAAAACAACTTTCTGCGGGCg                             <  2:1412490/150‑1 (MQ=255)
                                                                                                                       gTCTGGGATTAAATTACACCCCTGTGCCATTAGTCACTGTGACGGCCCAGCATAAACAGGGTGAAAGTGGCGAAAATCAAAATAACCTCGGGCTGAATCTTAATTACCGCTTTGGTGTACCGCTCAAAAAACAACTTTCTGCGGGCGAgg                          <  1:1857795/150‑1 (MQ=255)
                                                                                                                       gTCTGGGATTAAATTACACCCCTGTGCCATTAGTCACTGTGACGGCCCAGCATAAACAGGGTGAAAGTGGCGAAAATCAAAATAACCTCGGGCTGAATCTTAATTACCGCTTTGGTGTACCGCTCAAAAAACAACTTTCTGCGGGCGAgg                          >  1:1517818/1‑150 (MQ=255)
                                                                                                                       gTCTGGGATTAAATTACACCCCTGTGCCATTAGTCACTGTGACGGCCCAGCATAAACAGGGTGAAAGTGGCGAAAATCAAAATAACCTCGGGCTGAATCTTAATTACCGCTTTGGTGTACCGCTCAAAAAACAACTTTCTGCGGGCGAgg                          >  1:1481958/1‑150 (MQ=255)
                                                                                                                         cTGGGATTAAATTACACCCCTGTGCCATTAGTCACTGTGACGGCCCAGCATAAACAGGGTGAAAGTGGCGAAAATCAAAATAACCTCGGGCTGAATCTTAATTACCGCTTTGGTGTACCGCTCAAAAAACAACTTTCTGCGGGCGAGGtt                        >  1:1931705/1‑150 (MQ=255)
                                                                                                                          tGGGATTAAATTACACCCCTGTGCCATTAGTCACTGTGACGGCCCAGCATAAACAGGGTGAAAGTGGCGAAAATCAAAATAACCTCGGGCTGAATCTTAATTACCGCTTTGGTGTACCGCTCAAAAAACAACTTTCTGCGGGCGAGGTTg                       >  2:590132/1‑150 (MQ=255)
                                                                                                                            ggATTAAATTACACCCCTGTGCCATTAGTCACTGTGACGGCCCAGCATAAACAGGGTGAAAGTGGCGAAAATCAAAATAACCTCGGGCTGAATCTTAATTACCGCTTTGGTGTACCGCTCAAAAAACAACTTTCTGCGGGCGAGGTTGcc                     <  2:708822/150‑1 (MQ=255)
                                                                                                                            ggATTAAATTACACCCCTGTGCCATTAGTCACTGTGACGGCCCAGCATAAACAGGGTGAAAGTGGCGAAAATCAAAATAACCTCGGGCTGAATCTTAATTACCGCTTTGGTGTACCGCTCAAAAAACAACTTTCTGCGGGCGAGGTTGcc                     >  1:1248354/1‑150 (MQ=255)
                                                                                                                             gATTAAATTACACCCCTGTGCCATTAGTCACTGTGACGGCCCAGCATAAACAGGGTGAAAGTGGCGAAAATCAAAATAACCTCGGGCTGAATCTTAATTACCGCTTTGGTGTACCGCTCAAAAAACAACTTTCTGCGGGCGAGGTTGCCg                    >  2:1962995/1‑150 (MQ=255)
                                                                                                                               ttAAATTACACCCCTGTGCCATTAGTCACTGTGACGGCCCAGCATAAACAGGGTGAAAGTGGCGAAAATCAAAATAACCTCGGGCTGAATCTTAATTACCGCTTTGGTGTACCGCTCAAAAAACAACTTTCTGCGGGCGAGGTTGCCGaa                  <  2:626993/150‑1 (MQ=255)
                                                                                                                                 aaaTTACACCCCTGTGCCATTAGTCACTGTGACGGCCCAGCATAAACAGGGTGAAAGTGGCGAAAATCAAAATAACCTCGGGCTGAATCTTAATTACCGCTTTGGTGTACCGCTCAAAAAACAACTTTCTGCGGGCGAGGTTGCCGAAag                <  2:1909094/150‑1 (MQ=255)
                                                                                                                                  aaTTACACCCCTGTGCCATTAGTCACTGTGACGGCCCAGCATAAACAGGGTGAAAGTGGCGAAAATCAAAATAACCTCGGGCTGAATCTTAATTACCGCTTTGGTGTACCGCTCAAAAAACAACTTTCTGCGGGCGAGGTTGCCGAAagt               <  2:79437/150‑1 (MQ=255)
                                                                                                                                       cacCCCTGTGCCATTAGTCACTGTGACGGCCCAGCATAAACAGGGTGAAAGTGGCGAAAATCAAAATAACCTCGGGCTGAATCTTAATTACCGCTTTGGTGTACCGCTCAAAAAACAACTTTCTGCGGGCGAGGTTGCCGAAagtcagtc          >  1:1657102/1‑150 (MQ=255)
                                                                                                                                        acCCCTGTGCCATTAGTCACTGTGACGGCCCAGCATAAACAGGGTGAAAGTGGCGAAAATCAAAATAACCTCGGGCTGAATCTTAATTAAGTTTTTGGTTTACCGCTCAAAAAAAAAGTTTCTGAGTGAGAGGTTGCCGAAAGTCAGTCg         <  2:1468814/150‑1 (MQ=255)
                                                                                                                                           ccTGTGCCATTAGTCACTGTGACGGCCCAGCATAAACAGGGTGAAAGTGGCGAAAATCAAAATAACCTCGGGCTTAATCTTAATTACCGCTTTGGTGTACCGCTCAAAAAACAACTTTCTGCGGGCGAGGTTGCCGAAAGTCAGTCGTTa      <  1:704965/150‑1 (MQ=255)
                                                                                                                                           ccTGTGCCATTAGTCACTGTGACGGCCCAGCATAAACAGGGTGAAAGTGGCGAAAATCAAAATAACCTCGGGCTGAATCTTAATTACCGCTTTGGTGTACCGCTCAAAAAACAACTTTCTGCGGGCGAGGTTGCCGAAAGTCAGTCGTTa      <  1:499194/150‑1 (MQ=255)
                                                                                                                                            cTGTGCCATTAGTCACTGTGACGGCCCAGCATAAACAGGGTGAAAGTGGCGAAAATCAAAATAACCTCGGGCTGAATCTTAATTACCGCTTTGGTGTACCGCTCAAAAAACAACTTTCTGCGGGCGAGGTTGCCGAAAGTCAGTCGTTAc     >  2:15744/1‑150 (MQ=255)
                                                                                                                                               tgCCATTAGTCACTGTGACGGCCCAGCATAAACAGGGTGAAAGTGGCGAAAATCAAAATAACCTCGGGCTGAATCTTAATTACCGCTTTGGTGTACCGCTCAAAAAACAACTTTCTGCGGGCGAGGTTGCCGAAAGTCAGTCGTTACGTg  >  2:635185/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TCGGATGCGCATGCCGTTCTATCAACACCTCAATACCAGTGTCAGCCTAGAACAGTATTTTGGTGATCGTGTTGATTTGTTTAACTCTGGTACGGGTTATCACAATCCCGTCGCGTTGAGTCTGGGATTAAATTACACCCCTGTGCCATTAGTCACTGTGACGGCCCAGCATAAACAGGGTGAAAGTGGCGAAAATCAAAATAACCTCGGGCTGAATCTTAATTACCGCTTTGGTGTACCGCTCAAAAAACAACTTTCTGCGGGCGAGGTTGCCGAAAGTCAGTCGTTACGTG  >  NC_000913/1274482‑1274774

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