Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,873,917 G→C 6.5% H26D (CAT→GAT) ‡ cysC ← adenosine 5'‑phosphosulfate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,873,9170GC 6.5% 589.1 / 26.5 175H26D (CAT→GAT) ‡cysCadenosine 5'‑phosphosulfate kinase
Reads supporting (aligned to +/- strand):  ref base G (78/83);  new base C (4/7);  total (84/91)
Fisher's exact test for biased strand distribution p-value = 5.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.08e-02

Only 100 of 175 total aligned reads displayed.

CTGGCGTTCGGCGCGGTGTGGCGAGATAAATGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCAGTTCGAATGCACTGAATTCAGATGGCGCAGCAGTTGCCTGGCTAACTGGCTCGTGCACCATACCGGCACCCACGGTCACATTGCTCAGGCGATCGATAAAAATCAGCCCACCCGTCACCGGATTT  >  NC_000913/2873645‑2874187
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                              
           cgcgGTGTGGCGAGATAAATGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCTAACATTGTCGCCATCCAGCATATACGTACTTACGCCGACATTATGTAACGCCTCCTCCAGCTCCCCGGCGACCGTTGATTTACCGGACCCGGATACGACGGTACACCACATCCCTACACCACGAAGACCGTTGTGTAGCTcgc                                                                                                                                                                                                                                                               >  1:27071/1‑279 (MQ=255)
             cgGTGGGGCGAGAAAAATGCGTCCCGCACCACNNNTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTcgcg                                                                                                                                                                                                                                                              <  3:217295/278‑1 (MQ=255)
              ggTGTGGCGAGATAAATGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTcg                                                                                                                                                                                                                                                                <  2:677351/275‑1 (MQ=255)
                  tgGCGAGATAAATGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCgtt                                                                                                                                                                                                                                                            >  2:36162/1‑275 (MQ=255)
                          tAAATGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGa                                                                                                                                                                                                                                                                               >  1:48545/1‑248 (MQ=255)
                                   tCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGAAGCCGAGTTTCTGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAAAGGCCGGTAAACCACAGCACTACACCAAGATAACCGTGGTGTAGCCCGACTTGTTGCCCAGTGCCCGGa                                                                                                                                                                                                                                           >  4:679217/1‑275 (MQ=255)
                                     aGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTANTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACNCCACCACACCCCGCTGACCCGGGTGNAGCCCGCGGTGTTTCACCGTGACCggcagggc                                                                                                                                                                                                                                      >  3:132560/1‑272 (MQ=255)
                                               gTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAg                                                                                                                                                                                                                               >  2:412392/1‑275 (MQ=255)
                                                 ccGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTTTAGCTCGCGTTGTTGCACAGTGCCCGGATGGCTATGCCAgac                                                                                                                                                                                                                             >  2:578042/1‑275 (MQ=255)
                                                            catacattcccccACCTCAACGCAGCGACGGATNNNTTTTTACCGACCGGCCCCGGTAAAACCGAGATCGCTGCATCATCCGTGGCGAACCTCGTCGCCCTCCAGCAGCTACGTACTGCCGCCGAGGTTATGTACCGCCTCCTCCAGCGCCCCTGCGACCGTTGATTTACCGGTCCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATg                                                                                                                                                                                                              <  3:224789/268‑1 (MQ=255)
                                                              tCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCAACGCACTACACCCCGATGACCGTGGGGTaa                                                                                                                                                                                                                                                                    >  4:239016/1‑222 (MQ=255)
                                                                    tCGCCACTGCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAg                                                                                                                                                                                                           <  2:733928/274‑1 (MQ=255)
                                                                        catttCACCGACGTAACGGATATTCTCTTTACGATCGGAATCGCTAAAACCGAGATCGCTGCATAAGCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCa                                                                                                                                                                                                      <  4:319516/272‑1 (MQ=255)
                                                                                            ttttCTCTTTCCGATCGTCATCGCTAAACCCTAGATCCCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCCGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGTCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATcgcg                                                                                                                                                                              <  1:49399/277‑1 (MQ=255)
                                                                                                                        aCCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGAt                                                                                                                                                                                                                                          >  4:514639/1‑191 (MQ=255)
                                                                                                                         ccaagaTCGCTGCATAATCCGTGGCTAACTTTGTTGCCATCAAGCGGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGAGCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATTACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGaa                                                                                                                                                 <  1:104360/276‑1 (MQ=255)
                                                                                                                            agaTCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGa                                                                                                                                                                                                                                                >  4:754093/1‑181 (MQ=255)
                                                                                                                            agaTCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCt                                                                                                                                                                                                                                      <  2:514639/191‑1 (MQ=255)
                                                                                                                             gaTCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGa                                                                                                                                                                                                                                                                               >  6:274263/1‑149 (MQ=255)
                                                                                                                                cGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCgg                                                                                                                                                                                                                                            <  2:754093/181‑1 (MQ=255)
                                                                                                                                 gCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCgt                                                                                                                                                                                                                                                                           >  5:1030146/1‑149 (MQ=255)
                                                                                                                                 gCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCg                                                                                                                                                                                                                                                                            <  7:764958/148‑1 (MQ=255)
                                                                                                                                  cTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACGCCGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCgt                                                                                                                                                                                                                                                                           >  8:688466/1‑148 (MQ=255)
                                                                                                                                   tGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCgtgg                                                                                                                                                                                                                                                                         <  6:441313/149‑1 (MQ=255)
                                                                                                                                    gCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGGTCATCgtggt                                                                                                                                                                                                                                                                        <  6:274846/149‑1 (MQ=255)
                                                                                                                                    gCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGGTCAACgtggt                                                                                                                                                                                                                                                                        <  6:763130/149‑1 (MQ=255)
                                                                                                                                    gCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCgtggt                                                                                                                                                                                                                                                                        <  6:6437/149‑1 (MQ=255)
                                                                                                                                       tAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACACCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGGTCATCGTGGTGTa                                                                                                                                                                                                                                                                     <  6:400003/149‑1 (MQ=255)
                                                                                                                                        aaTCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTa                                                                                                                                                                                                                                                                     >  7:306512/1‑148 (MQ=255)
                                                                                                                                        aaTCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTa                                                                                                                                                                                                                                                                     >  7:435997/1‑148 (MQ=255)
                                                                                                                                         aTCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGGTCATCGTGGTGTAg                                                                                                                                                                                                                                                                    <  8:492319/148‑1 (MQ=255)
                                                                                                                                         aTCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGc                                                                                                                                                                                                                                                                   <  5:1077860/149‑1 (MQ=255)
                                                                                                                                           ccGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTc                                                                                                                                                                                                                                                                 <  6:868161/149‑1 (MQ=255)
                                                                                                                                           ccGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTc                                                                                                                                                                                                                                                                 >  6:105653/1‑149 (MQ=255)
                                                                                                                                               ggCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCACCAGAGCGGAAAGTGGCGACGACCCAGATCATTCAAGTcaa                                                                                                                               >  2:512735/1‑273 (MQ=255)
                                                                                                                                                gCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCgtt                                                                                                                                                                                                                                                            >  6:128681/1‑149 (MQ=255)
                                                                                                                                                   aaCATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCgttgtt                                                                                                                                                                                                                                                         <  5:596729/149‑1 (MQ=255)
                                                                                                                                                    aCATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGttgttg                                                                                                                                                                                                                                                        >  6:52531/1‑149 (MQ=255)
                                                                                                                                                        tGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGcac                                                                                                                                                                                                                                                     <  7:521461/148‑1 (MQ=255)
                                                                                                                                                          tCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGt                                                                                                                                                                                                                                                  >  5:552947/1‑149 (MQ=255)
                                                                                                                                                            gCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTg                                                                                                                                                                                                                                                 >  8:347440/1‑148 (MQ=255)
                                                                                                                                                             ccATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGa                                                                                                                                                                                                                                                <  8:435997/148‑1 (MQ=255)
                                                                                                                                                              cATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGAc                                                                                                                                                                                                                                               >  7:774004/1‑148 (MQ=255)
                                                                                                                                                                  cagcagATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGGTCATCGTGGTGTAgtgctgt                                                                                                                                                                                                                                                              >  6:169481/1‑123 (MQ=39)
                                                                                                                                                                  cagcagATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGa                                                                                                                                                                                                                                           <  8:400233/148‑1 (MQ=255)
                                                                                                                                                                     cagATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCAGTTCGAATGCACTGAATTCAGa                                                                                                         <  2:662351/275‑1 (MQ=255)
                                                                                                                                                                      agATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCAGTTCGAATGCACTGAATTCAGATGgcg                                                                                                    <  1:132560/279‑1 (MQ=255)
                                                                                                                                                                            gTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCa                                                                                                                                                                                                                                >  5:438875/1‑149 (MQ=255)
                                                                                                                                                                                tGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGa                                                                                                                                                                                                                                                >  6:958521/1‑129 (MQ=255)
                                                                                                                                                                                tGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGa                                                                                                                                                                                                                                                <  5:958521/129‑1 (MQ=255)
                                                                                                                                                                                  aCGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAgacgac                                                                                                                                                                                                                          <  6:438875/149‑1 (MQ=255)
                                                                                                                                                                                    gCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGt                                                                                                                                                                                                                        >  5:81952/1‑149 (MQ=255)
                                                                                                                                                                                    gCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGt                                                                                                                                                                                                                        >  6:603964/1‑149 (MQ=255)
                                                                                                                                                                                       gAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCa                                                                                                                                                                                                                                >  4:735757/1‑138 (MQ=255)
                                                                                                                                                                                        aGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATgg                                                                                                                                                                                                                                        >  2:341265/1‑129 (MQ=255)
                                                                                                                                                                                        aGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGtttt                                                                                                                                                                                                                     <  8:119104/148‑1 (MQ=255)
                                                                                                                                                                                        aGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGtttt                                                                                                                                                                                                                     <  8:306512/148‑1 (MQ=255)
                                                                                                                                                                                        aGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTc                                                                                                                                                                                                                    <  6:316233/149‑1 (MQ=255)
                                                                                                                                                                                         gTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCg                                                                                                                                                                                                                   >  5:950524/1‑149 (MQ=255)
                                                                                                                                                                                          tttATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGt                                                                                                                                                                                                                  >  5:653830/1‑149 (MQ=255)
                                                                                                                                                                                           ttATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTc                                                                                                                                                                                                                 <  5:173129/149‑1 (MQ=255)
                                                                                                                                                                                            tATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTAt                                                                                                                                                                                                                                    <  4:341265/129‑1 (MQ=255)
                                                                                                                                                                                             aTGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCAt                                                                                                                                                                                                               >  5:271772/1‑149 (MQ=255)
                                                                                                                                                                                                  aCGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAgc                                                                                                                                                                                                          >  5:540058/1‑149 (MQ=255)
                                                                                                                                                                                                  aCGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAg                                                                                                                                                                                                           >  7:271974/1‑148 (MQ=255)
                                                                                                                                                                                                     cctcctCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATg                                                                                                                                                                                                                                   >  8:530251/1‑121 (MQ=255)
                                                                                                                                                                                                     cctcctCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATg                                                                                                                                                                                                                                   <  7:530251/121‑1 (MQ=255)
                                                                                                                                                                                                     cctcctCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAgcgc                                                                                                                                                                                                        >  7:12915/1‑148 (MQ=255)
                                                                                                                                                                                                       tcctccAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGc                                                                                                                                                                                                             >  6:225425/1‑141 (MQ=255)
                                                                                                                                                                                                       tcctccAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGc                                                                                                                                                                                                             <  5:225425/141‑1 (MQ=255)
                                                                                                                                                                                                       tcctccAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCAGTTCGAATGCACTGAATTCAGATGGCGCAGCAGTTGCCTGGCTAACTGGCTCGTGCACCa                                                                   <  1:221327/279‑1 (MQ=255)
                                                                                                                                                                                                         ctccAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCatta                                                                                                                                                                                                   >  5:36899/1‑149 (MQ=255)
                                                                                                                                                                                                         ctccAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCAGTTCGAATGCACTGAATTCAGATGGCGCAGCAGTTGCCTGGCTAACTGGCTCGTGCAc                                                                     <  2:331198/275‑1 (MQ=255)
                                                                                                                                                                                                              gcgcCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGGTCATCGTGGTGTa                                                                                                                                                                                                                                                                     <  6:873221/78‑1 (MQ=38)
                                                                                                                                                                                                              gcgcCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGGTCATCGTGGTGTa                                                                                                                                                                                                                                                                     >  5:873221/1‑78 (MQ=38)
                                                                                                                                                                                                                       ggACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCAGTTCGAATGCACTGAATTCAGATGGCGCAGCAGTTGCCTGGCTAACTGGCTCGTGCACCATACCGGCACCCa                                                       <  2:557962/274‑1 (MQ=255)
                                                                                                                                                                                                                            gTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATcg                                                                                                                                                                                >  5:965516/1‑149 (MQ=255)
                                                                                                                                                                                                                             ttGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGTCGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATcg                                                                                                                                                                                <  8:668421/148‑1 (MQ=255)
                                                                                                                                                                                                                                     ccGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCAGTTCGAATGCACTGAATTCAGATGGCGCAGCAGTTGCCTGGCTAACTGGCTCGTGCACCATACCGGCACCCACGGTCACATTGCTc                                         >  2:304060/1‑275 (MQ=255)
                                                                                                                                                                                                                                          cccGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGAAAGTGGCGACGAACCAGAGCATTCAATTCCAGTTCGAATGCACTGAATTCAGATGGCGCAGCAGTTGCCTGGCTAACTGGCTCGTGCACCATACCGGCACCCACGGTCACATTGCTCAGGcg                                    >  4:480886/1‑274 (MQ=255)
                                                                                                                                                                                                                                            cGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGc                                                                                                                                                                 >  8:770011/1‑148 (MQ=255)
                                                                                                                                                                                                                                                  ggCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATCTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAg                                                                                                                                                           >  8:348657/1‑148 (MQ=255)
                                                                                                                                                                                                                                                   gCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTg                                                                                                                                                         <  5:824443/149‑1 (MQ=255)
                                                                                                                                                                                                                                                    ccGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTATTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTgg                                                                                                                                                        <  5:128681/149‑1 (MQ=255)
                                                                                                                                                                                                                                                     cGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTgg                                                                                                                                                        <  8:480326/148‑1 (MQ=255)
                                                                                                                                                                                                                                                     cGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGc                                                                                                                                                       <  6:81952/149‑1 (MQ=255)
                                                                                                                                                                                                                                                      ggTAAACCACAGCACTACACCACGATGACCGTGTTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCACCCCAGTGCGGAAAGTGGc                                                                                                                                                       >  7:751973/1‑148 (MQ=255)
                                                                                                                                                                                                                                                           aCCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGaa                                                                                                                                                 >  5:1041076/1‑149 (MQ=255)
                                                                                                                                                                                                                                                            ccACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATTCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGaa                                                                                                                                                 <  8:600642/148‑1 (MQ=255)
                                                                                                                                                                                                                                                            ccACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGaa                                                                                                                                                 <  7:348657/148‑1 (MQ=255)
                                                                                                                                                                                                                                                            ccACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGAGCGCCCCAGTGCGGAAAGTGGCGACGaa                                                                                                                                                 >  8:678681/1‑148 (MQ=255)
                                                                                                                                                                                                                                                             cacaGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAAcc                                                                                                                                               >  5:292688/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                 gCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCagag                                                                                                                                           <  5:798701/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                     tACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCattatt                                                                                                                                                                                                 >  4:559621/1‑91 (MQ=255)
                                                                                                                                                                                                                                                                     tACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTgg                                                                                                                                                        <  6:1036177/132‑1 (MQ=255)
                                                                                                                                                                                                                                                                      acacCACGGTCATCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCAtt                                                                                                                                       >  7:839304/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                          cacGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAAtt                                                                                                                                  >  6:568929/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                            cGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCCGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCAGTTCGAATCCACTGAATTCAGATGGCGCAGCAGTTGCCTGGCTAACTGCCTCGTGCCCCATACCGGCACCCCCGGTCACATTGCCCAGGCGATCGTAAAAAATCACCCCAACCGCACCCGGAttt  >  4:590685/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                            cGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCAGTTCGAATGCACTGAATTCAGATGGCGCAGCAGTTGCCTGGCTAACTGGCTCGTGCACCATACCGGCACCCACGGTCACATTGCTCAGGCGATCGAAAAAAATCAGccc                  >  2:386939/1‑259 (MQ=255)
                                                                                                                                                                                                                                                                              aTGACCGTGGTGTAGCTGGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCAGTTCGAATGCACTGAATTCAGATGGCGCAGCAGTTGCCTGGCTAACt                                                                                <  4:310776/195‑1 (MQ=255)
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                              
CTGGCGTTCGGCGCGGTGTGGCGAGATAAATGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCAGTTCGAATGCACTGAATTCAGATGGCGCAGCAGTTGCCTGGCTAACTGGCTCGTGCACCATACCGGCACCCACGGTCACATTGCTCAGGCGATCGATAAAAATCAGCCCACCCGTCACCGGATTT  >  NC_000913/2873645‑2874187

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