Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,362,956 C→A 12.7% S50S (TCC→TCA)  yhcD → putative outer membrane fimbrial subunit usher protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,362,9560CA12.7% 263.3 / 3.5 103S50S (TCC→TCA) yhcDputative outer membrane fimbrial subunit usher protein
Reads supporting (aligned to +/- strand):  ref base C (50/36);  new base A (15/2);  total (65/38)
Fisher's exact test for biased strand distribution p-value = 2.59e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 110 total aligned reads displayed.

TGTTAAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTT  >  NC_000913/3362808‑3363104
                                                                                                                                                    |                                                                                                                                                    
tGTTAAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTcc                                                                                                                                                      <  1:799515/149‑1 (MQ=255)
tGTTAAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTcc                                                                                                                                                      <  2:74771/149‑1 (MQ=255)
 ttaaaaaaaaaCGTTACTGGCCTACACCATTGTTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTccc                                                                                                                                                     <  1:1073471/146‑1 (MQ=255)
 gtaaaaaaaaaCGTTACTGGCCTACACCATTAGTGTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTATGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTccc                                                                                                                                                     <  1:451144/146‑1 (MQ=255)
  ttAAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCg                                                                                                                                                    >  1:1297895/1‑149 (MQ=255)
       aaaaaCGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAg                                                                                                                                               <  2:801269/149‑1 (MQ=255)
         aaaCGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATATTTCACTTTCAGcc                                                                                                                                             >  2:1346296/1‑149 (MQ=255)
                    gCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGcc                                                                                                                                  >  1:846784/1‑149 (MQ=255)
                    gCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCATAATATAACGCATTATTTTTCACGTTCAGCCCGTTTTTTGcc                                                                                                                                  >  2:1330908/1‑149 (MQ=255)
                     ccTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGAAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCa                                                                                                                                 >  1:722643/1‑149 (MQ=255)
                      cTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAg                                                                                                                                <  2:298969/149‑1 (MQ=255)
                      cTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAg                                                                                                                                >  1:662207/1‑149 (MQ=255)
                             aTTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGTAAACATTAAAATCACTAGGTTTATATCCTAATATATCGCAGTATTTTTCCCGTTCAGCCCGTTTTTAGCCAGGTGAata                                                                                                                         >  1:292444/1‑149 (MQ=255)
                                 gTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCa                                                                                                                     <  2:609753/149‑1 (MQ=255)
                                  ttttGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCAc                                                                                                                    >  1:71845/1‑149 (MQ=255)
                                  ttttGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCAAAATATAACCCATAAAATTAAAAAAAAAACCCTTTTTTTGCCAGGTGAATATTCAc                                                                                                                    >  2:1394925/1‑149 (MQ=255)
                                   tttGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCTGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGACCGTTTTTTGCCAGGTGAATATTCACt                                                                                                                   >  1:480903/1‑149 (MQ=255)
                                    ttGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCCCATTATTTTTAAAATTCATCACGTTTTTTGCCAGGTGAATATTCACtt                                                                                                                  >  2:1083374/1‑148 (MQ=255)
                                    ttGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCAAAATATAACCCATTAATTTAAAAATAAAACCCGTTTTTTGCCATGTGAATATTCACTg                                                                                                                  >  2:1046969/1‑149 (MQ=255)
                                     tGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGCTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGa                                                                                                                 >  1:960082/1‑149 (MQ=255)
                                       ccTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATa                                                                                                               <  2:485199/149‑1 (MQ=255)
                                        cTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAg                                                                                                              >  1:1273174/1‑149 (MQ=255)
                                           ttCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTTACTGATACTAt                                                                                                           <  1:1499650/149‑1 (MQ=255)
                                           attttgcccAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTAt                                                                                                           <  2:1347140/143‑1 (MQ=255)
                                                ccccAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATc                                                                                                          >  2:1379227/1‑145 (MQ=255)
                                                ccccAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTa                                                                                                      >  1:990921/1‑149 (MQ=255)
                                                ccccAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCAAGGTGAATATTCACTGATAGTATCAGTa                                                                                                      >  1:1559801/1‑149 (MQ=255)
                                                 cccAGCCAATGCAGATCGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTACATCCTAAGATCTCGCATTATTTTTCCAGTTCAGCGACTTTTTTGCCAGGTGAATATACACTGAGAGTAACAGTaa                                                                                                     >  1:13609/1‑149 (MQ=255)
                                                  ccAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTaaa                                                                                                    <  1:826518/149‑1 (MQ=255)
                                                            cAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGaaaaa                                                                                          <  1:1459133/149‑1 (MQ=255)
                                                            cAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGaaaaa                                                                                          >  1:1195151/1‑149 (MQ=255)
                                                              gATTGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGaaaaaaa                                                                                        <  2:169723/149‑1 (MQ=255)
                                                              gATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGaaaaaaa                                                                                        >  1:1023893/1‑149 (MQ=255)
                                                                 ggTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAgg                                                                                     <  2:1297895/149‑1 (MQ=255)
                                                                  gTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGc                                                                                    >  1:1566194/1‑149 (MQ=255)
                                                                    ttaGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCaa                                                                                  <  2:733410/148‑1 (MQ=255)
                                                                     taGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCATAATATAACCCATAAAATTAAAAAAAAACCCAGTTTTTTGCCAGGTCAATATTCACTGATAGTATCAGTATATGGCGAAAAAAAAGGCAAc                                                                                 >  2:94832/1‑149 (MQ=255)
                                                                     taGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCAAAATATAACCCATAAAATTAAAAAAAAAAACAGTTTTTTGCAAGGTGAATTTTCACTTATAGTATCAGTATATGGAGAAATAAAAGGCaaa                                                                                 >  2:1540319/1‑148 (MQ=255)
                                                                     taGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATACTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAAc                                                                                 >  1:1010470/1‑149 (MQ=255)
                                                                      agagATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACa                                                                                >  1:1231546/1‑149 (MQ=255)
                                                                       gagaTTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCTTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAACTGCAACAt                                                                               >  1:297674/1‑149 (MQ=255)
                                                                          aTTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACAtttc                                                                            >  1:1164277/1‑147 (MQ=255)
                                                                             gCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCAAAATATAACCCATTAAATTAAAAAAAAAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTAc                                                                         >  2:1368268/1‑149 (MQ=255)
                                                                               cGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTAcgc                                                                       >  1:1461735/1‑149 (MQ=255)
                                                                                gCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGAAACATTGCTAcgcg                                                                      >  1:315503/1‑149 (MQ=255)
                                                                                   gTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCtt                                                                   <  2:71845/149‑1 (MQ=255)
                                                                                      gATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGCTCCTAATATATCGCATTATTTTTCCAGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATATTATCAGTAAATGTCGAAAAAAAATGCAAAATTGCTTCGCTCTTtga                                                                >  1:315054/1‑149 (MQ=255)
                                                                                         tttGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTtgatga                                                             >  1:741998/1‑149 (MQ=255)
                                                                                         tttGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTtgatga                                                             >  1:1013453/1‑149 (MQ=255)
                                                                                              tagggAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATg                                                        <  1:918581/147‑1 (MQ=255)
                                                                                               cGGGAAACATTAAAATCACTAGGTGTAGATCATAATATAACGCATTAATTTTCACGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATgg                                                       >  2:1466853/1‑149 (MQ=255)
                                                                                                tgtataCATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCACGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGc                                                      <  1:117136/144‑1 (MQ=255)
                                                                                                gggAAACATTAAAATCACTAGGTGTAGATCAAAATATAACCCATAAAATTAAAAAAAAATCCCGTTTTTTGCCAGGTGAATATTCTCTGATAGTATCAGTTAATGGCGAAAATAAAGGCTACATTGCTACGCGCTTTGATGATTATggg                                                      >  2:637789/1‑148 (MQ=255)
                                                                                                 ggAAACATTAAAATCACTAGGTGTAGATCAAAATATAACCCATAAAATTAAAAAAAAAACACGTTTTTTGCAACGTGAATATTCACTCATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGCAAATGGCg                                                     >  2:1159021/1‑149 (MQ=255)
                                                                                                  taaaCATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGa                                                    <  1:1346296/148‑1 (MQ=255)
                                                                                                   aaaCATTAAAATCACTAGGTGTAGATCAAAATATAACCCATAAAATTACAATAACAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGTGAc                                                   >  2:1584914/1‑149 (MQ=255)
                                                                                                   aaaCATTAAAATAACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGAc                                                   <  2:1116826/149‑1 (MQ=255)
                                                                                                       aTTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACAttt                                               >  1:373874/1‑149 (MQ=255)
                                                                                                        ttAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACAtttt                                              >  1:327108/1‑148 (MQ=255)
                                                                                                        ttAAAATCACTAGGTGTAGATCCTAATATATAGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCAACTCGCTTTGATGAAAATGGCGACAtttt                                              >  1:334942/1‑148 (MQ=255)
                                                                                                        ttAAAATCACTAGGTGTAGATCAAAATATAACCCATAAAATTAAAAAAAAAACCCGTTTTTTGCCAGGTGAATATTAACTGATAGTATCAGTAAATGGCGAAAAAATAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTg                                              >  2:759624/1‑149 (MQ=255)
                                                                                                          aaaaTCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTc                                            <  1:800373/149‑1 (MQ=255)
                                                                                                           aaaTCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCTCTTTGATGAAAATGGATACATTTGTCt                                           >  1:290440/1‑149 (MQ=255)
                                                                                                            aaTCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCtt                                          <  2:846784/149‑1 (MQ=255)
                                                                                                            aaTCACTAGGTGTAGATACTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCtt                                          <  2:480903/149‑1 (MQ=255)
                                                                                                               cACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGAt                                       <  2:577925/149‑1 (MQ=255)
                                                                                                                 cTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATAAGTAATTGGCGAAAAAAAAGGCAACATTGCTACGAGCTTTGATGAAAATGGAGACATTTGTCTTGATCa                                     >  1:1080996/1‑149 (MQ=255)
                                                                                                                  tAGGTGTAGATCCTAATATATCTCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAg                                    <  1:1322781/149‑1 (MQ=255)
                                                                                                                   aGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAgg                                   <  1:1330908/149‑1 (MQ=255)
                                                                                                                     gtgtAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCATGTGAATATTCACTTATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCa                                 <  1:1394925/149‑1 (MQ=255)
                                                                                                                      tgtAGCTCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCAt                                >  1:954139/1‑149 (MQ=255)
                                                                                                                      tgtAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCAt                                >  1:453573/1‑149 (MQ=255)
                                                                                                                       gtAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCAtt                               >  1:809194/1‑149 (MQ=255)
                                                                                                                        tAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCAttt                              >  1:353679/1‑149 (MQ=255)
                                                                                                                        tAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCAttt                              >  1:63876/1‑149 (MQ=255)
                                                                                                                        tAGATCATAATATATCGAATTATTTTTCCCGTTCATCAAGTTTTTTGCGAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCAttt                              <  1:921295/149‑1 (MQ=255)
                                                                                                                         aGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTc                             >  1:408801/1‑149 (MQ=255)
                                                                                                                          taTCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCt                            <  2:712060/148‑1 (MQ=255)
                                                                                                                          gATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCt                            >  2:449944/1‑149 (MQ=255)
                                                                                                                             ccTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGca                         >  1:1237692/1‑149 (MQ=255)
                                                                                                                             ccTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGca                         >  2:530954/1‑149 (MQ=255)
                                                                                                                             ccTAATATATAGCATTATTTTTCCCGTTCAGCACGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGAGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTACTGca                         >  1:1401631/1‑149 (MQ=255)
                                                                                                                              cTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGcaa                        <  2:480844/149‑1 (MQ=255)
                                                                                                                               tAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGcaac                       >  2:1535913/1‑149 (MQ=255)
                                                                                                                                 atatatCGCATTATTTTTCCCGTTCATCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGcaacaa                     <  1:785430/149‑1 (MQ=255)
                                                                                                                                  tatatCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTTAATGGAGAAAAAAAAGGCAACATTGCTACGCGATTTGATGAAAATGGATACATTTGTCTTGATCCGGCATTTCTCCAACAAg                    >  1:152959/1‑149 (MQ=255)
                                                                                                                                   tttttggatttttTTTTCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGc                   <  2:815900/142‑1 (MQ=255)
                                                                                                                                   ttattgtattattTTTCACGTTCAGGCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGc                   <  2:315503/142‑1 (MQ=255)
                                                                                                                                     atCGCATTATTTTTCACGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCg                 >  2:190471/1‑149 (MQ=255)
                                                                                                                                         cATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCTACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGttt             >  2:1104190/1‑149 (MQ=255)
                                                                                                                                         cATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGttt             >  1:1095002/1‑149 (MQ=255)
                                                                                                                                          attattTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTa            >  1:1220421/1‑149 (MQ=255)
                                                                                                                                          attattTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTa            >  2:895474/1‑149 (MQ=255)
                                                                                                                                              tttttCCCGTTCAGCCCGTTTTTTGCCATGTGAATATTCACTGATAGTATCAGTATATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTaaaaa        <  1:1083374/149‑1 (MQ=255)
                                                                                                                                              tttttCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAGATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTaaaaa        <  1:1046969/149‑1 (MQ=255)
                                                                                                                                              tttttCCCGTTAAGCACGTTTTTTGCCAGTTGAATATTCACTGATAGTATCAGTAAATGGAGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTaaaaa        <  2:1121245/149‑1 (MQ=255)
                                                                                                                                               ttttCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAAt       >  1:1521402/1‑149 (MQ=255)
                                                                                                                                                 ttCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTc     <  2:152599/149‑1 (MQ=255)
                                                                                                                                                    ccGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCtt  <  1:449944/149‑1 (MQ=255)
                                                                                                                                                    ccGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTACTGCAACAAGCCGGTTTAAAAATTCCtt  <  1:1466853/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TGTTAAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTGATCGGGAAACATTAAAATCACTAGGTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTT  >  NC_000913/3362808‑3363104

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