Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,801,372 C→T 7.1% C400Y (TGT→TAT) ‡ thrS ← threonyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,801,3720CT 7.1% 364.7 / 14.3 113C400Y (TGT→TAT) ‡thrSthreonyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base C (54/51);  new base T (4/4);  total (58/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.08e-01

Only 100 of 113 total aligned reads displayed.

TCCAGCGCAACCGCCAGGTCCGCCTCAGCACGATCCCACATTTCGTCGCTGCCAATACGTTTTTCAGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCAGGCCATGCAGCGAACCTGACGGCTCGTTACGGTGGCAGCTACCAAACTCGGCCATACGC  >  NC_000913/1801234‑1801512
                                                                                                                                          |                                                                                                                                            
tcCAGCGCAACCGCCAGGTCCGCCTCAGCACGATCCCACATTTCGTCGCTGCCAATACGTTTTTCAGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTaa                                                                                                                                     >  4:20787/1‑148 (MQ=255)
    gcgcAACCGCCAGGTCCGCCTCAGCACGATCCCACATTTCGTCGCTGCCAATACGTTTTTCAGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTc                                                                                                                                 >  1:612339/1‑148 (MQ=255)
            gCCAGGTCCGCCTCAGCACGATCCCACATTTCGTCGCTGCCAATACGTTTTTCAGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGaa                                                                                                                         <  2:612339/148‑1 (MQ=255)
             ccAGGTCCGCCTCAGCACGATCCCACATTTCGTCGCTGCCAATACGTTTTTCAGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGt                                                                                                                                        <  4:296998/132‑1 (MQ=255)
             ccAGGTCCGCCTCAGCACGATCCCACATTTCGTCGCTGCCAATACGTTTTTCAGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGt                                                                                                                                        >  3:296998/1‑132 (MQ=255)
                  tCCGCCTCAGCACGATCCCACATTTCGTCGCTGCCAATACGTTTTTCAGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGtt                                                                                                                   <  1:1105580/148‑1 (MQ=255)
                        tCAGCACGATCCCACATTTCGTCGCTGCCAATACGTTTTTCAGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAg                                                                                                             >  1:343849/1‑148 (MQ=255)
                                aTCCCACATTTCGTCGCTGCCAATACGTTTTTCAGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGtt                                                                                                                   >  3:391440/1‑134 (MQ=255)
                                aTCCCACATTTCGTCGCTGCCAATACGTTTTTCAGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGtt                                                                                                                   <  4:391440/134‑1 (MQ=255)
                                aTCCCACATTTCGTCGCTGCCAATACGTTTTTCAGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGt                                                                                                                    <  3:448953/133‑1 (MQ=255)
                                aTCCCACATTTCGTCGCTGCCAATACGTTTTTCAGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGt                                                                                                                    >  4:448953/1‑133 (MQ=255)
                                aTCCCACATTTCGTCGCTGCCAATACGTTTTTCAGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACagaag                                                                                                     >  3:1130936/1‑148 (MQ=255)
                                 tCCCACATTTCGTCGCTGCCAATACGTTTTTCAGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATcg                                                                                                                             >  1:598873/1‑123 (MQ=255)
                                 tCCCACATTTCGTCGCTGCCAATACGTTTTTCAGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATcg                                                                                                                             <  2:598873/123‑1 (MQ=255)
                                 tCCCACATTTCGTCGCTGCCAATACGTTTTTCAGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACaga                                                                                                       <  1:610429/145‑1 (MQ=255)
                                 tCCCACATTTCGTCGCTGCCAATACGTTTTTCAGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACaga                                                                                                       >  2:610429/1‑145 (MQ=255)
                                    cacaTTTCGTCGCTGCCAATACGTTTTTCAGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACa                                                                                                         <  1:68904/140‑1 (MQ=255)
                                      caTTTCGTCGCTGCCAATACGTTTTTCAGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGtt                                                                                                                   >  4:264090/1‑128 (MQ=255)
                                      caTTTCGTCGCTGCCAATACGTTTTTCAGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGtt                                                                                                                   <  3:264090/128‑1 (MQ=255)
                                        tttCGTCGCTGCCAATACGTTTTTCAGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCg                                                                                                                                         >  4:110400/1‑104 (MQ=255)
                                        tttCGTCGCTGCCAATACGTTTTTCAGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCg                                                                                                                                         <  3:110400/104‑1 (MQ=255)
                                               gCTGCCAATACGTTTTTCAGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATcc                                                                                      >  4:1126341/1‑148 (MQ=255)
                                                 tGCCAATACGTTTTTCAGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTg                                                                                    <  4:1130936/148‑1 (MQ=255)
                                                      aTACGTTTTTCAGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAATTGCTGTACATATCATAGACTAAACGGATACATCCg                                                                                                                                         <  4:964257/90‑1 (MQ=255)
                                                      aTACGTTTTTCAGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGt                                                                                                                    <  3:774923/111‑1 (MQ=255)
                                                      aTACGTTTTTCAGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGt                                                                                                                    >  4:774923/1‑111 (MQ=255)
                                                      aTACGTTTTTCAGGACGAGTGGACAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCg                                                                                                                                         >  3:964257/1‑90 (MQ=255)
                                                       tACGTTTTTCAGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATgcg                                                                                              <  3:508327/132‑1 (MQ=255)
                                                           tttttCAGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAGAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCAt                                                                                                                               >  1:731541/1‑95 (MQ=255)
                                                           tttttCAGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATGTATCCGttta                                                                                                                                     >  2:257311/1‑87 (MQ=255)
                                                           tttttCAGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCAt                                                                                                                               <  2:731541/95‑1 (MQ=255)
                                                                 aGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGt                                                                                                                    >  3:933038/1‑100 (MQ=255)
                                                                 aGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACACATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGt                                                                                                                    <  4:933038/100‑1 (MQ=255)
                                                                   gACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCAcgcg                                                                  >  1:40850/1‑148 (MQ=255)
                                                                    aCGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACaga                                                                                                       >  3:156799/1‑110 (MQ=255)
                                                                    aCGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACaga                                                                                                       <  4:156799/110‑1 (MQ=255)
                                                                        gTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCa                                                             <  3:1126341/148‑1 (MQ=255)
                                                                        gTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCa                                                             <  2:40850/148‑1 (MQ=255)
                                                                         tGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCAg                                                            <  2:343849/148‑1 (MQ=255)
                                                                           gagaGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCa                                                             <  2:169269/145‑1 (MQ=255)
                                                                           gagaGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCa                                                             >  1:169269/1‑145 (MQ=255)
                                                                             gagTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCAGGCCa                                                        >  1:984780/1‑148 (MQ=255)
                                                                              agTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATg                                                                                                >  1:166636/1‑107 (MQ=255)
                                                                              agTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATg                                                                                                <  2:166636/107‑1 (MQ=255)
                                                                              agTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCAt                                                                                        <  2:515252/115‑1 (MQ=255)
                                                                              agTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCAGGCCAt                                                       >  2:357470/1‑148 (MQ=255)
                                                                                     cgacgaTCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCAGGCCATGCAGCGa                                                <  1:357470/148‑1 (MQ=255)
                                                                                       acgaTCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCAGGCCATGCAGCGAAc                                              >  3:1039285/1‑148 (MQ=255)
                                                                                          aTCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTa                                                                                                                                      <  1:385137/57‑1 (MQ=255)
                                                                                            cttctCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAAt                                                                                                                        >  4:418026/1‑69 (MQ=255)
                                                                                            cttctCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAAt                                                                                                                        <  3:418026/69‑1 (MQ=255)
                                                                                            cttctCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTg                                                                                    <  4:613333/105‑1 (MQ=255)
                                                                                              tctcGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAAtt                                                                                                                       <  2:809556/68‑1 (MQ=255)
                                                                                              tctcGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAAtt                                                                                                                       >  1:809556/1‑68 (MQ=255)
                                                                                              tctcGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGt                                                                                                            <  1:213594/79‑1 (MQ=255)
                                                                                              tctcGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGt                                                                                                            >  2:213594/1‑79 (MQ=255)
                                                                                                   aaGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCAGGCCATGCAGCGAACCTGACGGCTCGt                                  >  2:858792/1‑148 (MQ=255)
                                                                                                       cAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGCTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCAGGCCATGCAGCGAACCTGACGGCTCGTTACg                              >  3:699961/1‑148 (MQ=255)
                                                                                                         aaaGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCAGGCCATGCAGCGAACCTGACGGCTCGTTACGGt                            >  4:868031/1‑148 (MQ=255)
                                                                                                                 gTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCAGGCCATGCAGCGAACCTGACGGCTCGTTACGGTGGCAGCTa                    <  1:858792/148‑1 (MQ=255)
                                                                                                                   aCATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCAGGc                                                          <  2:747358/108‑1 (MQ=255)
                                                                                                                   aCATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCAGGCCATGCAGCGAACCTGACGGCTCGTTACGGTGGCAGCTAcc                  >  1:1084198/1‑148 (MQ=255)
                                                                                                                   aCATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCAGGCCATGCAGCGAACCTGACGGCTCGTTACGGTGGCAGCTAcc                  >  2:943968/1‑148 (MQ=255)
                                                                                                                    cATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCAGGCCATGCAGCGAACCTGACGGCt                                     <  1:876490/128‑1 (MQ=255)
                                                                                                                    cATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCAGGCCATGCAGCGAACCTGACGGCt                                     >  2:876490/1‑128 (MQ=255)
                                                                                                                        tCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCAt                                                                                        <  3:185402/73‑1 (MQ=255)
                                                                                                                         cATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATc                                                              >  4:737200/1‑98 (MQ=255)
                                                                                                                         cATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATc                                                              <  3:737200/98‑1 (MQ=255)
                                                                                                                         cATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCAGGCCATGCAGCGa                                                >  4:643590/1‑112 (MQ=255)
                                                                                                                         cATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCAGGCCATGCAGCGa                                                <  3:643590/112‑1 (MQ=255)
                                                                                                                         cATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCAGGCCATGCAGCGAACCTGa                                          >  1:187016/1‑118 (MQ=255)
                                                                                                                          aTAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCca                                                                          <  1:582435/85‑1 (MQ=255)
                                                                                                                          aTAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCca                                                                          >  2:582435/1‑85 (MQ=255)
                                                                                                                          aTAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCca                                                                          <  1:657368/85‑1 (MQ=255)
                                                                                                                          aTAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCca                                                                          >  2:657368/1‑85 (MQ=255)
                                                                                                                          aTAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCAGGCCATGCAGCGAACCTGACGGCTCGTTACGGTGGCAGCTACCAAACTCg           <  3:868031/148‑1 (MQ=255)
                                                                                                                           tAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTgg                                                                                   >  2:292854/1‑75 (MQ=255)
                                                                                                                           tAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTgg                                                                                   <  1:292854/75‑1 (MQ=255)
                                                                                                                           tAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCAGGCCATGCAGCGAACCTGACGGCTCGTTACGGTGGCAGCTACCAAACTCgg          <  3:338362/148‑1 (MQ=255)
                                                                                                                           tAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCAGGCCATGCAGCGAACCTGACGGCTCGTTACGGTGGCAGCTACCAAACTCgg          >  4:338362/1‑148 (MQ=255)
                                                                                                                            aGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCAGGCCATGCAGCGAACCTGACGGCTCGTTACGGTGGCAGCTACCAAACTCGGc         <  4:1039285/148‑1 (MQ=255)
                                                                                                                                ttaaCGGATGTATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTg                                                                                    >  1:23669/3‑69 (MQ=255)
                                                                                                                                ttaaCGGATGTATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTg                                                                                    <  2:23669/67‑1 (MQ=255)
                                                                                                                                 aaaCGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATg                                                                                                >  3:1038742/1‑56 (MQ=255)
                                                                                                                                 aaaCGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATg                                                                                                <  4:1038742/56‑1 (MQ=255)
                                                                                                                                 aaaCGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCAGGCCATGCAGCGAACCTGACGGCTCGTTACGGTGGCAGCTACCAAACTCgg          <  1:382596/142‑1 (MQ=255)
                                                                                                                                 aaaCGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCAGGCCATGCAGCGAACCTGACGGCTCGTTACGGTGGCAGCTACCAAACTCgg          >  2:382596/1‑142 (MQ=255)
                                                                                                                                  aaCGGATGTATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCAGGCCATGCAGCGa                                                >  1:292518/1‑103 (MQ=255)
                                                                                                                                  aaCGGATGTATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCAGGCCATGCAGCGa                                                <  2:292518/103‑1 (MQ=255)
                                                                                                                                  aaCGGATGTATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCAGGCCATGCAGCGAACCTGACGGCTCGTTACGGTGGc                         >  1:132287/1‑126 (MQ=255)
                                                                                                                                  aaCGGATGTATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCAGGCCATGCAGCGAACCTGACGGCTCGTTACGGTGGc                         <  2:132287/126‑1 (MQ=255)
                                                                                                                                   aCGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCAGGCCATGCAGCGAACCTGACGGCTCGTTACGGTGGCAGCTACCAAACTCGGCCATACGc  >  4:1127244/1‑148 (MQ=255)
                                                                                                                                    cGGATACATCCGTTAACTTCATCGCGAATTTGttctt                                                                                                                >  1:742260/1‑37 (MQ=255)
                                                                                                                                    cGGATACATCCGTTAACTTCATCGCGAATTTGttctt                                                                                                                <  2:742260/37‑1 (MQ=255)
                                                                                                                                    cGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCAgg                                                           >  4:136187/1‑90 (MQ=255)
                                                                                                                                    cGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCAgg                                                           <  3:136187/90‑1 (MQ=255)
                                                                                                                                        tACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATcc                                                                                      <  4:457137/59‑1 (MQ=255)
                                                                                                                                        tACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATcc                                                                                      >  3:457137/1‑59 (MQ=255)
                                                                                                                                         aCATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCAGGCCATGCAGCGAACCTGACGGCTCGTTa                                <  2:399946/112‑1 (MQ=255)
                                                                                                                                         aCATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCAGGCCATGCAGCGAACCTGACGGCTCGTTa                                >  1:399946/1‑112 (MQ=255)
                                                                                                                                          |                                                                                                                                            
TCCAGCGCAACCGCCAGGTCCGCCTCAGCACGATCCCACATTTCGTCGCTGCCAATACGTTTTTCAGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCAGGCCATGCAGCGAACCTGACGGCTCGTTACGGTGGCAGCTACCAAACTCGGCCATACGC  >  NC_000913/1801234‑1801512

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