Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,618,945 A→T 5.6% F406I (TTC→ATC)  uraA ← uracil permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,618,9450AT 5.6% 618.8 / 19.5 179F406I (TTC→ATC) uraAuracil permease
Reads supporting (aligned to +/- strand):  ref base A (77/92);  new base T (2/8);  total (79/100)
Fisher's exact test for biased strand distribution p-value = 1.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.96e-01

Only 100 of 184 total aligned reads displayed.

TCACAAAATCGCGAGAAGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCTGTGCTTTGTTGTAATCCACTTTC  >  NC_000913/2618796‑2619093
                                                                                                                                                     |                                                                                                                                                     
tCACAAAATCGCGAGAAGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGaa                                                                                                                                                       <  1:4118823/150‑1 (MQ=255)
  acaAAATCGCGAGAAGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGGGCAGCACACTGATCAGTTTGAAgt                                                                                                                                                     >  1:3854009/1‑149 (MQ=255)
  acaAAATCGCGAGAAGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGa                                                                                                                                                     >  2:2912400/1‑150 (MQ=255)
  acaAAATCGCGAGAAGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGa                                                                                                                                                     >  1:2037588/1‑150 (MQ=255)
  acaAAATCGCGAGAAGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGa                                                                                                                                                     >  1:2037617/1‑150 (MQ=255)
   caAAATCGCGAGAAGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGAt                                                                                                                                                    <  2:3740463/150‑1 (MQ=255)
    aaaaTCGCGAGAAGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATc                                                                                                                                                   <  1:924401/150‑1 (MQ=255)
     aaaTCGCGAGAAGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCa                                                                                                                                                  <  2:112569/150‑1 (MQ=255)
     aaaTCGCGAGAAGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCa                                                                                                                                                  <  2:1762176/150‑1 (MQ=255)
        tCGCGAGAAGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAAc                                                                                                                                               <  1:1347661/150‑1 (MQ=255)
            gagaAGTTTACCTCAGTATCCCATTATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCTCGTCTTCTGCGTCCATCTACACTTCTTCAGGGCGCAGCACACTGATCAGTTTGATCTTCAAACtgat                                                                                                                                           <  1:489410/150‑4 (MQ=255)
            gagaAGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATCTCCGCGTCTTCTGCGTCCAGCACCTCTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTaa                                                                                                                                           <  2:4129849/150‑1 (MQ=255)
            gagaAGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGATCTTCAAACtgat                                                                                                                                           <  1:83921/150‑4 (MQ=255)
            gagaAGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTaa                                                                                                                                           <  1:697897/150‑1 (MQ=255)
            gagaAGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTaa                                                                                                                                           <  1:2781960/150‑1 (MQ=255)
            gagaAGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTAGATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGATCTTCAAACtgat                                                                                                                                           <  1:1578509/150‑4 (MQ=255)
              gaAGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGc                                                                                                                                         <  1:3510538/150‑1 (MQ=255)
                aGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCg                                                                                                                                       >  1:2138720/1‑150 (MQ=255)
                aGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCg                                                                                                                                       >  1:3792658/1‑150 (MQ=255)
                  tttACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGAt                                                                                                                                     >  1:1006623/1‑150 (MQ=255)
                   ttACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATa                                                                                                                                    <  1:3496340/150‑1 (MQ=255)
                   ttACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATa                                                                                                                                    <  2:2592377/150‑1 (MQ=255)
                       cacaGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACcga                                                                                                                                <  2:2646748/150‑1 (MQ=255)
                         caGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACcgacg                                                                                                                              >  2:1222243/1‑150 (MQ=255)
                         caGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACcgacg                                                                                                                              >  2:1222189/1‑150 (MQ=255)
                          aGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACcgacga                                                                                                                             <  2:1483962/150‑1 (MQ=255)
                          aGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACcgacga                                                                                                                             <  2:1015500/150‑1 (MQ=255)
                          aGAAGCCCATGATGCTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTTTTTGACTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTACTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACcgacga                                                                                                                             <  2:1714351/150‑1 (MQ=255)
                                 cATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGc                                                                                                                      >  1:2197550/1‑150 (MQ=255)
                                    gatgTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTTTCTGTTATATAAGCGTCTTATGCGTACAGCACCACATCTTAAGGGCGCAGCACCAAGATAATTTTGAAGATAAAAATTAAGCCGATACCGAAAATGTTCGCCac                                                                                                                   >  1:354020/1‑149 (MQ=255)
                                      tgTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAgcg                                                                                                                 <  2:1314312/150‑1 (MQ=255)
                                      tgTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAgcg                                                                                                                 <  2:1852628/150‑1 (MQ=255)
                                       gTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAgcgc                                                                                                                >  1:2153312/1‑150 (MQ=255)
                                         ttGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCa                                                                                                              <  2:3321115/150‑1 (MQ=255)
                                         ttGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCa                                                                                                              <  1:3691841/150‑1 (MQ=255)
                                           gAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATa                                                                                                            <  2:3595799/150‑1 (MQ=255)
                                             aCCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATAcc                                                                                                          >  1:3018067/1‑150 (MQ=255)
                                                gggCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCttt                                                                                                      <  1:2451514/150‑1 (MQ=255)
                                                 ggCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCtttt                                                                                                     <  2:703318/150‑1 (MQ=255)
                                                 ggCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCtttt                                                                                                     <  1:1025133/150‑1 (MQ=255)
                                                  gcagcaACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGAGGGTCGCCAGCGCCATACCTTTTa                                                                                                    <  2:216005/150‑1 (MQ=255)
                                                  gcagcaACACTGCCCGGTCGGGACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTa                                                                                                    >  1:3584189/1‑150 (MQ=255)
                                                  gcagcaACACTGCCAGGTCGGTACCATATTTGTGTGTTATATCCGCGTATTCAGCGTCCAGCACCACTTCTTCAGGGCGCCGCACACTGATCAGTTTGAAGATAAAACTTAAGCCGATACCCACGATGGTAGCCAGCGCCAAACCTtata                                                                                                    >  1:3584225/1‑150 (MQ=255)
                                                                      gtACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGttt                                                                                >  1:3975829/1‑150 (MQ=255)
                                                                            tattTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTc                                                                          >  1:1912946/1‑150 (MQ=255)
                                                                            tattTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTc                                                                          <  2:3707372/150‑1 (MQ=255)
                                                                             attTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTcg                                                                         <  2:180612/150‑1 (MQ=255)
                                                                              cttGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTcgc                                                                        <  1:95381/149‑1 (MQ=255)
                                                                                tGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCgcgc                                                                      <  1:158594/150‑1 (MQ=255)
                                                                                tGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTGCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCgcgc                                                                      <  2:2853813/150‑1 (MQ=255)
                                                                                  tCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCa                                                                    <  1:3344394/150‑1 (MQ=255)
                                                                                   cTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCCGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCAc                                                                   <  1:1133470/150‑1 (MQ=255)
                                                                                    tGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACt                                                                  >  2:2505394/1‑150 (MQ=255)
                                                                                      ttATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGa                                                                <  2:3406219/150‑1 (MQ=255)
                                                                                      ttATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGa                                                                <  2:141700/150‑1 (MQ=255)
                                                                                      ttATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCAGATGTTAACCTACGCGCCACTGa                                                                >  2:1463700/1‑150 (MQ=255)
                                                                                       tataTCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGAc                                                               <  2:2731129/150‑1 (MQ=255)
                                                                                                 tcttctGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATc                                                     <  2:3850278/150‑1 (MQ=255)
                                                                                                   ttctGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCaa                                                   >  1:3900579/1‑150 (MQ=255)
                                                                                                     ctGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCaaaa                                                 >  1:2741053/1‑150 (MQ=255)
                                                                                                     ctGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCaaaa                                                 >  2:1940320/1‑150 (MQ=255)
                                                                                                      tGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAAt                                                <  2:3062822/150‑1 (MQ=255)
                                                                                                      tGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAAt                                                >  1:3768965/1‑150 (MQ=255)
                                                                                                        cGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCa                                              <  2:964276/150‑1 (MQ=255)
                                                                                                        cGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCa                                              <  2:1021049/150‑1 (MQ=255)
                                                                                                             aGCACCACTTCTTCAGGGCGCAGCACTCTGATCTGTTTGTAGATCAAAGTTATGCCGTTAGCGACGATGGTCGCCAGCGCCATAGCTTTTATCTCCGCCGCGCCGTTGTTTTCCTTCGCGCCTCTGACGGCGATGTTCTTATTCTCGGaa                                         <  2:86254/150‑1 (MQ=255)
                                                                                                             aGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCTCCAGCGCCATACCATTTAAACACCCCCGCCCCAATTTTAACGTACGCCCCACTGACCCCAAGAAGCAAAATCACGGa                                          >  1:2484300/1‑150 (MQ=255)
                                                                                                             aGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGaa                                         >  1:2688114/1‑150 (MQ=255)
                                                                                                             aGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGaa                                         >  1:2574449/1‑150 (MQ=255)
                                                                                                                accacTTCTTCCGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTc                                      <  1:2443773/150‑1 (MQ=255)
                                                                                                                 ccacTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCa                                     <  2:2257140/150‑1 (MQ=255)
                                                                                                                 ccacTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAAACACGGAAGTCa                                     >  2:3101924/1‑150 (MQ=255)
                                                                                                                  cacTTCTTCAGGGCGCAACACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAg                                    >  2:391781/1‑150 (MQ=255)
                                                                                                                   acTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAgg                                   >  1:1207140/1‑150 (MQ=255)
                                                                                                                    cttcttCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGa                                  >  1:445222/1‑150 (MQ=255)
                                                                                                                     ttcttcAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTGTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGag                                 >  2:2396711/1‑149 (MQ=255)
                                                                                                                        ttcAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCAGGGAAGTCAGGATCAg                              >  2:730719/1‑150 (MQ=255)
                                                                                                                        ttcAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAg                              >  1:2512902/1‑150 (MQ=255)
                                                                                                                        ttcAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAg                              >  1:3528342/1‑150 (MQ=255)
                                                                                                                        ttcAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAg                              >  2:40732/1‑150 (MQ=255)
                                                                                                                        ttcAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAg                              >  1:3188374/1‑150 (MQ=255)
                                                                                                                        ttcAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAg                              >  2:106336/1‑150 (MQ=255)
                                                                                                                         tcAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAgg                             >  2:3857133/1‑150 (MQ=255)
                                                                                                                         tcAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCAGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAgg                             >  2:208552/1‑150 (MQ=255)
                                                                                                                               cgcAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCtgt                       <  1:2369493/150‑1 (MQ=255)
                                                                                                                               cgcAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCtgt                       <  1:1694359/150‑1 (MQ=255)
                                                                                                                               cgcAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCtgt                       <  1:1694634/150‑1 (MQ=255)
                                                                                                                                gcagcaCACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGAGGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCtgtg                      <  2:3284587/150‑1 (MQ=255)
                                                                                                                                 cagcaCACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCTGTGc                     <  1:2780412/150‑1 (MQ=255)
                                                                                                                                 cagcaCACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCTGTGc                     <  2:2700614/150‑1 (MQ=255)
                                                                                                                                   gcaCACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCTGTGCtt                   >  1:1251180/1‑150 (MQ=255)
                                                                                                                                    cacaCTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCTGTGCttt                  <  1:1501792/150‑1 (MQ=255)
                                                                                                                                       acTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCTGTGCTttgtt               >  2:1297960/1‑150 (MQ=255)
                                                                                                                                       acTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCTGTGCTttgtt               >  2:217563/1‑150 (MQ=255)
                                                                                                                                          cttaaGTTTGATCTTCAAACTTAAGCCGATACCGACGATGGTCGCAAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCTGTGCTTTGTTGTa            >  2:4043879/5‑150 (MQ=255)
                                                                                                                                            tCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCTGTGCTTTGTTGTAAt          <  2:3404868/150‑1 (MQ=255)
                                                                                                                                                 ttGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCTGTGCTTTGTTGTAATCCACt     >  1:2966780/1‑150 (MQ=255)
                                                                                                                                                  tgtaGATCAAAGTTAAGCCGATACCGACGTTGGTCGCCTGCGCCATACGTTTTAACTGGGCCGCGCCGATGTTTACCTTCGCGCGACTGACGCCGATGATCAAATTCTGGGTTGTCAGGATCTGGTTCTGTGCTTTGTTGTATTCCAGtt    >  2:3097507/4‑150 (MQ=255)
                                                                                                                                                  tGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCTGTGCTTTGTTGTAATCCACtt    >  2:1665114/1‑150 (MQ=255)
                                                                                                                                                    aaGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCTGTGCTTTGTTGTAATCCACTttc  <  2:3975829/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCACAAAATCGCGAGAAGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACCTTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCTGTGCTTTGTTGTAATCCACTTTC  >  NC_000913/2618796‑2619093

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