Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,839,568 G→C 5.4% N142K (AAC→AAG)  nlpA ← lipoprotein‑28

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,839,5680GC 5.4% 317.5 / 7.1 92N142K (AAC→AAG) nlpAlipoprotein‑28
Reads supporting (aligned to +/- strand):  ref base G (47/40);  new base C (3/2);  total (50/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.35e-01

AGTACGCGCGGGAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTACCAGTTTATAACCATGCGCCTGATTATCTTGTTCAAG  >  NC_000913/3839419‑3839717
                                                                                                                                                     |                                                                                                                                                     
aGTACGCGCGGGAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAgg                                                                                                                                                       <  2:228736/150‑1 (MQ=255)
 gTACGCGCGGGAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGt                                                                                                                                                      >  2:158323/1‑150 (MQ=255)
   aCGCGCGGGAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTg                                                                                                                                                    <  1:1021354/150‑1 (MQ=255)
      cgcgGGAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTggtc                                                                                                                                                 >  1:1637558/1‑150 (MQ=255)
          ggAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAACCTTggtcgcgc                                                                                                                                             <  1:321918/150‑4 (MQ=255)
           gAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTggtcgggtc                                                                                                                                            >  2:980686/1‑150 (MQ=255)
              cTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCCTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGtt                                                                                                                                         >  2:1256445/1‑150 (MQ=255)
               tGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTg                                                                                                                                        >  1:819194/1‑150 (MQ=255)
               tGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTg                                                                                                                                        >  1:423010/1‑150 (MQ=255)
               tGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTg                                                                                                                                        <  1:387618/150‑1 (MQ=255)
               tGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTg                                                                                                                                        >  1:181405/1‑150 (MQ=255)
                gcgcCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTgg                                                                                                                                       <  1:176435/150‑1 (MQ=255)
                gcgcCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTgg                                                                                                                                       <  1:159861/150‑1 (MQ=255)
                 cgcCCCCTCGAGTTCCATAATCTGCAAATGAGGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTggg                                                                                                                                      >  2:989339/1‑150 (MQ=255)
                       cTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGc                                                                                                                                >  1:444203/1‑150 (MQ=255)
                        tCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGcc                                                                                                                               <  2:1460900/150‑1 (MQ=255)
                         cGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCa                                                                                                                              >  2:475493/1‑150 (MQ=255)
                         cGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCa                                                                                                                              >  2:237672/1‑150 (MQ=255)
                          gAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCAc                                                                                                                             <  2:106808/150‑1 (MQ=255)
                            gTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTg                                                                                                                           >  1:804208/1‑150 (MQ=255)
                               ccATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTcg                                                                                                                        >  1:1617276/1‑150 (MQ=255)
                               ccATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTcg                                                                                                                        >  1:238325/1‑150 (MQ=255)
                                cATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTcgc                                                                                                                       <  2:181405/150‑1 (MQ=255)
                                cATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTcgc                                                                                                                       <  1:1429628/150‑1 (MQ=255)
                                 aTAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGGAACAGTAAAAGTGCGCGACCAAGGTTggtcggg                                                                                                                                              >  1:908876/1‑126 (MQ=255)
                                 aTAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTggtcggg                                                                                                                                              <  2:908876/126‑1 (MQ=255)
                                 aTAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTcgcg                                                                                                                      <  1:158323/150‑1 (MQ=255)
                                  tAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCgcgc                                                                                                                     >  2:958003/1‑150 (MQ=255)
                                         cAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCttt                                                                                                              >  1:1476921/1‑150 (MQ=255)
                                         cAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCttt                                                                                                              >  2:579541/1‑150 (MQ=255)
                                             tGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTAtt                                                                                                          <  2:819194/150‑1 (MQ=255)
                                             tGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTAGCGCCTTCTTTTAtt                                                                                                          <  1:1630779/150‑1 (MQ=255)
                                               aCGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTg                                                                                                        <  1:813859/150‑1 (MQ=255)
                                                 gAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTgcg                                                                                                      <  1:1571309/150‑1 (MQ=255)
                                                          cGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGttt                                                                                             <  1:231050/150‑1 (MQ=255)
                                                           ggTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCct                                                                                                                   <  1:220434/127‑1 (MQ=255)
                                                           ggTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCct                                                                                                                   >  2:220434/1‑127 (MQ=255)
                                                           ggTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGGCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGtttt                                                                                            <  1:1576078/150‑1 (MQ=255)
                                                               atatCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGAtt                                                                                        >  2:1292423/1‑150 (MQ=255)
                                                                tatCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGAttt                                                                                       >  1:1088257/1‑150 (MQ=255)
                                                                tatCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGAttt                                                                                       >  1:1082535/1‑150 (MQ=255)
                                                                   ccAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTGGAtttttt                                                                                    >  2:760565/1‑150 (MQ=255)
                                                                     aGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTg                                                                                  <  1:67268/150‑1 (MQ=255)
                                                                       cgcgGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGa                                                                                <  1:1014085/150‑1 (MQ=255)
                                                                              ggtaATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAAccc                                                                         <  2:1363947/150‑1 (MQ=255)
                                                                                taataaCCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGc                                                                       >  1:1250981/1‑150 (MQ=255)
                                                                                 aataaCCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTTTTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGcc                                                                      >  2:1305484/1‑150 (MQ=255)
                                                                                       ccTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAggg                                                                <  1:1592454/150‑1 (MQ=255)
                                                                                       ccTTTTCCCTCATTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCTCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAggg                                                                <  2:756067/150‑1 (MQ=255)
                                                                                        cTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGa                                                               >  1:423266/1‑150 (MQ=255)
                                                                                           ttCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGa                                                            >  1:702958/1‑150 (MQ=255)
                                                                                            tCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGAc                                                           >  1:706012/1‑150 (MQ=255)
                                                                                             cccTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCAt                                                                    <  1:259550/140‑1 (MQ=255)
                                                                                             cccTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCAt                                                                    >  2:259550/1‑140 (MQ=255)
                                                                                             cccTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACa                                                          >  2:1344236/1‑150 (MQ=255)
                                                                                               ctctTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACaaa                                                        >  1:494822/1‑150 (MQ=255)
                                                                                                 ctTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACg                                                      <  1:763051/150‑1 (MQ=255)
                                                                                                 ctTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACg                                                      <  1:1248858/150‑1 (MQ=255)
                                                                                                     cAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTAtt                                                  <  2:201931/150‑1 (MQ=255)
                                                                                                       gagTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTc                                                >  1:1340046/1‑150 (MQ=255)
                                                                                                         gTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTccc                                              >  1:224620/1‑150 (MQ=255)
                                                                                                           aaTTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGGCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCAc                                            <  2:106792/150‑1 (MQ=255)
                                                                                                            aTTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCAcc                                           <  2:238325/150‑1 (MQ=255)
                                                                                                             ttAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGGCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCg                                          <  2:496629/150‑1 (MQ=255)
                                                                                                               aGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCt                                        >  2:50008/1‑150 (MQ=255)
                                                                                                                gCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTa                                       <  2:1402146/150‑1 (MQ=255)
                                                                                                                 cTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTAc                                      <  2:804208/150‑1 (MQ=255)
                                                                                                                      gtttttGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTAGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTACCAGtt                                 <  1:970363/149‑1 (MQ=255)
                                                                                                                      ctTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTACCAGtt                                 >  1:595075/1‑150 (MQ=255)
                                                                                                                      ctTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTACCAGtt                                 >  2:234066/1‑150 (MQ=255)
                                                                                                                        ttttGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTACCAGTTta                               <  1:1344236/150‑1 (MQ=255)
                                                                                                                            gCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTACCAGTTTATAAc                           >  1:764295/1‑150 (MQ=255)
                                                                                                                             cAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTACCAGTTTATAAcc                          <  2:1617276/150‑1 (MQ=255)
                                                                                                                             cAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTACCAGTTTATAAcc                          >  2:256730/1‑150 (MQ=255)
                                                                                                                             cAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTACCAGTTTATAAcc                          <  1:475493/150‑1 (MQ=255)
                                                                                                                             cAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTACCAGTTTATAAcc                          >  1:1509082/1‑150 (MQ=255)
                                                                                                                              aaCAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTACCAGTTTATAACCa                         <  1:1256445/150‑1 (MQ=255)
                                                                                                                                cAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTACCAGTTTATAACCATg                       <  1:237672/150‑1 (MQ=255)
                                                                                                                                   tAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTACCAGtt                                 <  1:794879/137‑1 (MQ=255)
                                                                                                                                   tAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTACCAGtt                                 >  2:794879/1‑137 (MQ=255)
                                                                                                                                   tAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTACCAGTTTATAACCATgcgc                    >  1:281703/1‑150 (MQ=255)
                                                                                                                                       acgacccgACCAACCTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTAcc                                     <  2:387214/123‑1 (MQ=255)
                                                                                                                                       acgacccgACCAACCTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTAcc                                     >  1:387214/7‑129 (MQ=255)
                                                                                                                                       aGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTACCAGTTTATAACCATGCGCCTGa                >  2:1205651/1‑150 (MQ=255)
                                                                                                                                         tGCGCGCCCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATACCCCGCCATAGGGAAGACAAACGGATTTCCCACCGCTACCAGTTTATAACCATGCGCCTGAtt              >  1:259695/1‑150 (MQ=255)
                                                                                                                                          gcgcGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTACCAGTTTATAACCATGCGCCTGATTa             >  1:752366/1‑150 (MQ=255)
                                                                                                                                          gcgcGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTACCAGTTTATAACCATGCGCCTGATTa             <  1:234066/150‑1 (MQ=255)
                                                                                                                                          acccgACCAACCTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTACCAGTTTATAACCATGCGCCTGATTa             >  1:1631024/4‑150 (MQ=255)
                                                                                                                                          acccgACCAACCTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTACCAGTTTATAACCATGCGCCTGATTa             >  1:513744/4‑150 (MQ=255)
                                                                                                                                               aCCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTACCAGTTTATAACCATGCGCCTGATTATCTTg        >  2:469677/1‑150 (MQ=255)
                                                                                                                                                ccAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTACCAGTTTATAACCATGCGCCTGATTATCTTGt       <  2:1188566/150‑1 (MQ=255)
                                                                                                                                                 cAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTACCAGTTTATAACCATGCGCCTGATTATCTTGtt      >  1:1259938/1‑150 (MQ=255)
                                                                                                                                                 cAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTACCAGTTTATAACCATGCGCCTGATTATCTTGtt      >  1:793055/1‑150 (MQ=255)
                                                                                                                                                     gTTGGTCGGGTCGTTGGGGATCCCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTACCAGTTTATAACGATGCGCCTGATTATCTTGTTCAAg  <  2:804152/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGTACGCGCGGGAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTTCAGAGTAATTAGCTTCTCTTTTTGCAACAGTAAAAGTGCGCGACCAAGGTTGGTCGGGTCGTTGGGGATCGCCACTGTCGCGCCTTCTTTTATTTGCGCCACGGTTTTGATTTTTTTGGAATAACCCGCCATAGGGAAGACAAACGTATTTCCCACCGCTACCAGTTTATAACCATGCGCCTGATTATCTTGTTCAAG  >  NC_000913/3839419‑3839717

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑