Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 121,012 G→A 5.5% G180G (GGC→GGT aroP ← aromatic amino acid transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913121,0120GA 5.5% 192.4 / 2.5 73G180G (GGC→GGTaroParomatic amino acid transporter
Reads supporting (aligned to +/- strand):  ref base G (35/33);  new base A (2/2);  total (37/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.88e-01

TATCAGCTTCTGCTGCGGTGATCCCCACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCATCTCGCCAAACACTTTAACGTTGGTCAGGTTGATGGCGTTAATCACCACAAAGAATACGGCGG  >  NC_000913/120867‑121160
                                                                                                                                                 |                                                                                                                                                    
tATCAGCTTCTGCTGCGGTGATCCCCACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGccgccg                                                                                                                                                   >  1:18934/1‑149 (MQ=255)
 aTCAGCTTCTGCTGCGGTGATCCCCACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGCTTGCTAACGGCGGCCCGCAGGCGACCGt                                                                                                                                                  <  1:1003231/149‑1 (MQ=255)
   cAGCTTCTGCTGCGGTGATCCCCACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTg                                                                                                                                                <  2:357979/149‑1 (MQ=255)
     gCTTCTGCTGCGGTGATCCCCACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGcc                                                                                                                                              <  2:618872/149‑1 (MQ=255)
      cctCTGCTGCGCTGATCCCCACCAGTTCCAGCCCACCTAACGAGACCATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCa                                                                                                                                             <  2:820314/147‑1 (MQ=255)
      cTTCTGCTGCGGTGATCCCCACCAGTTCCAGACCACCGAACTAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCa                                                                                                                                             <  1:1163533/149‑1 (MQ=255)
             tgcGGTGATCCCCACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATa                                                                                                                                      <  1:517616/149‑1 (MQ=255)
              gcGGTGATCCCCACCAGTTCAAGACCACCGAACGCGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACAAACCTGATCCCACAGGTTGCTAACGGACGCCTGCGGGCCGCCGTGGCCCCTGAATag                                                                                                                                     >  2:1100286/1‑149 (MQ=255)
               cGGTGATCCCCACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATagc                                                                                                                                    >  2:846599/1‑149 (MQ=255)
                ggTGATCCCCACCAGTGCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATagca                                                                                                                                   <  2:1096306/149‑1 (MQ=255)
                  tGATCCCCACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATAATCACCAGCCCGTTGAAGCAGTGCGGCATGAAACCACCCTGATCACAAAGGTTGCTAAAGGTCGCCTGCTGGCCGCCCTTGCCTCTGAATagcagc                                                                                                                                 >  2:85286/1‑149 (MQ=255)
                     tCCCCACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGcagccag                                                                                                                              >  1:195312/1‑149 (MQ=255)
                     tCCCCACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGcagccag                                                                                                                              >  2:112515/1‑149 (MQ=255)
                         caccaGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGtt                                                                                                                                                 >  1:363042/1‑126 (MQ=255)
                         caccaGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGtt                                                                                                                                                 <  2:363042/126‑1 (MQ=255)
                          accaGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGCAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAgccgcc                                                                                                                         >  1:122820/1‑149 (MQ=255)
                               ttCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCTGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCCCTGCATAGCAGCCAGCCGCCGAAga                                                                                                                    >  1:182416/1‑149 (MQ=255)
                                tCCAGACCACCGAACGAGCACATGATAATCGCCCTCCTCCTCCCCCGCCAGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAgat                                                                                                                   <  2:958590/149‑1 (MQ=255)
                                  ctgACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAgatga                                                                                                                 <  1:1129035/147‑1 (MQ=255)
                                   aGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCCCAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAgatgat                                                                                                                >  2:1080631/1‑149 (MQ=255)
                                       caccGAACTAGAACATGATTATCGCCATCATCATCACCAGCCCGGGGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACTGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGatcatc                                                                                                            <  2:1108257/149‑1 (MQ=255)
                                       caccGAACGAGAACATGATAATCGCCCTCATCATCACCCGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGatcatc                                                                                                            <  1:105933/149‑1 (MQ=255)
                                          cGACCGAGAACATGATAATCTCCATCATCATCCCCAGCCCGGTGAAGCCGAGCGGCAGGAAACCACCCTGATCCCAGAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCt                                                                                                         <  2:73794/149‑1 (MQ=255)
                                            aaCGAGAACATGATCAACGCCATCCTCATCACCAGCCCGGTGAAGCCGTGCGGCAGGCACCCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTac                                                                                                       <  1:695278/149‑1 (MQ=255)
                                            aaCGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGATGCCGGGCGGCAGGCAACCACCCTGATCCCACAGGTTCCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCGAGCCGCCGTAGATGATCATCGCTac                                                                                                       <  1:942314/149‑1 (MQ=255)
                                               gagaACATGATAATCGCCATCATCCTCCCCAGCCCGGTGAAGACGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCCTTGCCCCTGAATAGCACCCAGCCGCCGAAGATGATCATCGCTaccac                                                                                                    >  1:34136/1‑149 (MQ=255)
                                                 gaACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCCCCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACcg                                                                                                  >  1:859757/1‑149 (MQ=255)
                                                  aaCATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACcgc                                                                                                 <  1:350977/149‑1 (MQ=255)
                                                     aTGAGACTCGCCAGCATCATCACAAGCCCGGTGCAGCCGTGCGGCAGGACACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCCCCGCGAt                                                                                              <  2:667543/149‑1 (MQ=255)
                                                       gATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATaa                                                                                            >  1:259548/1‑149 (MQ=255)
                                                       gATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACCGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCAGCGCGATaa                                                                                            >  2:105306/1‑149 (MQ=255)
                                                           aTCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACttt                                                                                        <  2:1055627/149‑1 (MQ=255)
                                                            tCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACGGATTAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTg                                                                                       >  2:541957/1‑149 (MQ=255)
                                                            tCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCACCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTg                                                                                       >  2:591369/1‑149 (MQ=255)
                                                               ccATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCCGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATa                                                                                    >  1:199069/1‑149 (MQ=255)
                                                                catcatcatCACCAGCCCGGTGAAGCCGTGCGGCTGGAAACCACCCAGATCCCACAGCTTGCTAACGGTCGCCTGCGGGCCGCCGTTGGCACTGAATAGCAGCCAGCCGCCGAAGGTGCTCATCGCTACGCCAGCGATAACTTAGATaa                                                                                   >  2:687351/1‑149 (MQ=255)
                                                                catcatcatCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATaa                                                                                   <  1:739022/149‑1 (MQ=255)
                                                                 atcatcatcACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAAt                                                                                  <  1:885034/149‑1 (MQ=255)
                                                                   catcatcaCCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATgg                                                                                <  2:18934/149‑1 (MQ=255)
                                                                     tcatcaCCAGCCCGTTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTAGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATggaa                                                                              >  2:913330/1‑147 (MQ=255)
                                                                        tcaCCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAAc                                                                           >  2:1155402/1‑149 (MQ=255)
                                                                                ccGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTc                                                                   >  1:276619/1‑149 (MQ=255)
                                                                                       aGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCCCTGAATAGCAGCCAGCCGCCGACGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAACCCAGAACTCCAGCTCg                                                            >  1:1079516/1‑149 (MQ=255)
                                                                                       aGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCGGAAGATGATCATCGCTACCGCCGCGAGAACTTTGAGAATGGCATCCCAGACCTCCATCTCg                                                            >  1:494154/1‑149 (MQ=255)
                                                                                            tGCGGCAGGAAACCACCCTGATCACACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCAGTGAATAGCAGCCAGCCGCCGAAGATGATCATCTCTTCCACCGCGATTACTTTGTTAATGCCAAACCAGAACTCCATCTCGCCaaa                                                       >  1:306658/1‑149 (MQ=255)
                                                                                                gCAGGAACCCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGATGATGATCATCGCTCCCACCGGGATAACTTTGATAATGGCCAACCAGCAGTCCATCTGGCCAAACACt                                                   >  2:70439/1‑149 (MQ=255)
                                                                                                 cAGGAAACCACGCTGATACCACCGGTGGCTAACGGTCGCCGGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGAGGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCATCTCGCCAAACACtt                                                  <  2:642490/149‑1 (MQ=255)
                                                                                                  aGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCATCTCGCCAAACACttt                                                 >  2:467000/1‑149 (MQ=255)
                                                                                                    gAAACCAGGCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCAGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCTTCGCTACCTCGGCGATAACTTTGATTTTGGCATACCACAACTGCATCTCGACAAACACTTTaa                                               <  2:751312/149‑1 (MQ=255)
                                                                                                    gAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCCCCGCGCTAACTTTGATAATGGCACACCCGACCGCCCTGTCGCCGAACCCTTTaa                                               >  1:194321/1‑149 (MQ=255)
                                                                                                         caccCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTCAATAGCAGCCAGCAGCCGAAGATGATCATCGCTACCACCCCGATAACATTGATAATGTCAAAACAGAACTCCAACTCGCCAAACACTTTAACGtgt                                          >  2:444891/1‑147 (MQ=255)
                                                                                                           cccTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATCGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCATCTCGCCAAACACTTTAACGTCGGt                                        <  2:80288/149‑1 (MQ=255)
                                                                                                            ccTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGACTCGCAGCCAGCCGCCGAAGATGATCATCGCTTGCACCGCGATAACTTTGATATTGGCAACCCAGAACTCCATCTCGCCAGACTCTTTAACGTTGGTc                                       >  1:1145371/1‑149 (MQ=255)
                                                                                                            ccTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCATCTCGCCAAACACTTTAACGTTGGTc                                       >  1:320252/1‑149 (MQ=255)
                                                                                                               gATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCATCTCGCCAAACACTTTAACGTTGGTCAgg                                    <  1:846599/149‑1 (MQ=255)
                                                                                                                  cccACAGGTTGCTAACGGACGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCGACCACCGCGATAACTTTGATAATGGCACACCAGCACTCCATCTCGCGAAACGCTTTAACGTTGGTCAGGTTg                                 >  2:431542/1‑149 (MQ=255)
                                                                                                                   ccACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGCTAACTTTGATAATGGCAAACCAGAACTCCATCTCGCCAAACACTTTAACGTTGGTCAGGTTGa                                >  1:377143/1‑149 (MQ=255)
                                                                                                                   ccACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCATCTCGCCAAACACtt                                                  >  2:577505/1‑131 (MQ=255)
                                                                                                                   ccACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCATCTCGCCAAACACtt                                                  <  1:577505/131‑1 (MQ=255)
                                                                                                                       aGGTTGCTCACGGTCGCCTGCGTGCCTCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCATCTCGCCAAACTCTTTAACGTTTGTCAGGTTGATGGc                            <  2:527482/149‑1 (MQ=255)
                                                                                                                         gTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCATCTCGCAAAACACTTTAACGTTGGTCAGGTTGATGGCGt                          >  2:311812/1‑149 (MQ=255)
                                                                                                                         gTTGCTAACGGTCGCCTGCGGGACGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCATCTCGCCAAACACTTTAACGTTGGTCAGGTTGATGGCGt                          <  2:696379/149‑1 (MQ=255)
                                                                                                                            ggcaaCGGCGGCCCGCAGGCGACCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACAACCGCGATAACTTTGATAATGGCAAACCAGAACGCCAGGTCGCCAAACACGTTCACGTTGGTCAGGTGGATGGCGTTaa                       >  1:500000/4‑149 (MQ=255)
                                                                                                                             cTAACGGTCGCCTGCGGGCCGCCGTTGCCACTTAATATCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCATCTCGCCAAACACTTTAACGTTGGTCAGGTTGATGGCGTTAAt                      <  1:577396/149‑1 (MQ=255)
                                                                                                                             cTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCATCTCGCCAAACACTTTAACGTTGGTCAGGTTGATGGCGTTAAt                      <  2:122820/149‑1 (MQ=255)
                                                                                                                              tAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCATCTCGCCAAACACTTTAACGTTGGTCAGGTTGATGGCGTTAATc                     >  2:279539/1‑149 (MQ=255)
                                                                                                                              tAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCACAACTCAATCTCCACAAACACTTAAACGTTGGTCAGGTATATGGCGTTAATc                     >  2:4376/1‑149 (MQ=255)
                                                                                                                                aCGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCATCTCGCCAAACACTTTAAGGTTGGTCAGGTTGATGGCGTTAATcac                   <  2:561744/149‑1 (MQ=255)
                                                                                                                                aCGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCATCTCGCCAAACACTTTAACGTTGGTCAGGTTGATGGCGTTAATcac                   >  1:317217/1‑149 (MQ=255)
                                                                                                                                aCGGTCGCCTGCGGGCCACCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCATCTCGCCAAACACTTTAACGTTGGTCAGGTTGATGGCGTTAATcac                   <  1:591369/149‑1 (MQ=255)
                                                                                                                                   gTCGCCTGCGGGCCGCCGTTGCCACTTAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCATCTCGCCAAACACTTTAACGTTGGTCAGGTTGATGGCGTTAATcaccac                <  2:259548/149‑1 (MQ=255)
                                                                                                                                        cTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCATCTCGCCAAACACTTTAACGTTGGTCAGGTTGATGGCGTTAATCACCACAAAGa           <  1:912495/149‑1 (MQ=255)
                                                                                                                                         tGCGGGCCGCAGTTGCCAATGAATAGCAGCCCGCCGCCGAAGATGAGCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGATCTCCATCTCGCCAAACACTTTAACGTTGGTCAGGTTGATGGCGTTAATCACCACAAAGaa          <  2:777355/149‑1 (MQ=255)
                                                                                                                                         tGCGGGACGCCGGTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATACTGGCAAACCAGAACTCCATCTCGCCAAACACTTTAACGTTGGTCAGGTTGATGGCGTTAATCACCACAAAGaa          <  1:1039452/149‑1 (MQ=255)
                                                                                                                                               ccgccgTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCATCTCGCCAAACACTTTAACGTTGGTCGGGTTGATGGCGTTAATCACCACAAAGAATAcggc    >  1:632178/1‑149 (MQ=255)
                                                                                                                                                 gccgTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCCCCGCGATAACTTTGATAATGGCAAACCAGAACTCCATCTCGCCAACCACTTTAACGTTGGTGCGGTTGATGGCGTTAATCACCACAAAGAATAcggcgg  >  2:1039117/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
TATCAGCTTCTGCTGCGGTGATCCCCACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCATCTCGCCAAACACTTTAACGTTGGTCAGGTTGATGGCGTTAATCACCACAAAGAATACGGCGG  >  NC_000913/120867‑121160

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: