Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,588,479 T→A 6.8% N121Y (AAT→TAT)  ydeS ← putative fimbrial protein YdeS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,588,4790TA 6.8% 311.2 / 8.0 100N121Y (AAT→TAT) ydeSputative fimbrial protein YdeS
Reads supporting (aligned to +/- strand):  ref base T (49/44);  new base A (4/3);  total (53/47)
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 101 total aligned reads displayed.

GTTGCCGAAGCCCTCACCAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCT  >  NC_000913/1588331‑1588625
                                                                                                                                                    |                                                                                                                                                  
gTTGCCGAAGCCCTCACCAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATAttt                                                                                                                                                   >  2:840602/1‑150 (MQ=255)
 ttGCCGAAGCCCTCGCCAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTg                                                                                                                                                  <  2:243331/150‑1 (MQ=255)
 ttGCCGAAGCCCTCACCAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCAAATATg                                                                                                                                                  <  1:81484/150‑1 (MQ=255)
  tGCCGAAGCCCTCACCAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGc                                                                                                                                                 <  1:357546/150‑1 (MQ=255)
   gCCGAAGCCCTCACCAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTCCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCg                                                                                                                                                >  2:871251/1‑150 (MQ=255)
      gAAGCCCTCACCAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTcc                                                                                                                                             <  2:511872/150‑1 (MQ=255)
        aGCCCTCACCAGTCGCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCAAATATGCGTccgg                                                                                                                                           <  1:273419/150‑3 (MQ=255)
        aGCCCTCACCAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCaa                                                                                                                                           <  2:512963/150‑1 (MQ=255)
        aGCCCTCACCAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCAAATATGCGTccgg                                                                                                                                           <  2:1514671/150‑3 (MQ=255)
        aGCCCTCACCAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCAAATATGCGTccgg                                                                                                                                           >  1:728118/1‑148 (MQ=255)
          cccTCACCAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAata                                                                                                                                         >  2:579629/1‑150 (MQ=255)
           ccTCACCAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAatat                                                                                                                                        <  1:426937/150‑1 (MQ=255)
              cACCAGTCCCGGATTGACGGACTTCTTAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTc                                                                                                                                     <  1:1475657/150‑1 (MQ=255)
              cACCAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTc                                                                                                                                     >  1:1330476/1‑150 (MQ=255)
              cACCAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTc                                                                                                                                     >  2:1075427/1‑150 (MQ=255)
                ccAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCta                                                                                                                                   >  2:190846/1‑150 (MQ=255)
                   gTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATAc                                                                                                                                >  1:674814/1‑150 (MQ=255)
                     cccGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATAccc                                                                                                                              >  1:1464368/1‑150 (MQ=255)
                      ccGGATTGACGGACTTCTGAGGTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGAGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGAGGCATATTTGCGTCCAATTTTTCTATACCCa                                                                                                                             <  2:729664/150‑1 (MQ=255)
                           ttGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCg                                                                                                                        >  2:972525/1‑150 (MQ=255)
                            tGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTGAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCgg                                                                                                                       <  1:1047815/150‑1 (MQ=255)
                             gACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGa                                                                                                                      <  1:1075427/150‑1 (MQ=255)
                                   cTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGt                                                                                                                >  1:895317/1‑150 (MQ=255)
                                   cTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGt                                                                                                                <  2:1394926/150‑1 (MQ=255)
                                       tGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTg                                                                                                            <  1:552714/150‑1 (MQ=255)
                                          gTTGACTTCAGACGAGCGGAGTAAGTAAAAATATTGATCTGTTCTGGTAACAGTGGGATACAATGCATCCCGGCATGAAGATCATTAAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCTTATTGctt                                                                                                         >  2:1366628/1‑150 (MQ=255)
                                            tGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGcttcc                                                                                                       >  1:1399331/1‑150 (MQ=255)
                                              aCTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGcttcctt                                                                                                     <  2:1487487/150‑1 (MQ=255)
                                              aCTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGGCCAATATTTCTATACCCAAACCGGAGGCCGTATTGcttcctt                                                                                                     <  2:1356310/150‑1 (MQ=255)
                                               cTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTAACCGGACGAATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCCTATTGcttccttc                                                                                                    >  1:614404/1‑150 (MQ=255)
                                                 tCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCAt                                                                                                  <  2:923604/150‑1 (MQ=255)
                                                   aGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCCGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATcc                                                                                                >  1:1508546/1‑150 (MQ=255)
                                                    gACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCa                                                                                               >  1:1422560/1‑150 (MQ=255)
                                                    gACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCa                                                                                               >  2:321286/1‑150 (MQ=255)
                                                        aGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGttt                                                                                           >  2:67879/1‑150 (MQ=255)
                                                        aGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGttt                                                                                           <  1:1225196/150‑1 (MQ=255)
                                                         gCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTc                                                                                          <  1:1255215/150‑1 (MQ=255)
                                                          gggAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCAGTGCATCCCGGCATGATGATCATTCAGTTTGACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTca                                                                                         <  2:507452/149‑1 (MQ=255)
                                                          cGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTATca                                                                                         <  2:1422560/150‑1 (MQ=255)
                                                           ggAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTcaa                                                                                        >  2:1341360/1‑150 (MQ=255)
                                                            gAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTcaac                                                                                       <  1:1120779/150‑1 (MQ=255)
                                                                 aGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGcttcc                                                                                                       <  2:781722/129‑1 (MQ=255)
                                                                 aGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGcttcc                                                                                                       >  1:781722/1‑129 (MQ=255)
                                                                  ggCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGta                                                                                 >  1:857127/1‑150 (MQ=255)
                                                                   gCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGtat                                                                                >  1:933290/1‑150 (MQ=255)
                                                                    cAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGtatt                                                                               >  1:1245171/1‑150 (MQ=255)
                                                                     aaaaTATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGtatta                                                                              <  2:109727/150‑1 (MQ=255)
                                                                     aaaaTATTGTTCTGTTCTGGTACCAGTGGGATCCAATGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGtatta                                                                              <  2:104177/150‑1 (MQ=255)
                                                                          atTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTc                                                                         >  2:1272661/1‑150 (MQ=255)
                                                                           ttGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCt                                                                        <  1:106999/150‑1 (MQ=255)
                                                                              ttctgttctgGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCa                                                                     <  1:67879/150‑1 (MQ=255)
                                                                               tctgttctgGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCAt                                                                    <  2:1434013/150‑1 (MQ=255)
                                                                                 tgttctgGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATcc                                                                  >  1:830155/1‑150 (MQ=255)
                                                                                    tctgGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCa                                                               >  2:389832/1‑150 (MQ=255)
                                                                                     ctgGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCaa                                                              <  2:89636/150‑1 (MQ=255)
                                                                                      tgGTACCAGTGGGATCCCCTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATACTCAAt                                                             >  2:1363379/1‑150 (MQ=255)
                                                                                      tgGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAAt                                                             <  1:579629/150‑1 (MQ=255)
                                                                                       ggTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATa                                                            >  2:1404834/1‑150 (MQ=255)
                                                                                       ggTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATa                                                            >  2:475692/1‑150 (MQ=255)
                                                                                       ggTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATa                                                            <  2:447662/150‑1 (MQ=255)
                                                                                         tACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATAcc                                                          <  2:1464368/150‑1 (MQ=255)
                                                                                            cAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGtt                                                       <  1:1420223/150‑1 (MQ=255)
                                                                                            cAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTACAACAAAGTTTTATTTTCTGCAGCGTCAATCCCGta                                                       >  2:489410/1‑149 (MQ=255)
                                                                                              gTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAGTACCGTTaa                                                     >  2:526421/1‑150 (MQ=255)
                                                                                               tGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTaaa                                                    <  1:1341360/150‑1 (MQ=255)
                                                                                                gggATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAAt                                                   >  1:1439378/1‑150 (MQ=255)
                                                                                                gggATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAAt                                                   <  1:489968/150‑1 (MQ=255)
                                                                                                 ggATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATg                                                  <  2:1439378/150‑1 (MQ=255)
                                                                                                  gATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGc                                                 >  1:1185990/1‑150 (MQ=255)
                                                                                                    tCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCaa                                               >  1:681345/1‑150 (MQ=255)
                                                                                                    tCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCaa                                               >  2:589614/1‑150 (MQ=255)
                                                                                                     ccACTGCATCCCGTCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGTCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAAc                                              <  1:200487/150‑1 (MQ=255)
                                                                                                     ccACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTcaac                                                                                       >  2:486265/1‑109 (MQ=255)
                                                                                                     ccACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTcaac                                                                                       <  1:486265/109‑1 (MQ=255)
                                                                                                      cACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAAAGCAAcc                                             >  1:1182001/1‑150 (MQ=255)
                                                                                                       aCTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATCCCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAAccc                                            >  1:759800/1‑150 (MQ=255)
                                                                                                       aCTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAAccc                                            <  2:971708/150‑1 (MQ=255)
                                                                                                           cATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAAc                                        <  1:1053589/150‑1 (MQ=255)
                                                                                                             tCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAAccc                                      >  2:75527/1‑150 (MQ=255)
                                                                                                                      taatttCATTGAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTcc                             <  1:1429920/145‑1 (MQ=255)
                                                                                                                       aaGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTccc                            <  2:933290/150‑1 (MQ=255)
                                                                                                                       aaGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTccc                            <  1:475692/150‑1 (MQ=255)
                                                                                                                        aGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCt                           >  2:510232/1‑150 (MQ=255)
                                                                                                                           tCATTCAGTTTCACCGGACGCATATTTGCGCCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTtg                        >  1:377470/1‑150 (MQ=255)
                                                                                                                            cATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTtgc                       <  2:1376330/150‑1 (MQ=255)
                                                                                                                            cATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTtgc                       <  1:840602/150‑1 (MQ=255)
                                                                                                                                  gTTTCACCGGACGCATATTTGCGTCCAATATTTCTATCCCCAAACCGGAGGCCGGATTGCTTCCGTCCTCCAGTTTCAACAAAGTATTATTTTCTGCCTCCTCAATACCGTTAAATGCCACCCGAACCCCCGTTGGCCCTTTGCTGCAtt                 >  2:1471503/1‑150 (MQ=255)
                                                                                                                                  gTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCAtt                 >  1:1296026/1‑150 (MQ=255)
                                                                                                                                   tttCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTc                >  1:1496504/1‑150 (MQ=255)
                                                                                                                                      tattggACGCAAATATGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCGTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACt             >  1:782662/5‑150 (MQ=255)
                                                                                                                                      cACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCACTTTGCTGCATTCACt             <  1:190846/150‑1 (MQ=255)
                                                                                                                                        ttggACGCAAATATGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTa           >  1:1437798/3‑150 (MQ=255)
                                                                                                                                         tggACGCAAATATGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTaa          >  2:359981/2‑150 (MQ=255)
                                                                                                                                          ggACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAAc         >  2:750771/1‑150 (MQ=255)
                                                                                                                                             cGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTa      >  1:562840/1‑150 (MQ=255)
                                                                                                                                             cGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTa      <  2:1128431/150‑1 (MQ=255)
                                                                                                                                             cGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTa      >  2:1238867/1‑150 (MQ=255)
                                                                                                                                             cGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTa      >  1:172521/1‑150 (MQ=255)
                                                                                                                                                 tatTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCt  >  1:703749/1‑150 (MQ=255)
                                                                                                                                                 tatTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATACTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCt  >  1:883364/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GTTGCCGAAGCCCTCACCAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATCT  >  NC_000913/1588331‑1588625

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: ‑