Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 773,446 T→C 5.9% G135G (GGT→GGC cydB → cytochrome d terminal oxidase, subunit II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913773,4460TC 5.9% 245.4 / 2.6 83G135G (GGT→GGCcydBcytochrome d terminal oxidase, subunit II
Reads supporting (aligned to +/- strand):  ref base T (38/39);  new base C (2/3);  total (40/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTATGTGGCGATGATCCTCGTGCTGGCGTCTTTGTTCTTCCGTCCGGTCGGTTTTGACTACCGCTCCAAGATTGAAGAAACCCGCTGGCGTAACATGTGGGACTGGGGCATCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTC  >  NC_000913/773299‑773588
                                                                                                                                                   |                                                                                                                                              
cTATGTGGCGATGATCCTCGTGCTGGCGTCTTTGTTCTTCCGTCCGGTCGGTTTTGACTACCGCTCCAAGATTGAAGAAACCCGCTGGCGTAACATGTGGGACTGGGGCATCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGgt                                                                                                                                                <  4:222396/148‑1 (MQ=255)
cTATGTGGCGATGATCCTCGTGCTGGCGTCTTTGTTCTTCCGTCCGGTCGGTTTTGACTACCGCTCCAAGATTGAAGAAACCCGCTGGCGTAACATGTGGGACTGGGGCATCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGgt                                                                                                                                                <  6:222653/148‑1 (MQ=255)
   tgtgGCGATGATCCTCGTGCTGGCGTCTTTGTTCTTCCGTCCGGTCGGTTTTGACTACCGCTCCAAGATTGAAGAAACCCGCTGGCGTAACATGTGGGACTGGGGCATCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGgtgt                                                                                                                                              <  8:292302/147‑1 (MQ=255)
    gtgGCGATGATCCTCGTGCTGGCGTCTTTGTTCTTCCGTCAGGTCGGTTTTGACTACCGCTCCAAGATTGAATAAAACCGCTGGCGTAAAAAGTGGGACTCGGGCTTATTAATTTGTAGCTTCGTTCCTCCGGTGGTAATTGGTGTAGc                                                                                                                                           >  3:17519/1‑149 (MQ=255)
       gCGATGATCCTCGTGCTGGCGTCTTTGTTCTTCCGTCCGGTCGGTTTTGACTACCGCTCCAAGATTGAAGAAACCCGCTGGCGTAACATGTGGGACTGGGGCATCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGtt                                                                                                                                        <  5:116180/149‑1 (MQ=255)
        cGATGATCCTCGTGCTGGCGTCGTTGTTCTTCCGTCCGGTCGGTTTTGACTACCGCTCCAAGATTGAAGACACCCGCTGGCGTAACATGTGGGACTGGGGCATCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTc                                                                                                                                       <  6:338994/149‑1 (MQ=255)
         gatgatCCTCGTGCTGGCGTCTTTGTTCTTCCGTCCGGTCGGTTTTGACTACCGCTCCAAGATTGAAGAAACCCGCTGGCGTAACATGTGGGACTGGGGCATCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCg                                                                                                                                      >  3:374298/1‑149 (MQ=255)
         gatgatCCTCGTGCTGGCGTCTTTGTTCTTCCGTCCGGTCGGTTTTGACTACCGCTCCAAGATTGAAGAAACCCGCTGGCGTAACATGTGGGACTGGGGCATCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGGGTAGCGTTCg                                                                                                                                      <  6:85061/149‑1 (MQ=255)
               ccTCGTGCTGGCTTCTTTGTTCTTCCTTCCGGTCGGTTTTGACTACCGATCCAAGATTGAAGAAACCCGCTGGAGTAACATGTGGGACTGGGGCATCTTAATTGTTATATTCGTTCCGACGATGGTAATTGGCGTAGCGTTCGGTAAca                                                                                                                                >  7:364221/1‑148 (MQ=255)
               ccTCGTGCTGGCGTCTTTGTTCTTCCGTCCGGTCGGTTTTGACTACCGCTCCAAGATTGAAGAAACCCGCTGGCGTAACATGTGGGACTGGGGCATCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGCGTAGCGTTCGGTAAcc                                                                                                                                >  7:358908/1‑149 (MQ=255)
                 tCGTGCTGGCGTCTTTGTTCTTCCGTCCGGTCGGTTTTGACTACCGCTCCAAGATTGAAGAAACCCGCTGGCGTAACATGTGGGACTGGGGCATCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTg                                                                                                                              >  7:219728/1‑149 (MQ=255)
                   gTGCTGGCGTCTTTGTTCTTCCGTCCGGTCGGTTTTGACTACCGCTCCAAGATTGAAGAAACCCGCTGGCGTAACATGTGGGACTGGGGCATCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGtt                                                                                                                            >  8:351012/1‑149 (MQ=255)
                     gCTGGCGTCTTTGTTCTTCCGTCCGGTCGGTTTTGACTACCGCTCCAAGATTGAAGAAACCCGCTGGCGTAACATGTGGGACTGGGGCATCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTg                                                                                                                           >  4:38210/1‑148 (MQ=255)
                      cTGGCGTCTTTGTTCTTCCGTCCGGTCGGTTTTGACTACCGCTCCAAGATTGAAGAAACCCGCTGGCGTAACATGTGGGACTGGGGCATCTTCATTGGTAGTTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCa                                                                                                                         <  8:225796/149‑1 (MQ=255)
                         gCGTCTTTGTTCTTCCGTCCGGTCGGTTTTGACTACCGCTCCAAGATTGAAGAAACCCGCTGGCGTAACATGTGGGACTGGGGCATCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAggg                                                                                                                      >  4:229115/1‑149 (MQ=255)
                              tttGTTCTTCCGTCCGGTCGGTTTTGACTACCGCTCCAAGATTGAAGAAACCCGCTGGCGTAACATGTGGGACTGGGGCATCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTAc                                                                                                                 >  5:387976/1‑149 (MQ=255)
                                 gTTCTTCCGTCCGGTCGGTTTTGACTACCGCTCCAAGATTGAAGAAACCCGCTGGCGTAACATGTGGGACTGGGGCATCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGt                                                                                                              >  8:350414/1‑149 (MQ=255)
                                   tcttcCGTCCGGTCGGTTTTGACTACCGCTCCAAGATTGAAGAAACCCGCTGGCGTAACATGTGGGACTGGGGCATCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAAc                                                                                                                                 <  8:211377/128‑1 (MQ=255)
                                    cttcCGTCCGGTCGGTTTTGACTACCGCTCCAAGATTGAAGAAACCCGCTGGCGTAACATGTGGGACTGGGGCATCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTc                                                                                                            >  5:128437/1‑148 (MQ=255)
                                     ttcCGTCCGGTCGGTTTTGACTACCGCTCCAAGATTGAAGAAACCCGCTGGCGTAACATGTGGGACTGGGGCATCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCaa                                                                                                          >  8:178980/1‑149 (MQ=255)
                                       ccgtccgGTCGGTTTTGACTACCGCTCCAAGATTGAAGAAACCCGCTGGCGTAACATGTGGGACTGGGGCATCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAAc                                                                                                         <  1:368562/148‑1 (MQ=255)
                                        cgtccgGTCGGTTTTGACTACCGCTCCAAGATTGAAGAAACCCGCTGGCGTAACATGTGGGACTGGGGCATCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGt                                                                                                       >  6:95262/1‑149 (MQ=255)
                                           ccgGTCGGTTTTGACTACCGCTCCAAGATTGAAGAAACCCGCTGGCGTAACATGTGGGACTGGGGCATCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTtga                                                                                                    <  6:387976/149‑1 (MQ=255)
                                              gtcggtTTTGACTACCGCTCCAAGATTGAAGAAACCCGCTGGCGTAACATGTGGGACTGGGGCATCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTtgatg                                                                                                  >  5:121307/1‑148 (MQ=255)
                                                   ttttGACTACCGCTCCAAGATTGAAGAAACCCGCTGGCGTAACATGTGGGAGTGGGGGATCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATc                                                                                            >  5:26953/1‑149 (MQ=255)
                                                   ttttGACTACCGCTCCAAGATTGAAGAAACCCGCTGGCGTAACATGTGGGACTGGGGCATCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATc                                                                                            >  4:142365/1‑149 (MQ=255)
                                                     ttGACTACCGCTCCAAGATTGAAGAAACCCGCTGGCGTAACATGTGGGACTGGGGCATCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCt                                                                                           <  4:181572/148‑1 (MQ=255)
                                                     ttGACTACCGCTCCAAGATTGAAGAAACCCGCTGGCGTAACATGTGGGACTGGGGCATCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCt                                                                                           <  8:100279/148‑1 (MQ=255)
                                                       gACTACCGCTCCAAGATTGAAGAAACCCGCTGGCGTAACATGTGGGACTGGTGCATCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGTGCGTACCGTTCAACGTTGATGAATATCTGCg                                                                                        >  1:111549/1‑149 (MQ=255)
                                                             cGCTCCAAGATTGCAGAAACCCGCTGGCGTAACATGTGGGACTGGGGCATCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGta                                                                                  <  4:287585/149‑1 (MQ=255)
                                                              gCTCCAAGATTCACCCAACCCGCTGGCGGAACATGTGGGACTGGGGCATCTTCACTGGGAGCTTCGCTCCGCCGCTGGTAATTGGCGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGGTGATGAATATCTGCGTCTGtac                                                                                 <  8:383799/149‑1 (MQ=255)
                                                                tCCAAGATTGAAGAAACCCGCTGGCGTAACATGTGGGACTGGGGCATCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGtacta                                                                               >  6:162659/1‑149 (MQ=255)
                                                                         gaagaaACCCGCTGGCGTAACATGTGGGACTGGGGCATCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTa                                                                       <  7:178980/148‑1 (MQ=255)
                                                                          aagaaACCCGCTGGCGTAACATGTGGGACTGGGGCATCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTaa                                                                      <  7:350414/148‑1 (MQ=255)
                                                                          aagaaACCCGCTGGCGTAACATGTGGGACTGGGGCATCTTCAATGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACCCCGGTAAc                                                                     >  2:22603/1‑149 (MQ=255)
                                                                           agaaACCCGCTGGCGTAACATGTGGGACTGTGGCATCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAAct                                                                    >  4:289171/1‑149 (MQ=255)
                                                                           agaaACCCGCTGGCGTAACATGTGGGACTGGGGCATCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGTGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTtga                                                                                                    >  8:289951/1‑117 (MQ=255)
                                                                           agaaACCCGCTGGCGTAACATGTGGGACTGGGGCATCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGTGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTtga                                                                                                    <  7:289951/117‑1 (MQ=255)
                                                                             aaaCCCGCTGGCGTAACATGTGGGACTGGGGCATCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAActtc                                                                  <  5:193797/149‑1 (MQ=255)
                                                                             aaaCCCGCTGGCGTAACATGTGGGACTGGGGCATCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAActtc                                                                  >  6:304539/1‑149 (MQ=255)
                                                                              accccGCTGGCGTAACATGTGGGACTGGGGCATCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAActtct                                                                 >  7:255286/3‑149 (MQ=255)
                                                                               aCCCGCTGGCGTAACATGTGGGACTGGGGCATCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAActtctt                                                                >  8:273558/1‑149 (MQ=255)
                                                                                  cGCTGGCGTAACATGTGGGACTGGGGCATCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACttcttc                                                               >  5:180353/1‑147 (MQ=255)
                                                                                     tGGCGTAACATGTGGGACTGGGGCATCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGtt                                                          <  3:38210/149‑1 (MQ=255)
                                                                                       gCGTAACATGTGGGACTGGGGCATCTTCATTGCTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGTCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGc                                                        <  6:18760/149‑1 (MQ=255)
                                                                                         gTAACATGTGGGACTGGGGCATCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAActtct                                                                 >  2:5988/1‑138 (MQ=255)
                                                                                         gTAACATGTGGGACTGGGGCATCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAActtct                                                                 <  1:5988/138‑1 (MQ=255)
                                                                                         gTAACATGTGGGACTGGGGCATCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCtt                                                      <  5:304539/149‑1 (MQ=255)
                                                                                         gTAACATGTGGGACTGGGGCATCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCt                                                       >  5:291045/1‑148 (MQ=255)
                                                                                            acttgtgtgACTGTGGCATCTTCATTGGTAGCTTCGTTCCGCAGCTGGTAATTGGCGTAGCTTTCTGTAACCTTTTGCAGGGAGTACCGTTAAACGTTGATGAATATCTGCGTATGTAATACACAGGTAACTTCTTCCATTTGCTTAAc                                                   <  8:364221/146‑1 (MQ=255)
                                                                                            aCATGTGGGACTGGGGCATCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAAc                                                   <  7:351012/149‑1 (MQ=255)
                                                                                            aCATGTGGGACTGGGGCATCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGCGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTTATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAAc                                                   <  8:358908/149‑1 (MQ=255)
                                                                                              aTGTGGGACTGGGGCATCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTaa                                                                      >  4:136003/1‑128 (MQ=255)
                                                                                              aTGTGGGACTGGGGCATCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTaa                                                                      <  3:136003/128‑1 (MQ=255)
                                                                                               tgtgGGACTGGGGCATCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCg                                                >  4:371422/1‑149 (MQ=255)
                                                                                                  gggACTGGGGCATCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTc                                             >  7:78693/1‑149 (MQ=255)
                                                                                                   ggACTTGGTGATATTCATTTGTATATTCTTTCCGCCGCTTGTAATTTGTGTAGAGTTCGGTAACCTTTTGCAGGGCGTACAGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCg                                            <  4:17519/149‑1 (MQ=255)
                                                                                                    gACTGGGGCATCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCgg                                           <  3:142365/149‑1 (MQ=255)
                                                                                                     aCTGGGGCATCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTGCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGc                                          <  1:337326/149‑1 (MQ=255)
                                                                                                         gggCATCTTCATTGGGAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGGCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCctgc                                      <  5:129394/149‑1 (MQ=255)
                                                                                                            cATCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTgg                                   >  5:188381/1‑149 (MQ=255)
                                                                                                              tCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTgg                                   >  8:216952/1‑147 (MQ=255)
                                                                                                              tCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTgg                                   >  4:246596/1‑147 (MQ=255)
                                                                                                                ttCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACttcttc                                                               >  8:13051/1‑117 (MQ=255)
                                                                                                                ttCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACttcttc                                                               <  7:13051/117‑1 (MQ=255)
                                                                                                                 tCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTaca                                                                             <  2:142479/102‑1 (MQ=255)
                                                                                                                 tCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTaca                                                                             >  1:142479/1‑102 (MQ=255)
                                                                                                                  cATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCa                                 <  6:26953/145‑1 (MQ=255)
                                                                                                                   aTTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCgt                            <  6:188381/149‑1 (MQ=255)
                                                                                                                   aTTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCgt                            <  2:111549/149‑1 (MQ=255)
                                                                                                                   aTTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCgt                            <  3:49551/149‑1 (MQ=255)
                                                                                                                    ttGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCgtg                           >  2:336364/1‑149 (MQ=255)
                                                                                                                    ttGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCgtg                           <  8:250181/149‑1 (MQ=255)
                                                                                                                      ggTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCgtg                           <  5:189120/147‑1 (MQ=255)
                                                                                                                      ggTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCgtg                           >  6:189120/1‑147 (MQ=255)
                                                                                                                       gTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCTTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCCAGCAGGccgaac                         <  6:180353/148‑6 (MQ=255)
                                                                                                                       gTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCCAGCAGGccgaac                         <  3:177287/148‑6 (MQ=255)
                                                                                                                               gTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCgtggt                         <  4:239584/140‑1 (MQ=255)
                                                                                                                               gTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCgtggt                         >  3:239584/1‑140 (MQ=255)
                                                                                                                                ttCCGCCGCTGGTAATTGGTGTAGCGTTCGGTACCCTGGTGCAGGGCGTACCGGTCAACGTTGATGACGATCTGCGTCTGTCCTACACCGGGACCTTCTTCCAGTGGCTTAACCCGTTCGGACGGCTGGCAGGCGTGGTGAGCGTAgt                >  6:32118/1‑147 (MQ=255)
                                                                                                                                  ccgccgCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTATGga              <  6:121307/148‑1 (MQ=255)
                                                                                                                                  ccgccgCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGGCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGga              <  3:289171/148‑1 (MQ=255)
                                                                                                                                    gccgcTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGgatga           <  6:291045/149‑1 (MQ=255)
                                                                                                                                       gcTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCa        <  8:219728/149‑1 (MQ=255)
                                                                                                                                           gTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTAc    >  7:57306/1‑149 (MQ=255)
                                                                                                                                              aTTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTc  <  3:229115/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
CTATGTGGCGATGATCCTCGTGCTGGCGTCTTTGTTCTTCCGTCCGGTCGGTTTTGACTACCGCTCCAAGATTGAAGAAACCCGCTGGCGTAACATGTGGGACTGGGGCATCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTC  >  NC_000913/773299‑773588

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