Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 528,763 T→C intergenic (+409/+106) b0501 → / ← b0502 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913528,7630TC88.5% 531.3 / 53.2 209intergenic (+409/+106)b0501/b0502hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (10/14);  new base C (93/92);  total (103/106)
Fisher's exact test for biased strand distribution p-value = 5.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.58e-01

Only 100 of 210 total aligned reads displayed.

CCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTCA  >  NC_000913/528671‑528853
                                                                                            |                                                                                          
ccTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGt                                                                                            >  2:4346986/1‑93 (MQ=255)
 cTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTa                                                                                           >  1:2799888/1‑93 (MQ=255)
  ttCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGCAc                                                                                          >  2:6484871/1‑93 (MQ=255)
    ccGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGCAcgg                                                                                        <  2:2479490/93‑1 (MQ=255)
    ccGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGCAcgg                                                                                        >  1:1535115/1‑93 (MQ=255)
     cGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGCAcggc                                                                                       >  1:365810/1‑93 (MQ=255)
     cGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGCAcggc                                                                                       >  2:2337815/1‑93 (MQ=255)
     cGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGCAcggc                                                                                       >  2:4650396/1‑93 (MQ=255)
       gATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGCAcggcgg                                                                                     >  1:1141376/1‑93 (MQ=255)
       gATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGCAcggcgg                                                                                     <  1:5823681/93‑1 (MQ=255)
         ttAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGCACGGCGGcc                                                                                   >  1:1893467/1‑93 (MQ=255)
          tAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGCACGGCGGccc                                                                                  >  1:1789553/1‑93 (MQ=255)
           aaaaTAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGCACGGCGGCCCg                                                                                 >  1:974754/1‑93 (MQ=255)
            aaataaatTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGCACGGCGGCCCGa                                                                                <  2:789730/93‑1 (MQ=255)
            aaataaatTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGCACGGCGGCCCGa                                                                                >  1:5996466/1‑93 (MQ=255)
             aataaatTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGCACGGCGGCCCGAt                                                                               <  2:4631691/93‑1 (MQ=255)
             aataaatTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGCACGGCGGCCCGAt                                                                               >  1:5870205/1‑93 (MQ=255)
              ataaatTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATa                                                                              <  1:4731848/93‑1 (MQ=255)
              ataaatTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGCACGGCGGCCCGATa                                                                              <  1:3423445/93‑1 (MQ=255)
              ataaatTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGCACGGCGGCCCGATa                                                                              <  1:1034619/93‑1 (MQ=255)
               taaatTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAg                                                                             >  2:5483169/1‑93 (MQ=255)
                 aatTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGCACGGCGGCCCGATAGTc                                                                           >  1:2315153/1‑93 (MQ=255)
                  atTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTcc                                                                          >  2:2566843/1‑93 (MQ=255)
                   ttCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGCACGGCGGCCCGATAGTCCa                                                                         >  2:5942050/1‑93 (MQ=255)
                   ttCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGCACGGCGGCCCGATAGTCCa                                                                         <  1:3791263/93‑1 (MQ=255)
                    tCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGCACGGCGGCCCGATAGTCCAg                                                                        <  1:4639527/93‑1 (MQ=255)
                    tCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGCACGGCGGCCCGATAGTCCAg                                                                        <  1:3948634/93‑1 (MQ=255)
                      ttGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAggg                                                                      <  1:7387743/93‑1 (MQ=255)
                      ttGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGCACGGCGGCCCGATAGTCCAggg                                                                      >  2:1343173/1‑93 (MQ=255)
                       tGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGt                                                                     <  2:7033382/93‑1 (MQ=255)
                        gCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGtt                                                                    <  1:2132745/93‑1 (MQ=255)
                           gggAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGCACGGCGGCCCGATAGTCCAGGGTTCAg                                                                 >  2:944000/1‑93 (MQ=255)
                            ggAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGCACGGCGGCCCGATAGTCCAGGGTTCAGc                                                                >  2:5758725/1‑93 (MQ=255)
                            ggAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGCACGGCGGCCCGATAGTCCAGGGTTCAGc                                                                >  1:6332647/1‑93 (MQ=255)
                            ggAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGCACGGCGGCCCGATAGTCCAGGGTTCAGc                                                                >  1:7065955/1‑93 (MQ=255)
                             gaaagaaaGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGCACGGCGGCCCGATAGTCCAGGGTTCAGCg                                                               <  2:5170582/93‑1 (MQ=255)
                             gaaagaaaGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGCACGGCGGCCCGATAGTCCAGGGTTCAGCg                                                               >  2:2675810/1‑93 (MQ=255)
                             gaaagaaaGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGCACGGCGGCCCGATAGTCCAGGGTTCAGCg                                                               <  1:7004002/93‑1 (MQ=255)
                                agaaagaGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGCACGGCGGCCCGATAGTCCAGGGTTCAGCgggg                                                            <  2:1123318/93‑1 (MQ=255)
                                agaaagaGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGCACGGCGGCCCGATAGTCCAGGGTTCAGCgggg                                                            >  1:7462130/1‑93 (MQ=255)
                                 gaaagaGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGCACGGCGGCCCGATAGTCCAGGGTTCAGCggggg                                                           <  1:1165931/93‑1 (MQ=255)
                                 gaaagaGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGCACGGCGGCCCGATAGTCCAGGGTTCAGCggggg                                                           >  2:3006632/1‑93 (MQ=255)
                                 gaaagaGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGCACGGCGGCCCGATAGTCCAGGGTTCAGCggggg                                                           <  1:3368637/93‑1 (MQ=255)
                                   aagaGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCTCCAGCACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTa                                                         >  1:4891255/1‑93 (MQ=255)
                                   aagaGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGCACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTa                                                         >  1:2957956/1‑93 (MQ=255)
                                    agagGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGCACGGCGGCCCGATAGTCCAGGGATCAGCGGGGGAAc                                                        >  2:346980/1‑93 (MQ=255)
                                     gagGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGGACt                                                       >  1:2599332/1‑93 (MQ=255)
                                       ggAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGa                                                     <  1:902774/93‑1 (MQ=255)
                                       ggAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGCACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGa                                                     <  1:5037535/93‑1 (MQ=255)
                                         aaaTAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGCACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATa                                                   >  2:695400/1‑93 (MQ=255)
                                         aaaTAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGCACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATa                                                   <  2:5596661/93‑1 (MQ=255)
                                          aaTAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGCACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAg                                                  <  2:2788043/93‑1 (MQ=255)
                                           aTAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGt                                                 <  1:5662786/93‑1 (MQ=255)
                                           aTAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGCACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGt                                                 >  1:2288890/1‑93 (MQ=255)
                                            tAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGCACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTc                                                <  1:6945750/93‑1 (MQ=255)
                                              aaCCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGCACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTcgc                                              >  1:6489587/1‑93 (MQ=255)
                                              aaCCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGCACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTcgc                                              >  2:3408254/1‑93 (MQ=255)
                                              aaCCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGCACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTcgc                                              <  1:1446569/93‑1 (MQ=255)
                                                   ttAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGCACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGa                                         >  1:2997445/1‑93 (MQ=255)
                                                   ttAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGCACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGa                                         <  1:1224902/93‑1 (MQ=255)
                                                    tAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGCACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGAc                                        <  2:4845959/93‑1 (MQ=255)
                                                    tAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGCACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGAc                                        <  2:184633/93‑1 (MQ=255)
                                                      gCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGCACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTg                                      <  2:5470657/93‑1 (MQ=255)
                                                        ggAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGCACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCa                                    >  2:3301870/1‑93 (MQ=255)
                                                         gAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCAt                                   >  1:7427318/1‑93 (MQ=255)
                                                         gAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGCACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCAt                                   <  2:5427530/93‑1 (MQ=255)
                                                          aaaaCCAGCGGCACGCCAGGGTGACATGCCCCAGCACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATc                                  <  2:4797665/93‑1 (MQ=255)
                                                          aaaaCCAGCGGCACGCCAGGGTGACATGCCCCAGCACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATc                                  <  2:4259702/93‑1 (MQ=255)
                                                           aaaCCAGCGGCACGCCAGGGTGACATGCCCCAGCACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCg                                 <  2:2648252/93‑1 (MQ=255)
                                                           aaaCCAGCGGCACGCCAGGGTGACATGCCCCAGCACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCg                                 >  1:3545168/1‑93 (MQ=255)
                                                            aaCCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCga                                >  1:814448/1‑93 (MQ=255)
                                                            aaCCAGCGGCACGCCAGGGTGACATGCCCCAGCACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCga                                <  2:4502395/93‑1 (MQ=255)
                                                             aCCAGCGGCACGCCAGGGTGACATGCCCCAGCACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCgag                               >  1:4241717/1‑93 (MQ=255)
                                                              ccAGCGGCACGCCAGGGTGACATGCCCCAGCACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCgaga                              <  1:160353/93‑1 (MQ=255)
                                                                aGCGGCACGCCAGGGTGACATGCCCCAGCACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAAc                            <  1:2015574/93‑1 (MQ=255)
                                                                 gCGGCACGCCAGGGTGACATGCCCCAGCACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACt                           >  2:131143/1‑93 (MQ=255)
                                                                  cGGCACGCCAGGGTGACATGCCCCAGCACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTc                          <  2:6315984/93‑1 (MQ=255)
                                                                     cACGCCAGGGTGACATGCCCCAGCACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCAt                       <  2:3817526/93‑1 (MQ=255)
                                                                     cACGCCAGGGTGACATGCCCCAGCACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCAt                       <  2:4911137/93‑1 (MQ=255)
                                                                       cGCCAGGGTGACATGCCCCAGCACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATAc                     <  2:5273143/93‑1 (MQ=255)
                                                                        gCCAGGGTGACATGCCCCAGCACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACa                    >  2:6846600/1‑93 (MQ=255)
                                                                         ccAGGGTGACATGCCCCAGCACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACaa                   >  2:6691084/1‑93 (MQ=255)
                                                                         ccAGGGTGACATGCCCCAGCACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACaa                   <  1:1701693/93‑1 (MQ=255)
                                                                          cAGGGTGACATGCCCCAGCACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAAc                  >  2:5045/1‑93 (MQ=255)
                                                                            gggTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACtt                >  2:27732/1‑93 (MQ=255)
                                                                            gggTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACtt                <  1:4037163/93‑1 (MQ=255)
                                                                            gggTGACATGCCCCAGCACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACtt                <  2:5139274/93‑1 (MQ=255)
                                                                             ggTGACATGCCCCAGCACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACttt               >  2:7406426/1‑93 (MQ=255)
                                                                               tGACATGCCCCAGCACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACttttt             <  1:3141917/93‑1 (MQ=255)
                                                                                gACATGCCCCAGCACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTa            <  2:5375661/93‑1 (MQ=255)
                                                                                 aCATGCCCCAGCACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTAc           <  2:2138154/93‑1 (MQ=255)
                                                                                 aCATGCCCCAGCACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTAc           <  1:6667597/93‑1 (MQ=255)
                                                                                     gCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGaa       <  2:1199797/93‑1 (MQ=255)
                                                                                     gCCCCAGCACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGaa       <  1:2847405/93‑1 (MQ=255)
                                                                                     gCCCCAGCACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGaa       <  2:4018744/93‑1 (MQ=255)
                                                                                     gCCCCAGCACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGaa       <  1:4287903/93‑1 (MQ=255)
                                                                                     gCCCCAGCACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGaa       <  1:4106974/93‑1 (MQ=255)
                                                                                       cccAGCACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACg     >  1:6169262/1‑93 (MQ=255)
                                                                                        ccAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGt    >  1:4603680/1‑93 (MQ=255)
                                                                                          aGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTCa  <  2:5721195/93‑1 (MQ=255)
                                                                                            |                                                                                          
CCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTCA  >  NC_000913/528671‑528853

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