Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,302,827 G→A 6.1% intergenic (+14/‑72) oppA → / → oppB oligopeptide ABC transporter periplasmic binding protein/oligopeptide ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,302,8270GA 6.1% 265.7 / 5.2 81intergenic (+14/‑72)oppA/oppBoligopeptide ABC transporter periplasmic binding protein/oligopeptide ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (35/41);  new base A (3/2);  total (38/43)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

TAAGGATTCGGCCATTGTTCCTGTTTATTACTACGTGAATGCGCGTCTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCTTATTACTATTTCGT  >  NC_000913/1302681‑1302968
                                                                                                                                                  |                                                                                                                                             
tAAGGATTCGGCCATTGTTCCTGTTTATTACTACGTGAATGCGCGTCTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGc                                                                                                                                             <  3:1046032/149‑1 (MQ=255)
 aaGGATTCGGCCATTGTTCCTGTTTATTACTACGTGAATGCGCGTCTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTCATGGCAATACGTGGAGCa                                                                                                                                            >  3:845390/1‑149 (MQ=255)
  aGGATTCGGCCATTGTTCCTGTTTATTACTACGTGAATGCGCGTCTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAg                                                                                                                                           <  1:399212/149‑1 (MQ=255)
    gATTCGGCCATTGTTCCTGTTTATTACTACGTGAATGCGCGTCTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGa                                                                                                                                         <  3:75991/149‑1 (MQ=255)
    gATTCGGCCATTGTTCCTGTTTATTACTACGTGAATGCGCGTCTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGa                                                                                                                                         <  2:1125560/149‑1 (MQ=255)
         ggCCATTGTTCCTGTTTATTACTACGTGAATGCGCGTCTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTc                                                                                                                                    <  4:207644/149‑1 (MQ=255)
         ggCCATTGTTCCTGTTTATTACTACGTGAATGCGCGTCTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTc                                                                                                                                    >  3:556610/1‑149 (MQ=255)
           ccATTGTTCCTGTTTATTACTACGTGAATGCGCGTCTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCt                                                                                                                                  <  3:1319742/149‑1 (MQ=255)
            cATTGTTCCTGTTTATTACTACGTGAATGCGCGTCTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCt                                                                                                                                  >  1:830946/1‑148 (MQ=255)
            cATTGTTCCTGTTTATTACTACGTGAATGCGCGTCTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCt                                                                                                                                  <  2:830946/148‑1 (MQ=255)
            cATTGTTCCTGTTTATTACTACGTGAATGCGCGTCTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTg                                                                                                                                 >  2:494082/1‑149 (MQ=255)
             aTTGTTCCTGTTTATTACTACGTGAATGCGCGTCTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGc                                                                                                                                <  4:82379/149‑1 (MQ=255)
             aTTGTTCCTGTTTATTACTACGTGAATGCGCGTCTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGc                                                                                                                                >  1:664880/1‑149 (MQ=255)
             aTTGTTCCTGTTTATTACTACGTGAATGCGCGTCTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGc                                                                                                                                >  2:976725/1‑149 (MQ=255)
               tGTTCCTGTTTATTACTACGTGAATGCGCGTCTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTc                                                                                                                              <  2:664880/149‑1 (MQ=255)
               tGTTCCTGTTTATTACTACGTGAATGCGCGTCTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGAGCAGGACACTCCTGccc                                                                                                                              <  4:845390/149‑3 (MQ=255)
                     tGTTTATTACTACGTGAATGCGCGTCTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGgt                                                                                                                        >  2:1179031/1‑149 (MQ=255)
                       tttATTACTACGTGAATGCGCGTCTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGgtgt                                                                                                                      <  1:1420356/149‑1 (MQ=255)
                             actacGTGAATGCGCGTCTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGAtt                                                                                                                <  1:976725/149‑1 (MQ=255)
                              ctacGTGAATGCGCGTCTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGAttt                                                                                                               >  4:789458/1‑149 (MQ=255)
                               tacGTGAATGCGCGTCTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAgt                                                                                                                                       <  4:931315/124‑1 (MQ=255)
                               tacGTGAATGCGCGTCTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGAtttt                                                                                                              >  3:1151383/1‑149 (MQ=255)
                               tacGTGAATGCGCGTCTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGACAATACGTGGGGCAGGAgt                                                                                                                                       >  3:931315/1‑124 (MQ=255)
                                acGTGAATGCGCGTCTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGAttttt                                                                                                             >  3:1026172/1‑149 (MQ=255)
                                  gTGAATGCGCGTCTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGAttttt                                                                                                             >  3:356835/1‑147 (MQ=255)
                                  gTGAATGCGCGTCTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGAttttt                                                                                                             <  4:356835/147‑1 (MQ=255)
                                    gAATGCGCGTCTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCg                                                                                                         <  1:494082/149‑1 (MQ=255)
                                         cgcgTCTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTa                                                                                                    <  4:1026172/149‑1 (MQ=255)
                                          gcgTCTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTAc                                                                                                   >  4:1099394/1‑149 (MQ=255)
                                               tGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAg                                                                                              >  3:827732/1‑149 (MQ=255)
                                                ggTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAgg                                                                                             <  1:1179031/149‑1 (MQ=255)
                                                  tGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGca                                                                                           >  1:1074649/1‑149 (MQ=255)
                                                    aaaCCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTggg                                                                                                                                               >  4:1402613/1‑95 (MQ=255)
                                                    aaaCCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTggg                                                                                                                                               <  3:1402613/95‑1 (MQ=255)
                                                    aaaCCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGcaca                                                                                         >  2:49426/1‑149 (MQ=255)
                                                      aCCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAgg                                                                                       >  3:1245165/1‑149 (MQ=255)
                                                        cGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGa                                                                                                                  <  1:918297/120‑1 (MQ=255)
                                                        cGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGa                                                                                                                  >  2:918297/1‑120 (MQ=255)
                                                                 gtggCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAg                                                                            <  4:556610/149‑1 (MQ=255)
                                                                  tggCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGAtt                                                                                                                >  3:727570/1‑112 (MQ=255)
                                                                  tggCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGAtt                                                                                                                <  4:727570/112‑1 (MQ=255)
                                                                    gCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGa                                                                                                                                         >  1:190340/1‑85 (MQ=255)
                                                                    gCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGa                                                                                                                                         <  2:190340/85‑1 (MQ=255)
                                                                    gCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAg                                                                            >  3:1154483/1‑146 (MQ=255)
                                                                    gCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAg                                                                            <  4:1154483/146‑1 (MQ=255)
                                                                        taCCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTa                                                                             >  1:288404/1‑141 (MQ=255)
                                                                        taCCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTa                                                                             <  2:288404/141‑1 (MQ=255)
                                                                             gCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGAGCAGGACACTCCTGcccc                                                                                                                             >  4:702114/1‑88 (MQ=38)
                                                                              cAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTaaaa                                                               <  4:1151383/149‑1 (MQ=255)
                                                                               aaaGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAAt                                                              <  2:1074649/149‑1 (MQ=255)
                                                                                aaGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAg                                                                            >  3:418234/1‑134 (MQ=255)
                                                                                aaGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAg                                                                            <  4:418234/134‑1 (MQ=255)
                                                                                aaGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAAtt                                                             <  3:1099394/149‑1 (MQ=255)
                                                                                  gATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCg                                                                                                         <  4:75196/103‑1 (MQ=255)
                                                                                  gATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCg                                                                                                         >  3:75196/1‑103 (MQ=255)
                                                                                   aTCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTAt                                                          <  4:1245165/149‑1 (MQ=255)
                                                                                           gATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTc                                                  <  4:827732/149‑1 (MQ=255)
                                                                                                      taCCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACa                                                                                                  <  2:153417/90‑1 (MQ=255)
                                                                                                      taCCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACa                                                                                                  >  1:153417/1‑90 (MQ=255)
                                                                                                      taCCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGaa                                       <  3:789458/149‑1 (MQ=255)
                                                                                                         ccGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGAttt                                                                                                               >  3:636063/1‑74 (MQ=255)
                                                                                                         ccGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGAttt                                                                                                               <  4:636063/74‑1 (MQ=255)
                                                                                                          cGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGa                                   <  2:36521/149‑1 (MQ=255)
                                                                                                              atatGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTcc                               <  1:49426/149‑1 (MQ=255)
                                                                                                                     cATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCta                        >  3:486563/1‑149 (MQ=255)
                                                                                                                      aTTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTaaa                                                                >  1:541202/1‑108 (MQ=255)
                                                                                                                      aTTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTaaa                                                                <  2:541202/108‑1 (MQ=255)
                                                                                                                          tgAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTaaaa                                                               >  2:1393339/1‑105 (MQ=255)
                                                                                                                          tgAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTaaaa                                                               <  1:1393339/105‑1 (MQ=255)
                                                                                                                           gAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCtatttatt                  >  3:586592/1‑149 (MQ=255)
                                                                                                                             aGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAg                                                                                        >  3:639762/1‑77 (MQ=255)
                                                                                                                             aGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAg                                                                                        <  4:639762/77‑1 (MQ=255)
                                                                                                                              gCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGaa                                                                                               <  4:1340660/69‑1 (MQ=255)
                                                                                                                              gCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGaa                                                                                               >  3:1340660/1‑69 (MQ=255)
                                                                                                                                aCTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCttat             <  4:486563/149‑1 (MQ=255)
                                                                                                                                  tAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGaa                                                                                               <  3:265516/65‑1 (MQ=255)
                                                                                                                                  tAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGaa                                                                                               >  4:265516/1‑65 (MQ=255)
                                                                                                                                   aaTGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGa                                   >  4:326894/1‑124 (MQ=255)
                                                                                                                                   aaTGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGa                                   <  3:326894/124‑1 (MQ=255)
                                                                                                                                       gCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCTTATTACTAtt      >  4:943995/1‑149 (MQ=255)
                                                                                                                                           tACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCTTATTACTATTTCGt  >  1:119771/1‑149 (MQ=255)
                                                                                                                                             cGTGGAGCAGGACACTCCTGCCCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTc                 <  4:1408485/132‑1 (MQ=255)
                                                                                                                                             cGTGGAGCAGGACACTCCTGCCCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTc                 >  3:1408485/1‑132 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
TAAGGATTCGGCCATTGTTCCTGTTTATTACTACGTGAATGCGCGTCTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCTTATTACTATTTCGT  >  NC_000913/1302681‑1302968

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: