Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 721,955 A→G G820G (GGT→GGC kdpD ← sensory histidine kinase KdpD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913721,9550AG86.8% 188.9 / 11.5 68G820G (GGT→GGCkdpDsensory histidine kinase KdpD
Reads supporting (aligned to +/- strand):  ref base A (4/5);  new base G (27/32);  total (31/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGCGGTAATAGTGCCCCCGTGTACATCCACTATCGCCCGACAAATTGCCAGTCCAAGCCCTACCCCCGGTACTGCCGACTCTTTATTCCCGCGAGCAAACTTATCAAATATCGTCTGCTCCTGGCCTGGCGGAAGACCGGGGCCGTTATCCCAGACATCCAGTTGTAGATTTTCGCCCTCAACGTGGGCATCGATACCAATTTCGGCCTGCGCACCCGCATATTTCACCGCGTTCTCCAGCAGATTAATCAGCACCCGTTCAAAGAGTGG  >  NC_000913/721819‑722089
                                                                                                                                        |                                                                                                                                      
aCGCGGTAATAGTGCCCCCGTGTACATCCACTATCGCCCGACAAATTGCCAGTCCAAGCCCTACCCCCGGTACTGCCGACTCTTTATTCCCGCGAGCAAACTTATCAAATATCGTCTGCTCCTGGCCTGGCGGAAGGCCg                                                                                                                                     <  1:612346/140‑1 (MQ=255)
      tAATAGTGCCCCCGTGTACATCCACTATCGCCCGACAAATTGCCAGTCCAAGCCCTACCCCCGGTACTGCCGACTCTTTATTCCCGCGAGCAAACTTATCAAATATCGTCTGCTCCTGGCCTGGCGGAAGGCCGGGGCCg                                                                                                                               <  1:728938/140‑1 (MQ=255)
        aTAGTGCCCCCGTGTACATCCACTATCGCCCGACAAATTGCCAGTCCAAGCCCTACCCCCGGTACTGCCGACTCTTTATTCCCGCGAGCAAACTTATCAAATATCGTCTGCTCCTGGCCTGGCGGAAGGCCGGGGCCGtt                                                                                                                             <  1:800594/140‑1 (MQ=255)
            tGCCCCCGTGTACATCCACTATCGCCCGACAAATTGCCAGTCCAAGCCCTACCCCCGGTACTGCCGACTCTTTATTCCCGCGAGCAAACTTATCAAATATCGTCTGCTCCTGGCCTGGCGGAAGGCCGGGGCCGTTATcc                                                                                                                         <  1:1167129/140‑1 (MQ=255)
              cccccGTGTACATCCACTATCGCCCGACAAATTGCCAGTCCAAGCCCTACCCCCGGTACTGCCGACTCTTTATTCCCGCGAGCAAACTTATCAAATATCGTCTGCTCCTGGCCTGGCGGAAGGCCGGGGCCGTTATCCCa                                                                                                                       >  1:96163/1‑140 (MQ=255)
              cccccGTGTACATCCACTATCGCCCGACAAATTGCCAGTCCAAGCCCTACCCCCGGTACTGCCGACTCTTTATTCCCGCGAGCAAACTTATCAAATATCGTCTGCTCCTGGCCTGGCGGAAGGCCGGGGCCGTTATCCCa                                                                                                                       >  2:673714/1‑140 (MQ=255)
              cccccGTGTACATCCACTATCGCCCGACAAATTGCCAGTCCAAGCCCTACCCCCGGTACTGCCGACTCTTTATTCCCGCGAGCAAACTTATCAAATATCGTCTGCTCCTGGCCTGGCGGAAGGCCGGGGCCGTTATCCCa                                                                                                                       >  2:954604/1‑140 (MQ=255)
               ccccGTGTACATCCACTATCGCCCGACAAATTGCCAGTCCAAGCCCTACCCCCGGTACTGCCGACTCTTTATTCCCGCGAGCAAACTTATCAAATATCGTCTGCTCCTGGCCTGGCGGAAGGCCGGGGCCGTTATCCCAg                                                                                                                      <  2:834307/140‑1 (MQ=255)
                cccGTGTACATCCACTATCGCCCGACAAATTGCCAGTCCAAGCCCTACCCCCGGTACTGCCGACTCTTTATTCCCGCGAGCAAACTTATCAAATATCGTCTGCTCCTGGCCTGGCGGAAGGCCGGGGCCGTTATCCCa                                                                                                                       <  1:68452/138‑1 (MQ=255)
                cccGTGTACATCCACTATCGCCCGACAAATTGCCAGTCCAAGCCCTACCCCCGGTACTGCCGACTCTTTATTCCCGCGAGCAAACTTATCAAATATCGTCTGCTCCTGGCCTGGCGGAAGGCCGGGGCCGTTATCCCa                                                                                                                       >  2:68452/1‑138 (MQ=255)
                cccGTGTACATCCACTATCGCCCGACAAATTGCCAGTCCAAGCCCTACCCCCGGTACTGCCGACTCTTTATTCCCGCGAGCAAACTTATCAAATATCGTCTGCTCCTGGCCTGGCGGAAGGCCGGGGCCGTTATCCCAGa                                                                                                                     >  1:598065/1‑140 (MQ=255)
                        cATCCACTATCGCCCGACAAATTGCCAGTCCAAGCCCTACCCCCGGTACTGCCGACTCTTTATTCCCGCGAGCAAACTTATCAAATATCGTCTGCTCCTGGCCTGGCGGAAGGCCGGGGCCGTTATCCCAGACATCCAGt                                                                                                             <  2:966311/140‑1 (MQ=255)
                             aCTATCGCCCGACAAATTGCCAGTCCAAGCCCTACCCCCGGTACTGCCGACTCTTTATTCCCGCGAGCAAACTTATCAAATATCGTCTGCTCCTGGCCTGGCGGAAGGCCGGGGCCGTTATCCCAGACATCCAGTTGTAg                                                                                                        >  2:69589/1‑140 (MQ=255)
                              cTATCGCCCGACAAATTGCCAGTCCAAGCCCTACCCCCGGTACTGCCGACTCTTTATTCCCGCGAGCAAACTTATCAAATATCGTCTGCTCCTGGCCTGGCGGAAGGCCGGGGCCGTTATCCCAGACATCCAGTTGTAGa                                                                                                       >  2:504999/1‑140 (MQ=255)
                              cTATCGCCCGACAAATTGCCAGTCCAAGCCCTACCCCCGGTACTGCCGACTCTTTATTCCCGCGAGCAAACTTATCAAATATCGTCTGCTCCTGGCCTGGCGGAAGGCCGGGGCCGTTATCCCAGACATCCAGTTGTAGa                                                                                                       >  2:1150209/1‑140 (MQ=255)
                                 tCGCCCGACAAATTGCCAGTCCAAGCCCTACCCCCGGTACTGCCGACTCTTTATTCCCGCGAGCAAACTTATCAAATATCGTCTGCTCCTGGCCTGGCGGAAGGCCGGGGCCGTTATCCCAGACATCCAGTTGTAGAttt                                                                                                    >  1:372516/1‑140 (MQ=255)
                                   gCCCGACAAATTGCCAGTCCAAGCCCTACCCCCGGTACTGCCGACTCTTTATTCCCGCGAGCAAACTTATCAAATATCGTCTGCTCCTGGCCTGGCGGAAGGCCGGGGCCGTTATCCCAGACATCCAGTTGTAGATTTTc                                                                                                  <  1:69589/140‑1 (MQ=255)
                                   gCCCGACAAATTGCCAGTCCAAGCCCTACCCCCGGTACTGCCGACTCTTTATTCCCGCGAGCAAACTTATCAAATATCGTCTGCTCCTGGCCTGGCGGAAGGCCGGGGCCGTTATCCCAGACATCCAGTTGTAGATTTTc                                                                                                  >  1:426560/1‑140 (MQ=255)
                                   gCCCGACAAATTGCCAGTCCAAGCCCTACCCCCGGTACTGCCGACTCTTTATTCCCGCGAGCAAACTTATCAAATATCGTCTGCTCCTGGCCTGGCGGAAGACCGGGGCCGTTATCCCAGACATCCAGTTGTAGATTTTc                                                                                                  >  1:230272/1‑140 (MQ=255)
                                    cccGACAAATTGCCAGTCCAAGCCCTACCCCCGGTACTGCCGACTCTTTATTCCCGCGAGCAAACTTATCAAATATCGTCTGCTCCTGGCCTGGCGGAAGGCCGGGGCCGTTATCCCAGACATCCAGTTGTAGATTTTCg                                                                                                 >  2:69680/1‑140 (MQ=255)
                                         cAAATTGCCAGTCCAAGCCCTACCCCCGGTACTGCCGACTCTTTATTCCCGCGAGCAAACTTATCAAATATCGTCTGCTCCTGGCCTGGCGGAAGGCCGGGGCCGTTATCCCAGACATCCAGTTGTAGATTTTCGCCCTc                                                                                            >  1:697195/1‑140 (MQ=255)
                                             ttGCCAGTCCAAGCCCTACCCCCGGTACTGCCGACTCTTTATTCCCGCGAGCAAACTTATCAAATATCGTCTGCTCCTGGCCTGGCGGAAGGCCGGGGCCGTTATCCCAGACATCCAGTTGTAGATTTTCGCCCTCAACg                                                                                        <  1:673714/140‑1 (MQ=255)
                                             ttGCCAGTCCAAGCCCTACCCCCGGTACTGCCGACTCTTTATTCCCGCGAGCAAACTTATCAAATATCGTCTGCTCCTGGCCTGGCGGAAGGCCGGGGCCGTTATCCCAGACATCCAGTTGTAGATTTTCGCCCTCAACg                                                                                        <  1:504999/140‑1 (MQ=255)
                                             ttGCCAGTCCAAGCCCTACCCCCGGTACTGCCGACTCTTTATTCCCGCGAGCAAACTTATCAAATATCGTCTGCTCCTGGCCTGGCGGAAGGCCGGGGCCGTTATCCCAGACATCCAGTTGTAGATTTTCGCCCTCAACg                                                                                        <  1:1150209/140‑1 (MQ=255)
                                                     ccAAGCCCTACCCCCGGTACTGCCGACTCTTTATTCCCGCGAGCAAACTTATCAAATATCGTCTGCTCCTGGCCTGGCGGAAGGCCGGGGCCGTTATCCCAGACATCCAGTTGTAGATTTTCGCCCTCAACGTGGGCATc                                                                                <  1:69680/140‑1 (MQ=255)
                                                         gCCCTACCCCCGGTACTGCCGACTCTTTATTCCCGCGAGCAAACTTATCAAATATCGTCTGCTCCTGGCCTGGCGGAAGGCCGGGGCCGTTATCCCAGACATCCAGTTGTAGATTTTCGCCCTc                                                                                            >  2:1171912/1‑124 (MQ=255)
                                                         gCCCTACCCCCGGTACTGCCGACTCTTTATTCCCGCGAGCAAACTTATCAAATATCGTCTGCTCCTGGCCTGGCGGAAGGCCGGGGCCGTTATCCCAGACATCCAGTTGTAGATTTTCGCCCTc                                                                                            <  1:1171912/124‑1 (MQ=255)
                                                          cccTACCCCCGGTACTGCCGACTCTTTATTCCCGCGAGCAAACTTATCAAATATCGTCTGCTCCTGGCCTGGCGGAAGGCCGGGGCCGTTATCCCAGACATCCAGTTGTAGATTTTCGCCCTCAACGTGGGCATCGATAc                                                                           <  2:426560/140‑1 (MQ=255)
                                                          cccTACCCCCGGTACTGCCGACTCTTTATTCCCGCGAGCAAACTTATCAAATATCGTCTGCTCCTGGCCTGGCGGAAGGCCGGGGCCGTTATCCCAGACATCCAGTTGTAGATTTTCGCCCTCAACGTGGGCATCGATAc                                                                           <  2:483154/140‑1 (MQ=255)
                                                              aCCCCCGGTACTGCCGACTCTTTATTCCCGCGAGCAAACTTATCAAATATCGTCTGCTCCTGGCCTGGCGGAAGGCCGGGGCCGTTATCCCAGACATCCAGTTGTAGATTTTCGCCCTCAACGTGGGCATCGATACCAAt                                                                       >  1:811108/1‑140 (MQ=255)
                                                               cccccGGTACTGCCGACTCTTTATTCCCGCGAGCAAACTTATCAAATATCGTCTGCTCCTGGCCTGGCGGAAGGCCGGGGCCGTTATCCCAGACATCCAGTTGTAGATTTTCGCCCTCAACGTGGGCATCGATACCAAtt                                                                      >  2:1035091/1‑140 (MQ=255)
                                                                ccccGGTACTGCCGACTCTTTATTCCCGCGAGCAAACTTATCAAATATCGTCTGCTCCTGGCCTGGCGGAAGGCCGGGGCCGTTATCCCAGACATCCAGTTGTAGATTTTCGCCCTCAACGTGGGCATCGATACCAAttt                                                                     <  2:96163/140‑1 (MQ=255)
                                                                     gTACTGCCGACTCTTTATTCCCGCGAGCAAACTTATCAAATATCGTCTGCTCCTGGCCTGGCGGAAGACCGGGGCCGTTATCCCAGACATCCAGTTGTAGATTTTCGCCCTCAACGTGGGCATCGATACCAATTTCGGcc                                                                >  2:397168/1‑140 (MQ=255)
                                                                         tGCCGACTCTTTATTCCCGCGAGCAAACTTATCAAATATCGTCTGCTCCTGGCCTGGCGGAAGGCCGGGGCCGTTATCCCAGACATCCAGTTGTAGATTTTCGCCCTCAACGTGGGCATCGATACCAATTTCGGCCTgcg                                                            <  1:1035091/140‑1 (MQ=255)
                                                                            cGACTCTTTATTCCCGCGAGCAAACTTATCAAATATCGTCTGCTCCTGGCCTGGCGGAAGGCCGGGGCCGTTATCCCAGACATCCAGTTGTAGATTTTCGCCCTCAACGTGGGCATCGATACCAATTTCGGCCTGCGCAc                                                         <  2:372516/140‑1 (MQ=255)
                                                                              aCTCTTTATTCCCGCGAGCAAACTTATCAAATATCGTCTGCTCCTGGCCTGGCGGAAGGCCGGGGCCGTTATCCCAGACATCCAGTTGTAGATTTTCGCCCTCAACGTGGGCATCGATACCAATTTCGGCCTGCGCAccc                                                       <  2:625841/140‑1 (MQ=255)
                                                                                tctTTATTCCCGCGAGCAAACTTATCAAATATCGTCTGCTCCTGGCCTGGCGGAAGGCCGGGGCCGTTATCCCAGACATCCAGTTGTAGATTTTCGCCCTCAATGTGGGCATCGATACCAATTTCGGCCTGCGCACCCGc                                                     >  2:1085363/1‑140 (MQ=255)
                                                                                           gcgAGCAAACTTATCAAATATCGTCTGCTCCTGGCCTGGCGGAAGGCCGGGGCCGTTATCCCAGACATCCAGTTGTAGATTTTCGCCCTCAACGTGGGCATCGATACCAATTTCGGCCTGCGCACCCGCATATTTCACcg                                          <  2:358168/140‑1 (MQ=255)
                                                                                            cgAGCAAACTTATCAAATATCGTCTGCTCCTGGCCTGGCGGAAGGCCGGGGCCGTTATCCCAGACATCCAGTTGTAGATTTTCGCCCTCAACGTGGGCATCGATACCAATTTCGGCCTGCGCACCCGCATATTTCACcgc                                         >  2:842252/1‑140 (MQ=255)
                                                                                            cgAGCAAACTTATCAAATATCGTCTGCTCCTGGCCTGGCGGAAGGCCGGGGCCGTTATCCCAGACATCCAGTTGTAGATTTTCGCCCTCAACGTGGGCATCGATACCAATTTCGGCCTGCGCACCCGCATATTTCACcgc                                         >  2:48121/1‑140 (MQ=255)
                                                                                            cgAGCAAACTTATCAAATATCGTCTGCTCCTGGCCTGGCGGAAGGCCGGGGCCGTTATCCCAGACATCCAGTTGTAGATTTTCGCCCTCAACGTGGGCATCGATACCAATTTCGGCCTGCGCACCCGCATATTTCACcgc                                         >  1:295489/1‑140 (MQ=255)
                                                                                                cAAACTTATCAAATATCGTCTGCTCCTGGCCTGGCGGAAGACCGGGGCCGTTATCCCAGACATCCAGTTGTAGATTTTCGCCCTCAACGTGGGCATCGATACCAATTTCGGCCTGCGCACCCGCATATTTCACCGCGTtc                                     >  1:852591/1‑140 (MQ=255)
                                                                                                   aCTTATCAAATATCGTCTGCTCCTGGCCTGGCGGAAGGCCGGGGCCGTTATCCCAGACATCCAGTTGTAGATTTTCGCCCTCAACGTGGGCATCGATACCAATTTCGGCCTGCGCACCCGCATATTTCACCGCGTTCTcc                                  >  2:124662/1‑140 (MQ=255)
                                                                                                    cTTATCAAATATCGTCTGCTCCTGGCCTGGCGGAAGGCCGGGGCCGTTATCCCAGACATCCAGTTGTAGATTTTCGCCCTCAACGTGGGCATCGATACCAATTTCGGCCTGCGCACCCGCATATTTCACCGCGTTCTCca                                 <  1:1085363/140‑1 (MQ=255)
                                                                                                     ttATCAAATATCGTCTGCTCCTGGCCTGGCGGAAGGCCGGGGCCGTTATCCCAGACATCCAGTTGTAGATTTTCGCCCTCAACGTGGGCATCGATACCAATTTCGGCCTGCGCACCCGCATATTTCACCGCGTTCTCcag                                >  1:798386/1‑140 (MQ=255)
                                                                                                       aTCAAATATCGTCTGCTCCTGGCCTGGCGGAAGGCCGGGGCCGTTATCCCAGACATCCAGTTGTAGATTTTCGCCCTCAACGTGGGCATCGATACCAATTTCGGCCTGCGCACCCGCATATTTCACCGCGTTCTCcagca                              <  2:929740/140‑1 (MQ=255)
                                                                                                         cAAATATCGTCTGCTCCTGGCCTGGCGGAAGGCCGGGGCCGTTATCCCAGACATCCAGTTGTAGATTTTCGCCCTCAACGTGGGCATCGATACCAATTTCGGCCTGCGCACCCGCATATTTCACCGCGTTCTCCAGCAGa                            <  2:1147632/140‑1 (MQ=255)
                                                                                                         cAAATATCGTCTGCTCCTGGCCTGGCGGAAGGCCGGGGCCGTTATCCCAGACATCCAGTTGTAGATTTTCGCCCTCAACGTGGGCATCGATACCAATTTCGGCCTGCGCACCCGCATATTTCACCGCGTTCTCCAGCAGa                            <  1:842252/140‑1 (MQ=255)
                                                                                                          aaaTATCGTCTGCTCCTGGCCTGGCGGAAGACCGGGGCCGTTATCCCAGACATCCAGTTGTAGATTTTCGCCCTCAACGTGGGCATCGATACCAATTTCGGCCTGCGCACCCGCATATTTCACCGCGTTCTCCAGCAGAt                           <  1:397168/140‑1 (MQ=255)
                                                                                                          aaaTATCGTCTGCTCCTGGCCTGGCGGAAGACCGGGGCCGTTATCCCAGACATCCAGTTGTAGATTTTCGCCCTCAACGTGGGCATCGATACCAATTTCGGCCTGCGCACCCGCATATTTCACCGCGTTCTCCAGCAGAt                           <  2:230272/140‑1 (MQ=255)
                                                                                                          aaaTATCGTCTGCTCCTGGCCTGGCGGAAGACCGGGGCCGTTATCCCAGACATCCAGTTGTAGATTTTCGCCCTCAACGTGGGCATCGATACCAATTTCGGCCTGCGCACCCGCATATTTCACCGCGTTCTCCAGCAGAt                           <  2:735745/140‑1 (MQ=255)
                                                                                                              atCGTCTGCTCCTGGCCTGGCGGAAGGCCGGGGCCGTTATCCCAGACATCCAGTTGTAGATTTTCGCCCTCAACGTGGGCATCGATACCAATTTCGGCCTGCGCACCCGCATATTTCACCGCGTTCTCCAGCAGATTAAt                       >  1:924699/1‑140 (MQ=255)
                                                                                                                   cTGCTCCTGGCCTGGCGGAAGGCCGGGGCCGTTATCCCAGACATCCAGTTGTAGATTTTCGCCCTCAACGTGGGCATCGATACCAATTTCGGCCTGCGCACCCGCATATTTCATCGCGTTCTCCAGCAGATTAATCAGCa                  <  2:147997/140‑1 (MQ=255)
                                                                                                                   cTGCTCCTGGCCTGGCGGAAGGCCGGGGCCGTTATCCCAGACATCCAGTTGTAGATTTTCGCCCTCAACGTGGGCATCGATACCAATTTCGGCCTGCGCACCCGCATATTTCACCGCGTTCTCCAGCAGATTAATCAGCa                  <  2:697195/140‑1 (MQ=255)
                                                                                                                   cTGCTCCTGGCCTGGCGGAAGGCCGGGGCCGTTATCCCAGACATCCAGTTGTAGATTTTCGCCCTCAACGTGGGCATCGATACCAATTTCGGCCTGCGCACCCGCATATTTCACCGCGTTCTCCAGCAGATTAATCAGCa                  <  1:954604/140‑1 (MQ=255)
                                                                                                                    tGCTCCTGGCCTGGCGGAAGGCCGGGGCCGTTATCCCAGACATCCAGTTGTAGATTTTCGCCCTCAACGTGGGCATCGATACCAATTTCGGCCTGCGCACCCGCATATTTCACCGCGTTCTCCAGCAGATTAATCAGCa                  <  1:1085723/139‑1 (MQ=255)
                                                                                                                    tGCTCCTGGCCTGGCGGAAGGCCGGGGCCGTTATCCCAGACATCCAGTTGTAGATTTTCGCCCTCAACGTGGGCATCGATACCAATTTCGGCCTGCGCACCCGCATATTTCACCGCGTTCTCCAGCAGATTAATCAGCa                  >  2:1085723/1‑139 (MQ=255)
                                                                                                                     gCTCCTGGCCTGGCGGAAGGCCGGGGCCGTTATCCCAGACATCCAGTTGTAGATTTTCGCCCTCAACGTGGGCATCGATACCAATTTCGGCCTGCGCACCCGCATATTTCACCGCGTTCTCCAGCAGATTAATCAGCAcc                <  2:295489/140‑1 (MQ=255)
                                                                                                                      cTCCTGGCCTGGCGGAAGGCCGGGGCCGTTATCCCAGACATCCAGTTGTAGATTTTCGCCCTCAACGTGGGCATCGATACCAATTTCGGCCTGCGCACCCGCATATTTCACCGCGTTCTCCAGCAGATTAATCAGCAccc               >  1:874778/1‑140 (MQ=255)
                                                                                                                                ggcGGAAGGCCGGGGCCGTTATCCCAGACATCCAGTTGTAGATTTTCGCCCTCAACGTGGGCATCGATACCAATTTCGGCCTGCGCACCCGCATATTTCACCGCGTTCTCCAGCAGATTAAt                       <  2:672891/122‑1 (MQ=255)
                                                                                                                                ggcGGAAGGCCGGGGCCGTTATCCCAGACATCCAGTTGTAGATTTTCGCCCTCAACGTGGGCATCGATACCAATTTCGGCCTGCGCACCCGCATATTTCACCGCGTTCTCCAGCAGATTAAt                       >  1:672891/1‑122 (MQ=255)
                                                                                                                                ggcGGAAGGCCGGGGCCGTTATCCCAGACATCCAGTTGTAGATTTTCGCCCTCAACGTGGGCATCGATACCAATTTCGGCCTGCGCACCCGCATATTTCACCGCGTTCTCCAGCAGATTAATCAGCACCCGTTCAAagag     >  2:274749/1‑140 (MQ=255)
                                                                                                                                 gcGGAAGGCCGGGGCCGTTATCCCAGACATCCAGTTGTAGATTTTCGCCCTCAACGTGGGCATCGATACCAATTTCGGCCTGCGCACCCGCATATTTCACCGCGTTCTCCAGCAGATTAATCAGCACCCGTTCAAagag     <  2:145760/139‑1 (MQ=255)
                                                                                                                                 gcGGAAGGCCGGGGCCGTTATCCCAGACATCCAGTTGTAGATTTTCGCCCTCAACGTGGGCATCGATACCAATTTCGGCCTGCGCACCCGCATATTTCACCGCGTTCTCCAGCAGATTAATCAGCACCCGTTCAAagag     >  1:145760/1‑139 (MQ=255)
                                                                                                                                 gcGGAAGACCGGGGCCGTTATCCCAGACATCCAGTTGTAGATTTTCGCCCTCAACGTGGGCATCGATACCAATTTCGGCCTGCGCACCCGCATATTTCACCGCGTTCTCCAGCAGATTAATCAGCACCCGTTCAAag       <  1:450942/137‑1 (MQ=255)
                                                                                                                                 gcGGAAGACCGGGGCCGTTATCCCAGACATCCAGTTGTAGATTTTCGCCCTCAACGTGGGCATCGATACCAATTTCGGCCTGCGCACCCGCATATTTCACCGCGTTCTCCAGCAGATTAATCAGCACCCGTTCAAag       >  2:450942/1‑137 (MQ=255)
                                                                                                                                  cGGAAGGCCGGGGCCGTTATCCCAGACATCCAGTTGTAGATTTTCGCCCTCAACGTGGGCATCGATACCAATTTCGGCCTGCGCACCCGCATATTTCACCGCGTTCTCCAGCAGATTAATCAGCACCCGTTCAAAGAGTg   <  2:798386/140‑1 (MQ=255)
                                                                                                                                   ggAAGACCGGGGCCGTTATCCCAGACATCCAGTTGTAGATTTTCGCCCTCAACGTGGGCATCGATACCAATTTCGGCCTGTGCACCCGCATATTTCACCGCGTTCTCCAGCAGATTAATCAGCACCCGTTCAAAGAGTgg  <  2:852591/140‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                      
ACGCGGTAATAGTGCCCCCGTGTACATCCACTATCGCCCGACAAATTGCCAGTCCAAGCCCTACCCCCGGTACTGCCGACTCTTTATTCCCGCGAGCAAACTTATCAAATATCGTCTGCTCCTGGCCTGGCGGAAGACCGGGGCCGTTATCCCAGACATCCAGTTGTAGATTTTCGCCCTCAACGTGGGCATCGATACCAATTTCGGCCTGCGCACCCGCATATTTCACCGCGTTCTCCAGCAGATTAATCAGCACCCGTTCAAAGAGTGG  >  NC_000913/721819‑722089

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