Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,836,905 G→C 7.6% S278T (AGC→ACC) ‡ ynjB → putative ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,836,9050GC 7.6% 286.4 / 12.8 92S278T (AGC→ACC) ‡ynjBputative ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base G (40/44);  new base C (4/3);  total (44/48)
Fisher's exact test for biased strand distribution p-value = 7.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.38e-01

CAATATCTTGATGTGCTGCATCCGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCCTGCTCTCACCCGATGCGCAACTGC  >  NC_000913/1836757‑1837054
                                                                                                                                                    |                                                                                                                                                     
cAATATCTTGATGTGCTGCATCCGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGCCCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAgc                                                                                                                                                      >  1:950574/1‑150 (MQ=255)
    atCTTGATGTGCTGCATCCGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAgcggcg                                                                                                                                                  <  2:611961/150‑1 (MQ=255)
       ttGATGTGCTGCATCCGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGAtt                                                                                                                                               >  1:939997/1‑150 (MQ=255)
        tGATGTGCTGCATCCGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGAttt                                                                                                                                              >  1:1105253/1‑150 (MQ=255)
         gATGTGCTGCATCCGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCGCTGGCGAtttt                                                                                                                                             <  1:514158/150‑2 (MQ=255)
           tgtgCTGCATCCGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTAAACCCCGCGCATGCGCAGCAAAAAATCGCCGCTGGCGAtttttt                                                                                                                                           <  2:1184975/150‑4 (MQ=255)
            gtgCTGCATCCGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCGCTGGCGAttttttg                                                                                                                                          <  1:1325110/150‑5 (MQ=255)
             tgctgcATCCGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTg                                                                                                                                         <  2:985254/150‑1 (MQ=255)
             tgctgcATCCGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTg                                                                                                                                         >  2:560265/1‑150 (MQ=255)
               ctgcATCCGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCa                                                                                                                                       <  2:1283806/150‑1 (MQ=255)
               ctgcATCCGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCa                                                                                                                                       >  1:1576010/1‑150 (MQ=255)
               ctgcATCCGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCa                                                                                                                                       <  2:1189472/150‑1 (MQ=255)
                 gcATCCGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAg                                                                                                                                     <  1:670766/150‑1 (MQ=255)
                 gcATCCGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAg                                                                                                                                     >  1:565757/1‑150 (MQ=255)
                  cATCCGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACATTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGc                                                                                                                                    >  2:1533950/1‑150 (MQ=255)
                   aTCCGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCa                                                                                                                                   >  1:453710/1‑150 (MQ=255)
                     ccGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCCGCGGCGATTTGCCTGCAAGCAGt                                                                                                                                 >  1:482476/1‑150 (MQ=255)
                      cGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGtt                                                                                                                                >  1:1527872/1‑150 (MQ=255)
                          tttGTGGCGCGAAGTAAAGGATTTCCCGCCTTCACCCGCGCGGAGGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAg                                                                                                                            <  1:317194/150‑1 (MQ=255)
                            tgtgGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGtt                                                                                                                          <  2:242230/150‑1 (MQ=255)
                                  gcgAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCt                                                                                                                    >  1:1448590/1‑150 (MQ=255)
                                           aggaTTTCCCGCCTGCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGTCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAgg                                                                                                           <  1:1532281/150‑1 (MQ=255)
                                              aTTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGga                                                                                                        >  1:488010/1‑150 (MQ=255)
                                              aTTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGga                                                                                                        >  2:20300/1‑150 (MQ=255)
                                              aTTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGga                                                                                                        >  1:960496/1‑150 (MQ=255)
                                                ttCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGgatg                                                                                                      >  1:1720351/1‑150 (MQ=255)
                                                    cGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTg                                                                                                  >  2:750650/1‑150 (MQ=255)
                                                      ccTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTTCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGAACGCAGTTACAGTTTTGGCTTTCGCGCGGGGATGAATGGc                                                                                                >  2:313381/1‑150 (MQ=255)
                                                         tCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAAc                                                                                             >  1:1666769/1‑150 (MQ=255)
                                                         tCACCCACGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAAc                                                                                             <  1:1545287/150‑1 (MQ=255)
                                                          cACCCGCGCTGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACg                                                                                            <  2:960496/150‑1 (MQ=255)
                                                          cACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACg                                                                                            >  2:829289/1‑150 (MQ=255)
                                                           aCCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTACCCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGt                                                                                           >  2:91639/1‑150 (MQ=255)
                                                             ccGCGCGGATGGATGCTCTGCTGAAAGCCTGCTCATTGCGCCTGGCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGTCTTTCGCGAGGGGATGATTGGCAACGTGc                                                                                         <  1:84450/150‑1 (MQ=255)
                                                             ccGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGc                                                                                         <  1:1388777/150‑1 (MQ=255)
                                                             ccGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGc                                                                                         <  1:1678880/150‑1 (MQ=255)
                                                              cgcgcgGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGAtt                                                                                                   >  1:433490/1‑139 (MQ=255)
                                                              cgcgcgGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGAtt                                                                                                   <  2:433490/139‑1 (MQ=255)
                                                              cgcgcgGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCa                                                                                        >  2:938976/1‑150 (MQ=255)
                                                                cgcgGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCAtt                                                                                      <  2:857235/150‑1 (MQ=255)
                                                                cgcgGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCATGCAAGCAGTTACAGTTTTGGCTTACGCGAGGGGAAGATTGCCAACGTGCAtt                                                                                      >  2:144172/1‑150 (MQ=255)
                                                                 gcgGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCAttt                                                                                     <  2:1576010/150‑1 (MQ=255)
                                                                        gatgCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCAcc                                                                              <  1:560265/150‑1 (MQ=255)
                                                                        gatgCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCAcc                                                                              >  2:720143/1‑150 (MQ=255)
                                                                        gatgCTCTGCTGAAAGCCGGCACATTGCGCCGGTCGCTGACCTTTACCCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCAcc                                                                              <  2:950574/150‑1 (MQ=255)
                                                                          tgCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGAAGATTGGCAACGTGCATTTCGTCACCAt                                                                            <  1:288877/150‑1 (MQ=255)
                                                                          tgCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCAtt                                                                           >  1:1434425/1‑150 (MQ=255)
                                                                            ctctGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGATACAGTTTTGGCTTTCGCGAGGGGATGATTGGAAACGTGCATTTCGTCACCATTc                                                                          >  1:863948/1‑150 (MQ=255)
                                                                               tgctgAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTg                                                                       >  1:1192133/1‑150 (MQ=255)
                                                                                  tgAAAGCCGGCACATTGCGCCGGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCa                                                                    <  2:281865/150‑1 (MQ=255)
                                                                                     aaGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAAcg                                                                 <  2:453710/150‑1 (MQ=255)
                                                                                       gCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCTTGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGTTGATTGGCAACGTGCATTTCGGCTCCATTCCTGCCAAcgcg                                                               <  1:682008/150‑1 (MQ=255)
                                                                                       gCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAAcgcg                                                               >  1:1446552/1‑150 (MQ=255)
                                                                                        ccGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGa                                                              >  2:1075416/1‑150 (MQ=255)
                                                                                           gCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATg                                                           >  2:1218015/1‑150 (MQ=255)
                                                                                            cacaTTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGc                                                          <  1:938976/150‑1 (MQ=255)
                                                                                              caTTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAAACGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCa                                                        <  1:158390/150‑1 (MQ=255)
                                                                                                ttGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGt                                                      <  2:1527872/150‑1 (MQ=255)
                                                                                                 tGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGtg                                                     <  2:939997/150‑1 (MQ=255)
                                                                                                 tGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGtg                                                     <  2:1257124/150‑1 (MQ=255)
                                                                                                     ccTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGtgctgc                                                 >  1:1239717/1‑150 (MQ=255)
                                                                                                           gCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGaa                                           >  2:1188160/1‑150 (MQ=255)
                                                                                                           gCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGaa                                           <  1:1417740/150‑1 (MQ=255)
                                                                                                           gCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGAGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGaa                                           <  1:235650/150‑1 (MQ=255)
                                                                                                                  tttAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGtt                                    <  1:1582831/150‑1 (MQ=255)
                                                                                                                   ttAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTg                                   <  2:569412/150‑1 (MQ=255)
                                                                                                                       ccccGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGGAGTTACAGGTTTGGCTTTCGCGAGGGGATGATTGGCCACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGATGGTAGTTGcctt                               >  2:1367204/1‑148 (MQ=255)
                                                                                                                        cccGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAAt                              >  2:1692973/1‑150 (MQ=255)
                                                                                                                            cgcATGCGCAGCAAAACATCGCCAGCGGCGATTTGACTGCAAGCAGGTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCCACGGGCATTTCGTCACCATTCCTGCCTACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTcc                          <  1:316615/150‑1 (MQ=255)
                                                                                                                            cgcATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTcc                          <  2:1720351/150‑1 (MQ=255)
                                                                                                                             gcATGCGCAGCATAAAATCGCCAGCGGCGATTTGCCTTCAAGCAGTTACTGTTATGGCTTTCGCGAGGGTATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCTAATGCCAGTGCTGCGGCGAAGGTAGTTGACAATTTCCt                         <  1:219785/150‑1 (MQ=255)
                                                                                                                              cATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCCTg                        >  2:1542740/1‑150 (MQ=255)
                                                                                                                               aTGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCCTGc                       <  1:1218015/150‑1 (MQ=255)
                                                                                                                               aTGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCGTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCCTGc                       <  2:1280912/150‑1 (MQ=255)
                                                                                                                                      tcaAAAAATCGCCATCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCCTGCTCTCAcc                <  2:863948/149‑1 (MQ=255)
                                                                                                                                      gcaAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCCTGCTCTCAcc                <  2:1448590/150‑1 (MQ=255)
                                                                                                                                      gcaAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCCTGCTCTCAcc                <  1:169345/150‑1 (MQ=255)
                                                                                                                                      gcaAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCCTGCTCTCAcc                <  2:1666769/150‑1 (MQ=255)
                                                                                                                                      gcaAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCCTGCTCTCAcc                >  2:1012307/1‑150 (MQ=255)
                                                                                                                                      caggcaaaTCGCCGCTGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCCTGCTCTCAcc                >  2:680719/6‑150 (MQ=255)
                                                                                                                                      caggcaaaTCGCCGCTGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCCTGCTCTCAcc                >  2:806046/6‑150 (MQ=255)
                                                                                                                                       caAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCCTGCTCTCAccc               <  2:1192133/150‑1 (MQ=255)
                                                                                                                                        ggcaaaTCGCCGCTGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCCTGCTCTCACCCg              >  2:1691937/4‑150 (MQ=255)
                                                                                                                                        ggcaaaTCGCCGCTGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCCTGCTCTCACCCg              >  2:795453/4‑150 (MQ=255)
                                                                                                                                        aaaaaaTCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCCTGCTCTCACCCg              >  1:12509/1‑150 (MQ=255)
                                                                                                                                         aaaaaTCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCCTGCTCTCACCCGa             >  1:1660164/1‑150 (MQ=255)
                                                                                                                                          aaaaTCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCCTGCTCTCACCCGAt            >  1:740380/1‑150 (MQ=255)
                                                                                                                                            aaTCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCCTGCTCTCACCCGATgc          <  1:1075416/150‑1 (MQ=255)
                                                                                                                                              tCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCCTGCTCTCACCCGATgcgc        <  2:1105253/150‑1 (MQ=255)
                                                                                                                                               cGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGTGAAGGTAGTTGCCAATTTCCTGCTCTCACCCGATGCGCa       <  1:20300/150‑1 (MQ=255)
                                                                                                                                               cGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCCTGCTCTCACCCGATGCGCa       <  1:59299/150‑1 (MQ=255)
                                                                                                                                                  cAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCCTGCTCTCACCCGATGCGCAACt    >  2:1179709/1‑150 (MQ=255)
                                                                                                                                                  cAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCCTGCTCTCACCCGATGCGCAACt    >  2:355416/1‑150 (MQ=255)
                                                                                                                                                   aGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCCTGCTCTCACCCGATGCGCAACTg   <  1:829289/150‑1 (MQ=255)
                                                                                                                                                    gcggcgATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCCTGCTCTCACCCGATGCGCAACTGc  >  1:1352418/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CAATATCTTGATGTGCTGCATCCGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCCTGCTCTCACCCGATGCGCAACTGC  >  NC_000913/1836757‑1837054

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