Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,078,131 C→A 79.8% A228E (GCG→GAG)  yiiD → GNAT family putative N‑acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,078,1310CA79.8% 202.0 / 64.3 94A228E (GCG→GAG) yiiDGNAT family putative N‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (11/8);  new base A (39/36);  total (50/44)
Fisher's exact test for biased strand distribution p-value = 7.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.29e-01

GGTACGAACATATCCCGCTTAGTGAAAAAATGGGCGTGCGCATTCAGCAATATACCGGGCAAAAATTTATCACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGCGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGCCGATTAGCGGTAAACCTCATGCGGTAGCCGACCTCGGTGCCTTAAG  >  NC_000913/4077990‑4078275
                                                                                                                                             |                                                                                                                                                
ggTACGAACATATCCCGCTTAGTGAAAAAATGGGCGTGCGCATTCAGCAATATACCGGGCAAAAATTTATCACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGAGACGCTc                                                                                                                                           <  1:392112/149‑1 (MQ=255)
 gTACGAACATATCCCGCTTAGTGAAAAAATGGGCGTGCGCATTCAGCAATATACCGGGCAAAAATTTATCACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGAGACGCTCa                                                                                                                                          >  2:733027/1‑149 (MQ=255)
  tACGAACATATCCCGCTTAGTGAAAAAATGGGCGTGCGCATTCAGCAATATACCGGGCAAAAATTTATCACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGCGACGCTCAc                                                                                                                                         >  1:58156/1‑149 (MQ=255)
     gAACATATCCCGCTTAGTGAAAAAATGGGCGTGCGCATTCAGCAATATACCGGGCAAAAATTTATCACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGAGACGCTCACTgg                                                                                                                                      >  2:516095/1‑149 (MQ=255)
     gAACATATCCCGCTTAGTGAAAAAATGGGCGTGCGCATTCAGCAATATACCGGGCAAAAATTTATCACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGAGACGCTCACTgg                                                                                                                                      >  4:1338/1‑149 (MQ=255)
      aaCATATCCCGCTTAGTGAAAAAATGGGCGTGCGCATTCAGCAATATACCGGGCAAAAATTTATCACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGAGACGCTCACTGGt                                                                                                                                     >  2:92316/1‑149 (MQ=255)
           atCCCGCTTAGTGAAAAAATGGGCGTGCGCATTCAGCAATATACCGGGCAAAAATTTATCACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGAGACGCTCACTGGTTgggg                                                                                                                                >  2:331326/1‑149 (MQ=255)
             cccGCTTAGTGAAAAAATGGGCGTGCGCATTCAGCAATATACCGGGCAAAAATTTATCACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGAGACGCTCACTGGTTGGGGAc                                                                                                                              >  1:270244/1‑149 (MQ=255)
              ccGCTTAGTGAAAAAATGGGCGTGCGCATTCAGCAATATACCGGGCAAAAATTTATCACTACCATGCCAGAAACCGGCAATTAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGAGACGCTCACTGGTTGGGGACt                                                                                                                             >  3:312913/1‑149 (MQ=255)
              ccGCTTAGTGAAAAAATGGGCGTGCGCATTCAGCAATATACCGGGCAAAAATTTATCACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGAGACGCTCACTGGTTGGGGACt                                                                                                                             >  3:276013/1‑149 (MQ=255)
              ccGCTTAGTGAAAAAATGGGCGTGCGCATTCAGCAATATACCGGGCAAAAATTTATCACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGAGACGCTCACTGGTTGGGGACt                                                                                                                             >  3:823987/1‑149 (MQ=255)
               cGCTTAGTGAAAAAATGGGCGTGCGCATTCAGCAATATACCGGGCAAAAATTTATCACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGAGACGCTCACTGGTTGGGGACtt                                                                                                                            >  1:861713/1‑149 (MQ=255)
                  ttAGCGAAAGAATGGGCGTGCGGATACAGCAATATACCGGGCAAAAACTTATCGCTACCATGACAGAAACCGGCAATCAGAGTCCGCACCATACGCTGTTTCCCGGGAGTTTATTCTCACTGGAGACGCTCACTGGTTGGGGACTTATc                                                                                                                         <  2:549567/149‑1 (MQ=255)
                     gTGAAAAAATGGGCGTGCGCATTCAGCAATATACCGGGCAAAAATTTATCACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGAGACGCTCACTGGTTGGGGACTTATCTgg                                                                                                                      >  1:364567/1‑149 (MQ=255)
                     gTGAAAAAATGGGCGTGCGCATTCAGCAATATACCGGGCAAAAATTTATCACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGAGACGCTCACTGGTTGGGGACTTATCTgg                                                                                                                      <  1:331326/149‑1 (MQ=255)
                     gTGAAAAAATGGGCGTGCGCATTCAGCAATATACCGGGCAAAAATTTATCACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGAGACGCTCACTGGTTGGGGACTTATCTgg                                                                                                                      <  2:697315/149‑1 (MQ=255)
                      tGAAAAAATGGGCGTGCGCATTCAGCAATATACCGGGCAAAAATTTATCACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGAGACGCTCACTGGTTGGGGACTTATCTGGc                                                                                                                     <  1:733027/149‑1 (MQ=255)
                        aaaaaaTGGGCGTGCGCATTCAGCAATATACCGGGCAAAAATTTATCACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGCGACGCTCACTGGTTGGGGACTTATCTGGCTg                                                                                                                   >  2:879640/1‑149 (MQ=255)
                        aaaaaaTGGGCGTGCGCATTCAGCAATATACCGGGCAAAAATTTATCACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGAGACACTCACTGGTTGGGGACTTATCTGGCTg                                                                                                                   >  3:636245/1‑149 (MQ=255)
                          aaaaTGGGCGTGCGCATTCAGCAATATACCGGGCAAAAATTTATCACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGAGACGCTCACTGGTTGGGGACTTATCTGGCTGAt                                                                                                                 >  4:126231/1‑149 (MQ=255)
                             aTGGGCGTGCGCATTCAGCAATATACCGGGCAAAAATTTATCACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGAGACGCTCACTGGTTGGGGACTTATCTGGCTGAtgct                                                                                                              <  4:325975/149‑1 (MQ=255)
                             aTGGGCGTGCGCATTCAGCAATATACCGGGCAAAAATTTATCACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGAGACGCTCACTGGTTGGGGACTTATCTGGCTGAtgct                                                                                                              <  3:905373/149‑1 (MQ=255)
                              tGGGCGTGCGCATTCAGCAATATACCGGGCAAAAATTTATCACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGCGACGCTCACTGGTTGGGGACTTATCTGGCTGAtgctg                                                                                                             >  2:742855/1‑149 (MQ=255)
                               gggCGTGCGCATTCAGCAATATACCGGGCAAAAATTTATCACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGAGACGCTCACTGGTTGGGGACTTATCTGGCTGAtgctgc                                                                                                            <  2:93713/149‑1 (MQ=255)
                               gggCGTGCGCATTCAGCAATATACCGGGCAAAAATTTATCACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGAGACGCTCACTGGTTGGGGACTTATCTGGCTGAtgctgc                                                                                                            <  2:541064/149‑1 (MQ=255)
                                ggCGTGCGCATTCAGCAATATACCGGGCAAAAATTTATCACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGAGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCg                                                                                                           <  3:926300/149‑1 (MQ=255)
                                    tGCGCATTCAGCAATATACCGGGCAAAAATTTATCACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGAGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGaa                                                                                                       <  2:270244/149‑1 (MQ=255)
                                       gcATTCAGCAATATACCGGGCAAAAATTTATCACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGCGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGc                                                                                                    >  4:766347/1‑149 (MQ=255)
                                       gcATTCAGCAATATACCGGGCAAAAATTTATCACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGAGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGc                                                                                                    <  2:5848/149‑1 (MQ=255)
                                         aTTCAGCAATATACCGGGCAAAAATTTATCACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGAGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCa                                                                                                  >  4:141936/1‑149 (MQ=255)
                                         aTTCAGCAATATACCGGGCAAAAATTTATCACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGAGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCa                                                                                                  >  3:771203/1‑149 (MQ=255)
                                          ttCAGCAATATACCGGGCAAAAATTTATCACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGAGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCAc                                                                                                 >  1:176002/1‑149 (MQ=255)
                                          ttCAGCAATATACCGGGCAAAAATTTATCACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGAGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCAc                                                                                                 >  1:410812/1‑149 (MQ=255)
                                           tCAGCAATATACCGGGCAAAAATTTATCACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGCGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCAcc                                                                                                <  1:679580/149‑1 (MQ=255)
                                                 atatACCGGGCAAAAATTTATCACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGTGTTTATTCTCACTGGAGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTcggcg                                                                                          <  2:151151/149‑1 (MQ=255)
                                                 atatACCGGGCAAAAATTTATCACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGAGACGCTCACTGGTTGGGGACTTATCAGGCTGATGCTGCGTGAACGCCACCTcggcg                                                                                          >  3:753076/1‑149 (MQ=255)
                                                      ccGGGCAAAAATTTATCACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGAGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACg                                                                                     >  1:823438/1‑149 (MQ=255)
                                                       cGGGCAAAAATTTATCACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGAGACGCTCACTGGTTGGGGACTTATCTGGCTGAtgctgc                                                                                                            >  4:326152/1‑125 (MQ=255)
                                                       cGGGCAAAAATTTATCACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGAGACGCTCACTGGTTGGGGACTTATCTGGCTGAtgctgc                                                                                                            <  3:326152/125‑1 (MQ=255)
                                                         ggCAAAAATTTATCACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGAGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGatt                                                                                  <  2:823438/149‑1 (MQ=255)
                                                          gCAAAAATTTATCACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGAGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGatta                                                                                 >  3:456872/1‑149 (MQ=255)
                                                          gCAAAAATTTATCACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGAGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGatta                                                                                 >  4:31612/1‑149 (MQ=255)
                                                           cAAAAATTTATCACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGAGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGattat                                                                                <  1:849997/149‑1 (MQ=255)
                                                             aaaaTTTATCACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGAGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTc                                                                              >  1:74460/1‑149 (MQ=255)
                                                               aaTTTATCACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGAGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCtt                                                                            >  3:461960/1‑149 (MQ=255)
                                                               aaTTTATCACTACCATACCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGAGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCtt                                                                            <  4:753076/149‑1 (MQ=255)
                                                                aTTTATCACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGCGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTg                                                                           >  1:832347/1‑149 (MQ=255)
                                                                aTTTATCACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGAGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTg                                                                           <  1:100532/149‑1 (MQ=255)
                                                                 tttATCACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGCGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGc                                                                          <  2:719528/149‑1 (MQ=255)
                                                                 tttATCACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGAGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGc                                                                          >  1:179976/1‑149 (MQ=255)
                                                                  ttATCACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGAGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCg                                                                         >  1:878028/1‑149 (MQ=255)
                                                                        cTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGAGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGaa                                                                                       <  4:394163/129‑1 (MQ=255)
                                                                        cTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGAGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGaa                                                                                       >  3:394163/1‑129 (MQ=255)
                                                                         tACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGCGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATgcgc                                                                  >  4:765893/1‑149 (MQ=255)
                                                                              tGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGAGACGCTCACTGGTTGGGGACTTATCTgg                                                                                                                      <  1:207855/92‑1 (MQ=255)
                                                                              tGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGAGACGCTCACTGGTTGGGGACTTATCTgg                                                                                                                      >  2:207855/1‑92 (MQ=255)
                                                                                 cAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGCGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGc                                                          <  3:765893/149‑1 (MQ=255)
                                                                                 cAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGAGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGc                                                          >  4:445732/1‑149 (MQ=255)
                                                                                 cAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGAGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGc                                                          >  3:687517/1‑149 (MQ=255)
                                                                                 cAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGAGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGc                                                          <  4:745585/149‑1 (MQ=255)
                                                                                   gAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGAGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTa                                                        <  4:276013/149‑1 (MQ=255)
                                                                                    aaaCCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGAGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTAc                                                       >  3:583338/1‑149 (MQ=255)
                                                                                     aaCCGGCAATCAGAATCCGCACCATACACTGTTTGCCGGGAGTTTATTCTCACTGGAGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACa                                                      >  4:101408/1‑149 (MQ=255)
                                                                                        cGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGAGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCa                                                   >  2:460648/1‑149 (MQ=255)
                                                                                           cAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGCGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGc                                                >  1:129681/1‑149 (MQ=255)
                                                                                            aaTCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGCGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGcc                                               >  3:408297/1‑149 (MQ=255)
                                                                                              tCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGAGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGCCGa                                             >  1:589676/1‑149 (MQ=255)
                                                                                              tCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGAGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGCCGa                                             <  4:583338/149‑1 (MQ=255)
                                                                                               cAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGAGACGCTCACTGGTAGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGCCGAt                                            <  4:636245/149‑1 (MQ=255)
                                                                                                aGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGCGACGCTCAc                                                                                                                                         <  4:497468/55‑1 (MQ=255)
                                                                                                aGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGCGACGCTCAc                                                                                                                                         >  3:497468/1‑55 (MQ=255)
                                                                                                aGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGAGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGCCGAtt                                           <  1:447143/149‑1 (MQ=255)
                                                                                                  aaTCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGAGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGCCGATTAg                                         >  3:706483/1‑149 (MQ=255)
                                                                                                   aTCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGAGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGCCGATTAGc                                        <  2:364567/149‑1 (MQ=255)
                                                                                                    tCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGAGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGCCGATTAGCg                                       >  4:14927/1‑149 (MQ=255)
                                                                                                         aCCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGAGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGCCGATTAGCGGTaaa                                  <  3:445732/149‑1 (MQ=255)
                                                                                                          ccATACGCTGTTTGCCGGGAGTTTATTCTCACTGGCGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTc                                                                                              >  4:468015/1‑88 (MQ=255)
                                                                                                          ccATACGCTGTTTGCCGGGAGTTTATTCTCACTGGCGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTc                                                                                              <  3:468015/88‑1 (MQ=255)
                                                                                                             tACGCTGTTTGCCGGGAGTTTATTCTCACTGGAGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGCCGATTAGCGGTAAACCTc                              >  3:331347/1‑149 (MQ=255)
                                                                                                               cGCTGTTTGCCGGGAGTTTATTCTCACTGGAGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGCCGATTAGCGGTAAACCCCAt                            <  2:540939/149‑1 (MQ=255)
                                                                                                                        ccGGGAGTTTATTCTCACTGGCGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGCCGATTAGCGGTAAACCTCATGCGGTAGcc                   <  2:58156/149‑1 (MQ=255)
                                                                                                                        ccGGGAGTTTATTCTCACTGGAGACGCTCACTGGGTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGCCGATTAGCGGTAAACCTCATGCGGTAGcc                   <  2:461797/149‑1 (MQ=255)
                                                                                                                           ggAGTTTATTCTCACTGGCGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGCCGATTAGCGGTAAACCTCATGCGGTAGCCGAc                >  3:325032/1‑149 (MQ=255)
                                                                                                                           ggAGTTTATTCTCACTGGAGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGCCGATTAGCGGTAAACCTCATGCGGTAGCCGAc                >  2:682507/1‑149 (MQ=255)
                                                                                                                             aGTTTATTCTCACTGGAGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGCCGATTAGCGGTAAACCTCATGCGGTAGCCGACCt              <  1:516095/149‑1 (MQ=255)
                                                                                                                               tttATTCTCACTGGCGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGCCGATTAGCGGTAAACCTCATGCGGTAGCCGACCTCg            <  4:408297/149‑1 (MQ=255)
                                                                                                                                 tATTCTCACTGGAGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGCCGATTAGCGGTAAACCTCATGCGGTAGCCGACCTCGGt          <  3:882781/149‑1 (MQ=255)
                                                                                                                                 tATTCTCACTGGAGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGCCGATTAGCGGTAAACCTCATGCGGTAGCCGACCTCGGt          <  3:730622/149‑1 (MQ=255)
                                                                                                                                 tATTCTCACTGGAGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGCCGATTAGCGGTAAACCTCATGCGGTAGCCGACCTCGGt          <  3:65998/149‑1 (MQ=255)
                                                                                                                                  aTTCTCACTGGAGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGCCGATTAGCGGTAAACCTCATGCGGTAGCCGACCTCGGTg         <  4:461960/149‑1 (MQ=255)
                                                                                                                                   ttCTCACTGGCGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGCCGATTAGCGGTAAACCTCATGCGGTAGCCGACCTCGGTGc        <  2:570660/149‑1 (MQ=255)
                                                                                                                                    tctcACTGGAGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGCCGATTAGCgg                                      <  4:446243/118‑1 (MQ=255)
                                                                                                                                    tctcACTGGAGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGCCGATTAGCgg                                      >  3:446243/1‑118 (MQ=255)
                                                                                                                                         cTGGAGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGCCGATTAGCGGTAAACCTCATGCGGTAGCCGACCTCGGTGCCTTAAg  <  2:229809/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                
GGTACGAACATATCCCGCTTAGTGAAAAAATGGGCGTGCGCATTCAGCAATATACCGGGCAAAAATTTATCACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGCGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGCCGATTAGCGGTAAACCTCATGCGGTAGCCGACCTCGGTGCCTTAAG  >  NC_000913/4077990‑4078275

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