Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,443,2770AC41.4% 159.9 / 31.0 87V274G (GTA→GGA) secYpreprotein translocase membrane subunit
Reads supporting (aligned to +/- strand):  ref base A (48/3);  new base C (0/36);  total (48/39)
Fisher's exact test for biased strand distribution p-value = 2.50e-21
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCG‑CCCCCGAACCATGACGCGATGGTCGCCGGGAAC‑AGAATAATACTGGAAGCGAAG‑ATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAAATACTAAT  >  NC_000913/3443129‑3443423
                                                                                                                                                        |                                                                                                                                                  
aGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCG‑CCCCCGAACCATGACGCGATGGTCGCCGGGAAC‑AGAATAATACTGGAAGCGAAG‑ATTGCCGGGATTAcc                                                                                                                                                  >  1:968307/1‑151 (MQ=255)
  aCGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCG‑CCCCCGAACCATGACGCGATGGTCGCCGGGAAC‑AGAATAATACTGGAAGCGAAG‑ATTGCCGGGATTAcccc                                                                                                                                                >  2:711805/1‑151 (MQ=255)
   cGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCG‑CCCCCGAACCATGACGCGATGGTCGCCGGGAAC‑AGAATAATACTGGAAGCGAAG‑ATTGCCGGGATTAccccc                                                                                                                                               <  2:240994/151‑1 (MQ=255)
        agagTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCG‑CCCCCGAACCATGACGCGATGGTCGCCGGGAAC‑AGAATAATACTGGAAGCGAAG‑ATTGCCGGGATTACCCCCGCCat                                                                                                                                          >  2:264886/1‑151 (MQ=255)
          agTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCG‑CCCCCGAACCATGACGCGATGGTCGCCGGGAAC‑AGAATAATACTGGAAGCGAAG‑ATTGCCGGGATTACCCCCGCCatat                                                                                                                                        >  1:381378/1‑151 (MQ=255)
           gTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCG‑CCCCCGAACCATGACGCGATGGTCGCCGGGAAC‑AGAATAATACTGGAAGCGAAG‑ATTGCCGGGATTCCCCCCGCCatat                                                                                                                                        <  1:308646/150‑1 (MQ=255)
           gTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCG‑CCCCCGAACCATGACGCGATGGTCGCCGGGAAC‑AGAATAATACTGGAAGCGAAG‑ATTGCCGGGATTACCCCCGCCATAtt                                                                                                                                       >  1:50699/1‑151 (MQ=255)
           gTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCG‑CCCCCGAACCATGACGCGATGGTCGCCGGGAAC‑AGAATAATACTGGAAGCGAAG‑ATTGCCGGGATTACCCCCGCCATAtt                                                                                                                                       >  1:70633/1‑151 (MQ=255)
                acaTAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCG‑CCCCCGAACCATGACGCGATGGTCGCCGGGAAC‑AGAGTAATACTGGAAGCGAAG‑ATTGCCGGGATTACCCCCGCCATATTCACtt                                                                                                                                  >  2:10009/1‑151 (MQ=255)
                acaTAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCG‑CCCCCGAACCATGACGCGATGGTCGCCGGGAAC‑AGAATAATACTGGAAGCGAAG‑ATTGCCGGGATTACCCCCGCCATATTCACtt                                                                                                                                  >  2:402641/1‑151 (MQ=255)
                  aTAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCG‑CCCCCGAACCATGACGCGATGGTCGCCGGGAAC‑AGAATAATACTGGAAGCGAAG‑ATTGCCGGGATTACCCCCGCCATATTCACTTTc                                                                                                                                >  1:682032/1‑151 (MQ=255)
                   aaaggCGGTGCCCCAGGCTAAAAATACAGCAAATTTGTTGTCCGCCAGTCCCAC‑CCAGTCCCC‑CCCCCAAACCAGGCCGCGTGGGTCGCCGGAAAC‑AAAAAAATACTGGAAGCGAAG‑ATGGCGGGGATTCCCCCCGCCATATTCACTTTCa                                                                                                                               <  2:324661/147‑1 (MQ=255)
                      aGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCG‑CCCCCGAACCATGACGCGATGGTCGCCGGGAAC‑AGAATAATACTGGAAGCGAAG‑ATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCg                                                                                                                            >  1:621697/1‑151 (MQ=255)
                             ccccGGGCTAAAAATACAGAGAATTTTTTGCCACCCATTCCCAA‑CAAGCCCCC‑CCCCCAAACCAGGCGGCGGGGGCCGCGGGAAAA‑AAAATAATAGGGGAAGAAAAA‑ATGGGGGGTATTCCCCCCGCCATATTCACTTTCAGCGGTaaa                                                                                                                       <  2:629282/148‑1 (MQ=255)
                              cccagccgcaaaaCACAGCAAATTTGTTGCCCCCCTTTCCCAC‑CCAGTCCCC‑CCCCCAACCCATGACGCGTTGGCCCCCGGAAAA‑AAAAAAATACGGGAAGAGAAA‑ATGGCGGGTATTCCCCCCGCCATATTCACTTTCAGCGGTaaa                                                                                                                       <  2:689816/138‑1 (MQ=255)
                               ccaggtaacAAATACAGCAAATTTGTTGTCAGCAATTTCAAC‑CAAGTCCCC‑CCCCCAAACCAGGACGCGAGGGTCGCGGGAAAA‑AAATTAATACGGGAAGCGAAG‑ATGGCGGGAATTACCCCCGCCATATTCACTTTCAGCGGTAAAtgtg                                                                                                                   <  1:1132552/143‑1 (MQ=255)
                                        aaaaaaagggAATTTTTTGCCAGCCATTCCCAC‑CCAGTCCCG‑CCCCCAACCCAGGACGCGATGGTCGCGGGGAAA‑AAAATAATTCTGGAAGTAAAA‑ATTGCGGGGATACCCCCCGCCATATTCACTTTCAGCGGTAAAtgtgt                                                                                                                  <  2:1139167/134‑1 (MQ=255)
                                              agaaaatttgttgCCAGCCAGTCCCAA‑CCAGTCCCG‑CCCCCGAACCATGACGCGATGGTCGCCGGAAAC‑AAAAAAATACTGGAAGCGAAG‑TTGGCCGGTATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTgcagc                                                                                                      <  1:639000/142‑1 (MQ=255)
                                                 gAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCG‑CCCCCGAACCATGACGCGATGGTCGCCGGGAAC‑AGAATAATACTGGAAGCGAAG‑ATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCAt                                                                                                    >  1:1219557/1‑148 (MQ=255)
                                                  aattttttgtCAGCCAGTCCCAA‑CCAGTCCCC‑CCCCCGAACCATGCCGCGATGGTCGCGGGAAAC‑AAAATAATACGGGAAGCGAAG‑ATTGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGa                                                                                                 <  1:1115125/144‑1 (MQ=255)
                                                    attttttgtcgCCCGTTCCAA‑CCCGGCCCCCCCCCCGAACCATGACGCGATGGTCGCGGGGAAA‑AAAATAATAGGGGAGGCGAAG‑AGGGCGGGATTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGac                                                                                                <  2:70651/140‑1 (MQ=255)
                                                    aTTGTTGTCAGCCAGTTCCAA‑CCAGTACCG‑CCCCCGAACCATGACGCGATGGTCGCCGGGAAC‑AGAATAATACTGGAAGCGAAG‑ATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGacac                                                                                              >  1:453155/1‑151 (MQ=255)
                                                      ttgttgccgCCCGTTCCAACCCGGTCCCC‑CCCCCAACCCAGGACGGGTGGGTCGCGGGAAAA‑AAAAAAATAGTGGAAGGAAAG‑ATTGCGGGGTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTgcagca                                                                                                     <  2:1116620/135‑1 (MQ=255)
                                                        ttgcCAGCCAGTCCCAC‑CAAGTCCCC‑CCCCCGAACCATGACGCGTTGGTCGCGGGAAAA‑AAAATAATCCTGAAAGCGAAG‑ATTGCGGGGATCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACacgac                                                                                           <  1:300387/150‑1 (MQ=255)
                                                           ccgcccAGTCCCAC‑CCATCCCCC‑CCCCCGAACCAGGACGGGATGGTCGCGGGAAAA‑AAAATAATACTGGAAGCGAAG‑TTTGCGGGGTTCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACAcgacga                                                                                         <  2:195581/145‑1 (MQ=255)
                                                            cAGCCAGTTCCAA‑CCAGTACCG‑CCCCCGAACCATGACGCGATGGTCGCCGGGAAC‑AGAATAATACTGGAAGCGAAG‑ATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCt                                                                                      >  1:953475/1‑151 (MQ=255)
                                                             aGCCAGTTCCAA‑CCAGTACCG‑CCCCCGAACCATGACGCGATGGTCGCCGGGAAC‑AGAATAATACTGGAAGCGAAG‑ATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCtt                                                                                     >  1:921520/1‑151 (MQ=255)
                                                               ccAGTTCCAA‑CCAGTACCG‑CCCCCGAACCATGACGCGATGGTCGCCGGGAAC‑AGAATAATACTGGAAGCGAAG‑ATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGc                                                                                   >  1:676255/1‑151 (MQ=255)
                                                               ccAGTTCCAA‑CCAGTACCG‑CCCCCGAACCATGACGCGATGGTCGCCGGGAAC‑AGAATAATACTGGAAGCGAAG‑ATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGc                                                                                   >  2:48110/1‑151 (MQ=255)
                                                                 aGTTCCAA‑CCAGTACCG‑CCCCCGAACCATGACGCGATGGTCGCCGGGAAC‑AGAATAATACTGGAAGCGAAG‑ATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTg                                                                                 >  1:339714/1‑151 (MQ=255)
                                                                 aGTTCCAA‑CCAGTACCG‑CCCCCGAACCATGACGCGATGGTCGCCGGGAAC‑AGAATAATACTGGAAGCGAAG‑ATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTg                                                                                 >  1:1007399/1‑151 (MQ=255)
                                                                 aGTTCCAA‑CCAGTACCG‑CCCCCGAACCATGACGCGATGGTCGCCGGGAAC‑AGAATAATACTGGAAGCGAAG‑ATTGCCGGGATTACACCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTg                                                                                 >  1:374124/1‑151 (MQ=255)
                                                                  attaccaacagtcaCCC‑CCCCCCACCCATGTCGCGAGGGTC‑CCGGGAACAAGAAAAATCCGGGAGGCAAAT‑TTGGCGGGTATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCt                                                                                  <  2:106535/137‑1 (MQ=255)
                                                                     ccAA‑CCAGCCCCG‑CCCCAAACGCGTGACGCGATGGCCGCGGGAAAC‑AAAAAAATGCGGGAAGCAAGG‑ATTGCGGGGTTTCCCGCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTg                                                                                    <  1:189968/144‑1 (MQ=255)
                                                                      caagccagccccG‑CCCCCGCACCGGGGCGCAAGGGGCGCGAGGAGA‑AGAAAAATAATGG‑AGGGAAGTTTTGCTGGAATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGAc                                                                               <  1:636251/139‑1 (MQ=255)
                                                                        a‑CCAGTCCCC‑CCCCAGAACAATGCCGCATGGGTCGCAGGAAAC‑AGAAAAATATTGGAAGGGAAG‑ATTGCCGGATTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTc                                                                          <  2:978050/151‑1 (MQ=255)
                                                                          ccagccccG‑CCCCCAACCCATGCCGGGATGGTCGCGGGAAAA‑AAAATAATACTGGAAGCAAAG‑ATTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTc                                                                          <  1:310729/144‑1 (MQ=255)
                                                                          ccAGTCCCG‑CCCCCGAACCATGACGCGATGGTCGCCGGGAAC‑AGAATAATACTGGAAGAGAAG‑ATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTc                                                                          <  2:1076250/150‑1 (MQ=255)
                                                                          ccAGTCCCC‑CCCCAGACCCATGAGGGGATGGTCGCGGGAAAA‑AAAATAAAACTGGAAGCGAAA‑ATTGCGGGATTTCCCCCCGCCATATTCACTTTCAGCTGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTc                                                                          <  2:639731/150‑1 (MQ=255)
                                                                           cAGTACCG‑CCCCCGAACCATGACGCGATGGTCGCCGGGAAC‑AGAATAATACTGGAAGCGAAG‑ATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTcgc                                                                        >  1:715737/1‑151 (MQ=255)
                                                                              tACCG‑CCCCCGAACCATGACGCGATGGTCGCCGGGAAC‑AGAATAATACTGGAAGCGAAG‑ATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGt                                                                      >  1:876779/1‑150 (MQ=255)
                                                                                cccccccccAAACCATGGCGGGAGGGTCCCCGGAAAC‑AAAAAAATACTGGAAGAAAAG‑TTGGCGGGTATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTa                                                                     <  2:936377/146‑1 (MQ=255)
                                                                                       ccGAACCATGACGCGATGGTCGCCGGGAAC‑AGAATAATACTGGAAGCGAAG‑ATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTAcc                                                              >  2:113097/1‑150 (MQ=255)
                                                                                          aaCCATGACGCGATGGTCGCCGGGAAC‑AGAATAATACTGGAAGCGAAG‑ATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACcaca                                                           >  2:848027/1‑150 (MQ=255)
                                                                                          aaCCATGACGCGATGGTCGCCGGGAAC‑AGAATAATACTGGAAGCGAAG‑ATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACcaca                                                           >  1:1088517/1‑150 (MQ=255)
                                                                                                  cgcgATGGTCGCCGGGAAC‑AGAATAATACTGGAAGCGAAG‑ATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGt                                                  >  1:1063001/1‑151 (MQ=255)
                                                                                                   gggATGGTCGCGGGAAAA‑AAAAAAATACTGGAAGGGAAG‑ATTGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATgc                                                       <  1:109003/143‑1 (MQ=255)
                                                                                                   gcgATGGTCGCCGGGAAC‑AGAATAATACTGGAAGCGAAG‑ATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGtt                                                 >  1:7078/1‑151 (MQ=255)
                                                                                                     gtgggTCGCGGGAAAC‑AAAAAAATACGGAAAGCGAAG‑TTGGCCGGGTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGtt                                                 <  1:1102229/146‑1 (MQ=255)
                                                                                                     gATGGTCGCCGGGAAC‑AGAATAATACTGGAAGCGAAG‑ATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTg                                                >  1:647701/1‑150 (MQ=255)
                                                                                                       tGGTCGCCGGGAAC‑AGAATAATACTGGAAGCGAAG‑ATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGAcc                                             >  1:410816/1‑151 (MQ=255)
                                                                                                       tGGTCGCCGGGAAC‑AGAATAATACTGGAAGCGAAG‑ATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGAcc                                             >  2:106725/1‑151 (MQ=255)
                                                                                                        ggTCGCGGGAAAC‑AAAATAATACGGGAAGCAAAG‑ATTGCGGGATTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCa                                            <  1:3540/151‑1 (MQ=255)
                                                                                                             gcgggaaaaaaaaaaataCGGAAAGCGAAG‑TTTGCCGGGTTCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTg                                                <  2:192921/129‑1 (MQ=255)
                                                                                                               ggaaaC‑AAAATAATACGGGAAGCGAAG‑TTTGCCGGTATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCAc                                           <  2:333946/142‑1 (MQ=255)
                                                                                                                ggaaaaaaaataatACTGGAAGCAAAG‑ATTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTc                                       <  1:855599/141‑1 (MQ=255)
                                                                                                                    c‑AGAATAATACTGGAAGCGAAG‑ATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAAt                                >  1:787222/1‑151 (MQ=255)
                                                                                                                    c‑AGAATAATACTGGAAGCGAAG‑ATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAAt                                >  1:1151665/1‑151 (MQ=255)
                                                                                                                    c‑AGAATAATACTGGAAGCGAAG‑ATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAAt                                >  1:851935/1‑151 (MQ=255)
                                                                                                                      ataataatACGTGAAGCAAAG‑ATTGCCGGCATGCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATa                               <  2:47715/149‑1 (MQ=255)
                                                                                                                      aaaataatACTGGAAGCAAAG‑TTTGCCGGAATCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAAt                                <  2:403823/148‑1 (MQ=255)
                                                                                                                      aaaaaattccTGGAAGCGAAG‑ATGGCGGGGATCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTcaaca                                   <  2:234500/138‑1 (MQ=255)
                                                                                                                      aGAATAATACTGGAAGCGAAG‑ATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAAt                                >  2:613931/1‑150 (MQ=255)
                                                                                                                        aataatACGGGAAGCGAAG‑ATTGCGGGGATCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATaca                             <  2:1062983/151‑1 (MQ=255)
                                                                                                                         ataataCTGGAAGCGAAG‑ATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATaca                             >  1:497114/1‑150 (MQ=255)
                                                                                                                           atgccggaaaGCGAAG‑ATGGCCAGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATacaaca                          <  1:653557/143‑1 (MQ=255)
                                                                                                                           aataCTGGAAGCGAAG‑ATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATacaac                           >  2:31604/1‑150 (MQ=255)
                                                                                                                                 ggAAGCGAAG‑ATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAAc                    >  1:400689/1‑151 (MQ=255)
                                                                                                                                 ggAAGAGAAG‑ATTGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGaa                     <  1:1057665/150‑1 (MQ=255)
                                                                                                                                 gaaaGCGAAG‑ATTGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAAc                    <  2:42890/149‑1 (MQ=255)
                                                                                                                                   aaGCGAAG‑ATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGt                  >  2:24307/1‑151 (MQ=255)
                                                                                                                                   aaGCGAAG‑ATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGt                  >  2:371638/1‑151 (MQ=255)
                                                                                                                                   aaGCAAAG‑ATTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGt                  <  1:192599/151‑1 (MQ=255)
                                                                                                                                   aaGAGATG‑ATTGCCGGGTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGaa                     <  2:371549/148‑1 (MQ=255)
                                                                                                                                     gCGAAG‑ATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTc                 >  2:1127973/1‑150 (MQ=255)
                                                                                                                                     gCGAAG‑ATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTc                 >  2:620793/1‑150 (MQ=255)
                                                                                                                                      cGAAG‑ATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCAc               >  1:523824/1‑151 (MQ=255)
                                                                                                                                      cGAAG‑ATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCAc               >  1:1169251/1‑151 (MQ=255)
                                                                                                                                       gAAG‑ATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCAc               >  2:1148671/1‑150 (MQ=255)
                                                                                                                                        aaG‑ATTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACt              <  1:465475/150‑1 (MQ=255)
                                                                                                                                         aG‑ATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGc            >  2:819330/1‑151 (MQ=255)
                                                                                                                                         aG‑ATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGc            >  1:152030/1‑151 (MQ=255)
                                                                                                                                          g‑ATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGc            >  2:742544/1‑150 (MQ=255)
                                                                                                                                               gCCGGGTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCaa          <  2:319238/148‑1 (MQ=255)
                                                                                                                                                  gggATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaatact     >  1:313157/1‑150 (MQ=255)
                                                                                                                                                   ggATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaatactaa   <  2:1031240/151‑1 (MQ=255)
                                                                                                                                                     aTTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaatactaat  >  1:687809/1‑150 (MQ=255)
                                                                                                                                                        |                                                                                                                                                  
AGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCG‑CCCCCGAACCATGACGCGATGGTCGCCGGGAAC‑AGAATAATACTGGAAGCGAAG‑ATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAAATACTAAT  >  NC_000913/3443129‑3443423

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑