Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,901,405 C→T L272L (CTG→CTA bglH ← carbohydrate‑specific outer membrane porin, cryptic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,901,4050CT96.0% 261.7 / 0.9 74L272L (CTG→CTAbglHcarbohydrate‑specific outer membrane porin, cryptic
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base T (34/37);  total (36/38)
Fisher's exact test for biased strand distribution p-value = 6.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01

GTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGATTTCAATTTTCGGCGCACCGTGTTTACCCACCCAGAAATCAGCCTCTGGC  >  NC_000913/3901261‑3901552
                                                                                                                                                |                                                                                                                                                   
gTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGcgcg                                                                                                                                               <  1:188936/151‑1 (MQ=255)
 tttGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAgcgc                                                                                                                                                >  1:603790/1‑149 (MQ=255)
   tGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGat                                                                                                                                             >  1:1068307/1‑150 (MQ=255)
       gTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCa                                                                                                                                         >  1:59739/1‑150 (MQ=255)
            aTAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAAtttt                                                                                                                                    >  1:133449/1‑150 (MQ=255)
              aaCGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTTAc                                                                                                                                  <  2:1205097/150‑1 (MQ=255)
                cGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCg                                                                                                                               <  1:578078/151‑1 (MQ=255)
                   ccACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGAc                                                                                                                            >  2:949918/1‑151 (MQ=255)
                     aCTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCg                                                                                                                          <  1:482436/151‑1 (MQ=255)
                      cTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCg                                                                                                                          <  2:313203/150‑1 (MQ=255)
                      cTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCg                                                                                                                          >  2:443509/1‑150 (MQ=255)
                          ccATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACttt                                                                                                                     <  1:1171485/151‑1 (MQ=255)
                              tAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTcc                                                                                                                   >  2:289738/1‑149 (MQ=255)
                              tAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTTACCCGGACCGACTTTcc                                                                                                                   >  2:1259009/1‑149 (MQ=255)
                                aGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGtt                                                                                                               <  1:595637/151‑1 (MQ=255)
                                        ttATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGa                                                                                                        <  2:94319/150‑1 (MQ=255)
                                          aTCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGAcc                                                                                                      <  1:177303/150‑1 (MQ=255)
                                          aTCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCt                                                                                                     >  2:616615/1‑151 (MQ=255)
                                           tCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCt                                                                                                     >  1:449529/1‑150 (MQ=255)
                                               aTAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACAcc                                                                                                <  2:1037706/151‑1 (MQ=255)
                                                 aaCGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCg                                                                                              >  2:316274/1‑151 (MQ=255)
                                                  aCGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGc                                                                                             <  2:95671/151‑1 (MQ=255)
                                                   cGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCt                                                                                            <  1:697963/151‑1 (MQ=255)
                                                   cGGGATATCTTTATAGCGTAAATCAATGGTATTGGAATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGc                                                                                             >  2:1120564/1‑150 (MQ=255)
                                                      gATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCg                                                                                         <  2:760975/151‑1 (MQ=255)
                                                      gATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCg                                                                                         >  2:80931/1‑151 (MQ=255)
                                                            tttATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATc                                                                                    >  2:111609/1‑150 (MQ=255)
                                                            tttATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATc                                                                                    >  2:807614/1‑150 (MQ=255)
                                                               ataGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTa                                                                                <  2:281534/151‑1 (MQ=255)
                                                               ataGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTa                                                                                <  2:323954/151‑1 (MQ=255)
                                                                   cGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCt                                                                            >  1:978978/1‑151 (MQ=255)
                                                                    gTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTg                                                                           <  2:187488/151‑1 (MQ=255)
                                                                      aaaTCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCg                                                                         >  1:918219/1‑151 (MQ=255)
                                                                       aaTCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGt                                                                        <  1:834848/151‑1 (MQ=255)
                                                                       aaTCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGt                                                                        <  2:1133800/151‑1 (MQ=255)
                                                                       aaTCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGt                                                                        >  2:977656/1‑151 (MQ=255)
                                                                        aTCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGtt                                                                       <  2:52559/151‑1 (MQ=255)
                                                                           aaTGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTc                                                                    <  1:388938/151‑1 (MQ=255)
                                                                              ggtattggtattAATCTTCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCa                                                                  <  1:1101379/150‑1 (MQ=255)
                                                                                 attggtattAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTc                                                               >  2:307150/1‑150 (MQ=255)
                                                                                   tggtattAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCaa                                                             <  1:443498/150‑1 (MQ=255)
                                                                                      tattAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCAt                                                         >  1:532421/1‑151 (MQ=255)
                                                                                       attAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATc                                                        <  2:1088247/151‑1 (MQ=255)
                                                                                          aaTCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTgg                                                     >  1:648925/1‑151 (MQ=255)
                                                                                            tCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCCTCGTCATCTTCGCGAACTAGCGCGACATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGTGGCATCAGTCCGCTGCGTTTTCCAGTCAAGCATCTGGAt                                                   >  2:948264/1‑151 (MQ=255)
                                                                                            tCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGAt                                                   >  2:156779/1‑151 (MQ=255)
                                                                                             ctgctgTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTAGCGAACTAGCGCGATATCAATTTTACCCGAACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGTTGCGTTTTCCAGTCAAGCTTCTGGAtt                                                  >  1:79663/1‑151 (MQ=255)
                                                                                                  gtttgtttTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGATTTCaa                                              <  2:555759/150‑1 (MQ=255)
                                                                                                   tttgtttTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGATTTCAAtt                                            >  1:213184/1‑151 (MQ=255)
                                                                                                     tgtttTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGATTTCAAtttt                                          >  1:282410/1‑151 (MQ=255)
                                                                                                     tgtttTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGATTTCAAtttt                                          >  2:925535/1‑151 (MQ=255)
                                                                                                      gtttTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGATTTCAATTTTc                                         >  2:260835/1‑151 (MQ=255)
                                                                                                         ttGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGATTTCAATTTTCgg                                       <  2:449540/150‑1 (MQ=255)
                                                                                                          tGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGATTTCAATTTTCGgc                                      <  2:533446/150‑1 (MQ=255)
                                                                                                            cAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGATTTCAATTTTCGgcgc                                    >  1:254276/1‑150 (MQ=255)
                                                                                                            cAGGCTGCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGATTTCAATTTTCGgcgc                                    >  2:265619/1‑150 (MQ=255)
                                                                                                                  gCGATCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGATTTCAATTTTCGGCGCACCgtgt                             >  1:658835/1‑151 (MQ=255)
                                                                                                                     aTCGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGATTTCAATTTTCGGCGCACCGTGTTTa                          <  2:1132260/151‑1 (MQ=255)
                                                                                                                       cGTAATCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGATTTCAATTTTCGGCGCACCGTGTTTacc                        >  2:1213183/1‑151 (MQ=255)
                                                                                                                          aaTCATCGACGTCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGATTTCAATTTTCGGCGCACCGTGTTTacccac                     <  2:258742/151‑1 (MQ=255)
                                                                                                                          aaTCATCGACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGATTTCAATTTTCGGCGCACCGTGTTTacccac                     <  2:624609/151‑1 (MQ=255)
                                                                                                                                 gACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGATTTCAATTTTCGGCGCACCGTGTTTACCCACCCAGaaa              <  2:1065797/151‑1 (MQ=255)
                                                                                                                                 gACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGATTTCAATTTTCGGCGCACCGTGTTTACCCACCCAGaaa              <  1:1019419/151‑1 (MQ=255)
                                                                                                                                 gACATCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGATTTCAATTTTCGGCGCACCGTGTTTACCCACCCAGaa               >  1:605081/1‑150 (MQ=255)
                                                                                                                                     tCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGATTTCAATTTTCGGCGCACCGTGTTTACCCACCCAGAAATc            >  1:470163/1‑149 (MQ=255)
                                                                                                                                     tCTTCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGATTTCAATTTTCGGCGCACCGTGTTTACCCACCCAGAAATc            >  1:915504/1‑149 (MQ=255)
                                                                                                                                        tCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGATTTCAATTTTCGGCGCACCGTGTTTACCCACCCAGAAATCAGcc        <  2:735409/150‑1 (MQ=255)
                                                                                                                                        tCGCGAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGATTTCAATTTTCGGCGCACCGTGTTTACCCACCCAGAAATCAGcc        <  2:893510/150‑1 (MQ=255)
                                                                                                                                         cgcgAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGATTTCAATTTTCGGCGCACCGTGTTTACCCACCCAGAAATCAGCctc      <  2:258905/151‑1 (MQ=255)
                                                                                                                                          gcgAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGATTTCAATTTTCGGCGCACCGTGTTTACCCACCCAGAAATCAGCctct     >  2:1213033/1‑151 (MQ=255)
                                                                                                                                          gcgAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGATTTCAATTTTCGGCGCACCGTGTTTACCCACCCAGAAATCAGCctct     <  1:157203/151‑1 (MQ=255)
                                                                                                                                           cgAACTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGATTTCAATTTTCGGCGCACCGTGTTTACCCACCCAGAAATCAGCCTCTg    <  2:203612/151‑1 (MQ=255)
                                                                                                                                             aaCTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGATTTCAATTTTCGGCGCACCGTGTTTACCCACCCAGAAATCAGCCTCTgg   <  1:705735/150‑1 (MQ=255)
                                                                                                                                             aaCTAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGATTTCAATTTTCGGCGCACCGTGTTTACCCACCCAGAAATCAGCCTCTGgc  <  1:844736/151‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
GTTTGCCGTGACATAACGACCACTTACCATTAAGGTGGCTTTATCCCATAACGGGATATCTTTATAGCGTAAATCAATGGTATTGGTATTAATCTGCTGTTTGTTTTGCAGGCTGCGATCGTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGATTTCAATTTTCGGCGCACCGTGTTTACCCACCCAGAAATCAGCCTCTGGC  >  NC_000913/3901261‑3901552

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