Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,612,353 G→A 5.5% T300M (ACG→ATG)  glsB ← glutaminase 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,612,3530GA 5.5% 236.0 / 2.5 73T300M (ACG→ATG) glsBglutaminase 2
Reads supporting (aligned to +/- strand):  ref base G (30/39);  new base A (2/2);  total (32/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.01e-01

TTGTTCGATGACTGGCGCGCCTTTCTCCAGGCAATACTGACGCTCTTTCATGCCCAGACGAAAGCGAGAGCGAAATTTTGAACGGGATAAACGCGCGAAGAGTGGATCGAGAGACTGCATTAATAAACCGAACGCCCTAACTGTTTCGTCAATTGTTCAAGAACGGCAATACCCGCAAGCGAGTTACCTGCATCATCCAGTTCCGGACTCCAGACAGCGATGGCCATTTCATGCGGAACAATCGCCACAATACCGCCACCAACGCCAGATTTCGCCGGTAGCCCCACCCGCCAG  >  NC_000913/1612206‑1612499
                                                                                                                                                   |                                                                                                                                                  
ttgttgGATGACTGGCGCGCCTTTGTCCAGGCAATACTGACGCTCTTTCATGCCAAGACGAAAGCGAGAGCGAAATTTTGAACGGGATAAACGCGCGAAGAGTGGATCGAGAGACTGCATTAATAAACCGAACGCCCTAACTGTTTCg                                                                                                                                                    <  1:111145/148‑1 (MQ=255)
ttgttCGATGACTGGCGCGCCTTTCTCCAGGCAATACTGACGCTCTTTCATGCCCAGACGAAAGCGAGAGCGAAATTTTGAACGGGATAAACGCGCGAAGAGTGGATCGAGAGACTGCATTAATAAACCGAACGCCCTAACTGTTTCg                                                                                                                                                    >  3:4010/1‑148 (MQ=255)
 tgttCGATGACTGGCGCGCCTTTCTCCAGGCAATACTGACGCTCTTTCATGCCCAGACGAAAGCGAGAGCGAAATTTTGAACGGGATAAACGCGCGAAGAGTGGATCGAGAGACTGCATTAATAAACCGAACGCCCTAACTGTTTCGt                                                                                                                                                   >  7:92468/1‑148 (MQ=255)
        tGACTGGCGCGCCTTTCTCCAGGCAATACTGACGCTCTTTCATGCCCAGACGATAGCGAGAGCGAAATTTTGAACGGGATAAACGCGCGAAGAGTGGATCGAGAGACTGCATTAATAAACCGAACGCCCTAACTGTTTCGTCAATTGt                                                                                                                                            <  3:410205/148‑1 (MQ=255)
         caCTGGCGCGCCTTTCTCCAGGCAATACTGACGCTCGTTCCTGCCCAGACGAAAGCGCGAGCGAAATTTTGAACTGGCTAAACGCGCGAAGAGTGTATCGATAGTCTGCATTAATAAACCGAACGCCATAACTGTTTCGTCAATTGt                                                                                                                                            <  4:5821/146‑1 (MQ=255)
          aCTGGCGCGCCTTTCTCCAGGCAATACTGACGCTCTTTCATGCCCAGACGAAAGCGAGAGCGAAATTTTGAACGGGATAAACGCGCGAAGAGTGGATCGAGAGACTGCATTAATAAACCGAACGCCCTAACTGTTTCGTCAATTGtt                                                                                                                                           <  7:32337/147‑1 (MQ=255)
             ggCGCGCCTTTCTCCAGGCAATACTGACGCTCTTTCATGCCCAGACGAAAGCGAGAGCGAAATTTTGAACGGGATAAACGCGCGAAGAGTGGATCGAGAGACTGCATTAATAAACCGAACGCCCTAACTGTTTCGTCAATTGTTc                                                                                                                                          >  1:64983/1‑145 (MQ=255)
              gcgcgcCTTTCTCCAGGCAATACTGACGCTCTTTCATGCCCAGACGAAAGCGAGAGCGAAATTTTGAACGGGATAAACGCGCGAAGAGTGGATCGAGAGACTGCATTAATAAACCGAACGCCCTAACTGTTTCGTCAATTGTTCAAGa                                                                                                                                      <  7:335291/148‑1 (MQ=255)
              gcgcgcCTTTCTCCAGGCAATACTGACGCTCTTTCATGCCCAGACGAAAGCGAGAGCGAAATTTTGAACGGGATAAACGCGCGAAGAGTGGATCGAGAGACTGCATTAATAAACCGAACGCCCTAACTGTTTCGTCAATTGTTCAAGa                                                                                                                                      <  7:361713/148‑1 (MQ=255)
               cgcgcCTTTCTCCAGGCAATACTGACGCTCTTTCATGCCCAGACGAAAGCGAGAGCGAAATTTTGAACGGGATAAACGCGCGAAGAGTGGATCGAGAGACTGCATTAATAAACCGAACGCCCTAACTGTTTCGTCAATTGTTCAAGaa                                                                                                                                     >  2:352786/1‑148 (MQ=255)
                 cgcCTTTCTCCAGGCAATACTGACGCTCTTTCATGCCCAGACGAAAGCGAGAGCGAAATTTTGAACGGGATAAACGCGCGAAGAGTGGATCGAGAGACTGCATTAATAAACCGAACGCCCTAACTGTTTCGTCAATTGTTCAAGAACg                                                                                                                                   <  8:7354/148‑1 (MQ=255)
                  gcCTTTCTCCAGGCAATACTGACGCTCTTTCATGCCCAGACGAAAGCGAGAGCGAAATTTTGAACGGGATAAACGCGCGAAGAGTGGATCGAGAGACTGCATTAATAAACCGAACGCCCTAACTGTTTCGTCAATTGTTCAAGAACgg                                                                                                                                  <  4:136568/148‑1 (MQ=255)
                    cTTTCTCCAGGCAATACTGACGCTCTTTCATGCCCAGACGAAAGCGAGAGCGAAATTTTGAACGGGATAAACGCGCGAAGAGTGGATCGAGAGACTGCATTAATAAACCGAACGCCCTAACTGTTTCGTCAATTGTTCAAGAACGGCa                                                                                                                                <  4:4010/148‑1 (MQ=255)
                      ttCTCCAGGCAATACTGACGCTCTTTCATGCCCAGACGAAAGCGAGAGCGAAATTTTGAACGGGATAAACGCGCGAAGAGTGGATCGAGAGACTGCATTAATAAACCGAACGCCCTAACTGTTTCGTCAATTGTTCAAGAACGGCaa                                                                                                                               >  7:353229/1‑147 (MQ=255)
                       tctcCAGGCAATACTGACGCTCTTTCATGCCCAGACGAAAGCGAGAGCGAAATTTTGAACGGGATAAACGCGCGAAGAGTGGATCGAGAGACTGCATTAATAAACCGAACGCCCTAACTGTTTCGTCAATTGTTCAAGAACGGCAATa                                                                                                                             >  4:288222/1‑148 (MQ=255)
                          ccAGGCAATACAGACGCTCTTTCATGCCCAGACGAAAGCGAGAGCGAAATTTTGAACGGGATAAACGCGCGAAGAGTCGATCGAGACACTGCATTAATAAACCGAACGCCCTAACTGTTTCGTCAATTGTTCAAGAACGGCAATAccc                                                                                                                          >  4:15565/1‑148 (MQ=255)
                               cAATACTGACGCTCTTTCATGCCCAGACGAAAGCGAGAGCGAAATTTTGAACGGGATAAACGCGCGAAGAGTGGATCGAGAGACTGCATTAATAAACCGAACGCCCTACCTGTTTCGTCAATTGTTCAAGAACGGCAATACCCGCAAg                                                                                                                     >  3:329744/1‑148 (MQ=255)
                                 aTACTGACGCTCTTTCATGCCCAGACGAAAGCGAGAGCGAAATTTTGAACGGGATAAACGCGCGAAGAGTGGATCGAGAGACTGCATTAATAAACCGAACGCCCTAACTGTTTCGTCAATTGTTCAAGAACGGCAATACCCGCAAGCg                                                                                                                   <  2:64983/148‑1 (MQ=255)
                                  tACTGACGCTCTTTCATGCCCAGACGAAAGCGAGAGCGAAATTTTGAACGGGATAAACGCGCGAAGAGTGGATCGAGAGACTGCATTAATAAACCGAACGCCCTAACTGTTTCGTCAATTGt                                                                                                                                            >  2:6500/1‑122 (MQ=255)
                                  tACTGACGCTCTTTCATGCCCAGACGAAAGCGAGAGCGAAATTTTGAACGGGATAAACGCGCGAAGAGTGGATCGAGAGACTGCATTAATAAACCGAACGCCCTAACTGTTTCGTCAATTGt                                                                                                                                            <  1:6500/122‑1 (MQ=255)
                                    cTGACGCTCTTTCATTCCCAGCCGAAAGCGAGAGCGAAATTTTGAACGGGATAAACGCGCGAAGAGTCGATCGAGAGACTGCTGTTATAAACCGACCGCCCTAACTGTTTCGTCAATTGTTCAAGAACGGCAATACCCGCAagagag                                                                                                                 <  6:21643/147‑1 (MQ=255)
                                    cTGACGCTCTTTCATGCCCAGACGAAAGCGAGAGCGAAATTTTGAACGGGATAAACGCGCGAAGAGTGGATCGAGAGACTGCATTAATAAACCGAACGCCCTAACTGTTTCGTCAATTGTTCAAGAACGGCa                                                                                                                                <  1:365253/132‑1 (MQ=255)
                                    cTGACGCTCTTTCATGCCCAGACGAAAGCGAGAGCGAAATTTTGAACGGGATAAACGCGCGAAGAGTGGATCGAGAGACTGCATTAATAAACCGAACGCCCTAACTGTTTCGTCAATTGTTCAAGAACGGCa                                                                                                                                >  2:365253/1‑132 (MQ=255)
                                          ctctTTCATGCCCAGACGAAAGCGAGAGCGAAATTTTGAACGGGATAAACGCGCGAAGAGTGGATCGAGAGACTGCATTAATAAACCGAACGCCCTAACTGTTTCGTCAATTGTTCAAGAACGGCAATACCCGCAAGCGAGTTACCTg                                                                                                          >  2:10919/1‑148 (MQ=255)
                                            atttCATGCCCAGACGAAAGCGAGAGAGAAATTTTTAACGGTATTAACGCGAGAAGAGTGGATCGAGAGACTGCATTAATAAACCGAACGCCCTAACTGTTTCGTCAATTGTTCAAGAACGGCAATACCCGCAAGCGAGTTACCTGc                                                                                                         <  8:353229/146‑1 (MQ=255)
                                              ttCATGCCCAGACGAAAGCGAGAGCGAAATTTTGAACGGGATAAACGCGCGAAGAGTGGATCGAGAGACTGCATTAATAAACCGAACGCCCTAACTGTTTCGTCAATTGTTCAAGAACGGCAATACCCGCAAGCGAGTTACCTGcat                                                                                                       >  7:369384/1‑147 (MQ=255)
                                                cATGCCCAGACGAAAGCGAGAGCGAAATTTTGAACGGGATAAACGCGCGAAGAGTGGATCGAGAGACTGCATTAATAAACCGAACGCCCTAACTGTTTCGTCAATTGTTCAAGAACgg                                                                                                                                  <  3:370652/118‑1 (MQ=255)
                                                cATGCCCAGACGAAAGCGAGAGCGAAATTTTGAACGGGATAAACGCGCGAAGAGTGGATCGAGAGACTGCATTAATAAACCGAACGCCCTAACTGTTTCGTCAATTGTTCAAGAACgg                                                                                                                                  >  4:370652/1‑118 (MQ=255)
                                                cATGCCCAGACGAAAGCGAGAGCGAAATTTTGAACGGGATAAACGCGCGAAGAGTGGATCGAGAGACTGCATTAATAAACCGAACGCCCTAACTGTTTCGTCAATTGTTCAAGAACGGCAATACCCGCAAGCGAGTTACCTGcatcat                                                                                                    >  6:247869/1‑148 (MQ=255)
                                                     ccAGACGAAAGCGAGAGCGAAATTTTGAACGGGATAAACGCGCGAAGAGTGGATCGAGAGACTGCATTAATAAACCGAACGCCCTAACTGTTTCATCAATTGTTCAAGAACGGCAATACCCGCAAGCGAGTTACCTGcatcat                                                                                                    >  7:321338/1‑143 (MQ=255)
                                                     ccAGACGAAAGCGAGAGCGAAATTTTGAACGGGATAAACGCGCGAAGAGTGGATCGAGAGACTGCATTAATAAACCGAACGCCCTAACTGTTTCATCAATTGTTCAAGAACGGCAATACCCGCAAGCGAGTTACCTGcatcat                                                                                                    <  8:321338/143‑1 (MQ=255)
                                                         aCGAAAGCGAGAGCGAAATTTTGAACGGGATAAACGCGCGAAGAGTGGATCGAGAGACTGCATTAATAAACCGAACGCCCTAACTGTTTCGTCAATTGTTCAAGAACGGCAATACCCGCAAGCGAGTTACCTGCATCATCCAGTTCCg                                                                                           >  2:38267/1‑148 (MQ=255)
                                                          cGAAAGCGAGAGCGAAATTTTGAACGGGATAAACGCGCGAAGAGTGGATCGAGAGACTGCATTAATAAACCGAACGCCCTAACTGTTTCGTCAATTGTTCAAGAACGGCa                                                                                                                                <  3:156505/110‑1 (MQ=255)
                                                          cGAAAGCGAGAGCGAAATTTTGAACGGGATAAACGCGCGAAGAGTGGATCGAGAGACTGCATTAATAAACCGAACGCCCTAACTGTTTCGTCAATTGTTCAAGAACGGCa                                                                                                                                >  4:156505/1‑110 (MQ=255)
                                                            aaaGCGAGAGCGAAATTTTGAACGGGATAAACGCGCGAAGAGTGGATCGAGAGACTGCATTAATAAACCGAACGCCCTAACTGTTTCGTCAATTGTTCAAGAACGGCAATACCCGCAAGCGAGTTACCTGc                                                                                                         >  8:19866/1‑131 (MQ=255)
                                                            aaaGCGAGAGCGAAATTTTGAACGGGATAAACGCGCGAAGAGTGGATCGAGAGACTGCATTAATAAACCGAACGCCCTAACTGTTTCGTCAATTGTTCAAGAACGGCAATACCCGCAAGCGAGTTACCTGc                                                                                                         <  7:19866/131‑1 (MQ=255)
                                                            aaaGCGAGAGCGAAATTTTGAACGGGATAAACGCGCGAAGAGTGGATCGAGAGACTGCATTAATAAACCGAACGCCCTAACTGTTTCATCAATTGTTCAAGAACGGCAATACCCGCAAGCGAGTTACCTGCATCATCCAGTTCCGGAc                                                                                        >  6:269086/1‑148 (MQ=255)
                                                                 gagaGCGAAATTTTGAACGGGATAAACGCGCGAAGAGTGGATCGAGAGACTGCATTAATAAACCGAACGCCCTAACTGTTTCGTCAATTGTTCAAGAACGGCAATACCCGCAAGCGAGTTACCTGCATCATCCAGTTCCGGACTCCAg                                                                                   <  5:415296/148‑1 (MQ=255)
                                                                   gagCGAAATTTTGAACGGGATAAACGCGCGAAGAGTGGATCGAGAGACTGCATTAATAAACCGAACGCCCTAACTGTTTAGTCAATTGTTCAAGAACGGCAATACCCGCAAGCGAGTTACCTGCATCATCCAGTTCCGGACTCCAGa                                                                                  >  6:46157/1‑147 (MQ=255)
                                                                    agCGAAATTTTGAACGGGATAAACGCGCGAAGAGTGGATCGAGAGACTGCATTAATATACCGAAAGCCCTAACTGTTTCGTCAATTGTTCAAGAACGGCAATACCCGCAACCGAGTTAAATGCATCATCCAGTTCCGGACTCCAGa                                                                                  >  1:109223/1‑146 (MQ=255)
                                                                       gAAATTTTGAACGGGATAAACGCGCGAAGAGTGGATCGAGAGACTGCATTAATAAACCGAACGCCCTAACTGTTTTGTCAATTGTTCAAGAACGGCAATACCCGCAAGCGAGTTACCTGCATCATCCAGTTCCGGACTCCAGACAGCg                                                                             <  5:46157/148‑1 (MQ=255)
                                                                       gAAATTTTGAACGGGATAAACGCGCGAAGAGTGGATCGAGAGACTGCATTAATAAACCGAACGCCCTAACTGTTTCGTCAATTGTTCAAGAACGGCAATACCCGCAAGCGAGTTACCTGCATCATCCAGTTCCGGACTCCAGACAGCg                                                                             <  4:280182/148‑1 (MQ=255)
                                                                       gAAATTTTGAACGGGATAAACGCGCGAAGAGTGGATCGAGAGACTGCATTAATAAACCGAACGCCCTAACTGTTTCATCAATTGTTCAAGAACGGCAATACCCGCAAGCGAGTTACCTGCATCATCCAGTTCCGGACTCCAGACAGCg                                                                             <  5:269086/148‑1 (MQ=255)
                                                                        aaaTTTTGAACGGGATAAACGCGCGAAGAGTGGATCGAGAGACTGCATTAATAAACCGAACGCCCTAACTGTTTCGTCAATTGTTCAAGAACGGCAATACCCGCAAGCGAGTTACCTGCATCATCCAGTTCCGGACTCCAGACAGCGa                                                                            <  8:369384/148‑1 (MQ=255)
                                                                         aaTTTTGAACGGTATAAACGCGCGAAGAGTGTATCGAGAGACTGCATTAATAAACCGAACGCCCTAACTGTTTCGTCAATTGTTCAAGAACGGCAATACCCGCAAGCGAGTTACCTGCATCATCCAGTTCCGGACTCCAGACAGCGAt                                                                           <  7:53533/148‑1 (MQ=255)
                                                                              tGAACGGGATAAACGCGCGAAGAGTGGATCGAGAGACTGCATTAATAAACCGAACGCCCTAACTGTTTCGTCAATTGTTCAAGAACGGCAATACCCGCAAGCGAGTTACCTGCATCATCCAGTTCCGGACTCCAGACAGCGATGGcc                                                                       <  3:288222/147‑1 (MQ=255)
                                                                                     gATAAACGCGCGAAGAGTGGATCGAGAGACTGCATTAATAAACCGAACGCCCTAACTGTTTCGTCAATTGTTCAAGAACGGCAATACCCGCAAGCGAGTTACCTGCATCATCCAGTTCCGGACTCCAGACAGCGATGGCCATTTCATg                                                               <  1:38267/148‑1 (MQ=255)
                                                                                         aaCGCGCGAAGAGTGGATCGAGAGACTGCATTAATAAACCGAACGCCCTAACTGTTTCGTCAATTGTTCAAGAACGGCAATACCCGCAAGCGAGTTAcc                                                                                                            >  2:204499/1‑99 (MQ=255)
                                                                                         aaCGCGCGAAGAGTGGATCGAGAGACTGCATTAATAAACCGAACGCCCTAACTGTTTCGTCAATTGTTCAAGAACGGCAATACCCGCAAGCGAGTTAcc                                                                                                            <  1:204499/99‑1 (MQ=255)
                                                                                         aaCGCGCGAAGAGTGGATCGAGAGACTGCATTAATAAACCGAACGCCCTAACTGTTTCGTCAATTGTTCAAGAACGGCAATACCCGCAAGCGAGTTACCTGCATCATCCAGTTCCGGACTCCAGACAGCGATGGc                                                                        >  5:100132/1‑135 (MQ=255)
                                                                                         aaCGCGCGAAGAGTGGATCGAGAGACTGCATTAATAAACCGAACGCCCTAACTGTTTCGTCAATTGTTCAAGAACGGCAATACCCGCAAGCGAGTTACCTGCATCATCCAGTTCCGGACTCCAGACAGCGATGGc                                                                        <  6:100132/135‑1 (MQ=255)
                                                                                             cgcgAAGAGTGGATCGAGAGACTGCATTAATAAACCGAACGCCCTAACTGTTTCGTCAATTGTTCAAGAACGGCAATACCCGCAAGCGAGTTACCTGCATCATCCAGTTCCGGACTCCAGACAGCGATGGCCATTTCATGCGGAACaa                                                       <  5:247869/148‑1 (MQ=255)
                                                                                               cgAAGAGTGGATCGAGAGACTGCATTAATAAACCGAACGCCCTAACTGTTTCGTCAATTGTTCAAGAACGGCAATACCCGCAAGCGAGTTACCTGCATCATCCAGTTCCGGACTCCAGACa                                                                                <  7:42511/121‑1 (MQ=255)
                                                                                               cgAAGAGTGGATCGAGAGACTGCATTAATAAACCGAACGCCCTAACTGTTTCGTCAATTGTTCAAGAACGGCAATACCCGCAAGCGAGTTACCTGCATCATCCAGTTCCGGACTCCAGACa                                                                                >  8:42511/1‑121 (MQ=255)
                                                                                                 cagagTGCATCGAGAGACTGCATTAATACACCGAACGCCCTAACTGTTTCGTCAATTGTTCAAGAACGGCAATACCCGCAAGCGAGTTACCTGCATCATCCAGTTCCGGACTCCAGACAGCGATGGCCATTTCATGCGGAACAATCGc                                                   <  8:92468/147‑1 (MQ=255)
                                                                                                     gTGGATCGAGAGACTGCATTAATAAACCGAACGCCCTAACTGTTTCGTCAATTGTTCAAGAACGGCAATACCCGCAAGCGAGTTACCTGCATCATCCAGTTCCGGACTCCAGACAGCGATGGCCATTTCATGCGGAACAATCGCcaca                                               <  4:329744/148‑1 (MQ=255)
                                                                                                         aTCGAGAGACTGCATTAATAAACCGAACGCCCTAACTGTTTCGTCAATTGTTCAAGAACGGCAATACCCGCAAGCGAGTTACCTGCATCATCCAGTTCCGGACTCCAGACAGCGATGGCCATTTCATGCGGAACAATCGCCACAATAc                                           <  1:352786/148‑1 (MQ=255)
                                                                                                                      aTTAATAAACCGAACGCCCTAACTGTTTCGTCAATTGTTCAAGAACGGCAATACCCGCAAGCGAGTTACCTGcat                                                                                                       <  5:819/75‑1 (MQ=255)
                                                                                                                      aTTAATAAACCGAACGCCCTAACTGTTTCGTCAATTGTTCAAGAACGGCAATACCCGCAAGCGAGTTACCTGcat                                                                                                       >  6:819/1‑75 (MQ=255)
                                                                                                                        taataaACCGAACGCCCTAACTGTTTCGTCAATTGTTCAAGAACGGCAATACCCGCAAGCGAGTTACCTGCATCATCCAGTTCCGGACTCCAGACAGCGATGGCCATTTCATGCGGAACAATCGCCACAATACCGCCACCAACGc                               >  3:137625/1‑145 (MQ=255)
                                                                                                                                   aCGCCCTAACTGGTTCGTCAATTGTTCAAGAACGGCAATACCCGCAAGCGAGTTACCTGCATCATCCAGTTCCGGACTCCAGACAGCGATGGCCATTTCATGCGGAACAATCGCCACAATACCGCCACCAACGCCAGATTTCGCCGGt                 <  2:187521/148‑1 (MQ=255)
                                                                                                                                     gCCCTAACTGTTTCGTCAATTGTTCAAGAACGGCAATACCCGCAAGCGAGTTACCTGCATCATCCAGTTCCGGACTCCAGACAGCGATGGCCATTTCATGCGGAACAATCGCCACAATACCGCCACCAACGCCAGATTTCGCCGGTAg               <  3:15565/148‑1 (MQ=255)
                                                                                                                                        cTAACTGTTTCGTCAATTGTTCAAGAACGGCAATACCCGCAAGCGAGTTACCTGCATCATCCAGTTCCGGACTCCAGACAGCGATGGCCATTTCATGCGGAACAATCGCCACAATACCGCCACCAACGCCAGATTTCGCCGGTAGccc            >  7:410232/1‑148 (MQ=255)
                                                                                                                                         tAACTGTTTCGTCAATTGTTCAAGAACGGCAATACCCGCAAGCGAGTTACCTGCATCATCCAGTTCCGGACTCCAGACAGCGATGGCCATTTCATGCGGAACAATCGCCACAATACCGCCACCACCGCCAGATTTCGCCGGTAGcccc           >  6:50206/1‑148 (MQ=255)
                                                                                                                                         tAACTGTTTCGTCAATTGTTCAAGAACGGCAATACCCGCAAGCGAGTTACCTGCATCATCCAGTTCCGGACTCCAGACAGCGATGGCCATTTCATGCGGAACAATCGCCACAATACCGCCACCAACGCCAGATTTCGCCGGTAGccc            >  3:179434/1‑147 (MQ=255)
                                                                                                                                          aaCTGTTTCGTCAATTGTTCAAGAACGGCAATACCCGCAAGCGAGTTACCTGCATCATCCAGTTCCGGACTCCAGACAGCg                                                                             <  7:3326/81‑1 (MQ=255)
                                                                                                                                          aaCTGTTTCGTCAATTGTTCAAGAACGGCAATACCCGCAAGCGAGTTACCTGCATCATCCAGTTCCGGACTCCAGACAGCg                                                                             >  8:3326/1‑81 (MQ=255)
                                                                                                                                          aaCTGTTTCGTCAATTGTTCAAGAACGGCAATACCCGCAAGCGAGTTACCTGCATCATCCAGTTCCGGACTCCAGACAGCGATGGCCATTTc                                                                  >  4:95560/1‑92 (MQ=255)
                                                                                                                                          aaCTGTTTCGTCAATTGTTCAAGAACGGCAATACCCGCAAGCGAGTTACCTGCATCATCCAGTTCCGGACTCCAGACAGCGATGGCCATTTc                                                                  <  3:95560/92‑1 (MQ=255)
                                                                                                                                           aCTGTTTCGTCAATTGTTCAAGAACGGCAATACCCGCAAGCGAGTTACCTGCATCATCCAGTTCCGGACTCCAGACAGCGATGGCCATTTCATGCGGAACAATCGCCACAATACCGCCACCAACGCCAGATTTCGCCGGTAGcccc           <  4:179434/146‑1 (MQ=255)
                                                                                                                                                ttCGTCAATTGTTCAAGAACGGCAATACCCGCAAGCGAGTTACCTGCATCATCCAGTTCCGGACTCCAGACAGCGATGGCCATTTCATGCGGAACAATCGCCACAATACCGCCACCaa                                  >  7:56072/1‑118 (MQ=255)
                                                                                                                                                ttCGTCAATTGTTCAAGAACGGCAATACCCGCAAGCGAGTTACCTGCATCATCCAGTTCCGGACTCCAGACAGCGATGGCCATTTCATGCGGAACAATCGCCACAATACCGCCACCaa                                  <  8:56072/118‑1 (MQ=255)
                                                                                                                                                 tCGTCAATTGTTCAAGAACGGCAATACCCGCAAGCGAGTTACCTGCATCATCCAGTTCCGGACTCCAGACAGCGATGGCCATTTCATGCGGAACAATCGCCACAATACCGCCACCAACGCCAGATTTCGCCGGTAGCCCCACCCGcc    <  1:10919/147‑1 (MQ=255)
                                                                                                                                                  cGTCAATTGTTCAAGAACGGCAATACCCGCAAGCGAGTTACCTGCATCATCCAGTTCCGGACTCCAGACAGCGATGGCCATTTCATGCGGAACAATCGCCACAATACCGccacca                                   <  5:61170/115‑1 (MQ=255)
                                                                                                                                                  cGTCAATTGTTCAAGAACGGCAATACCCGCAAGCGAGTTACCTGCATCATCCAGTTCCGGACTCCAGACAGCGATGGCCATTTCATGCGGAACAATCGCCACAATACCGccacca                                   >  6:61170/1‑115 (MQ=255)
                                                                                                                                                   gTCAATTGTTCAAGAACGGCAATACCCGCAAGCGAGTTACCTGCATCATCCAGTTCCGGACTCCAGACAGCGATGGCCATTTCATGCGGAACAATCGCCACAATACCGCCACCAACGCCAGATTTCGCCGGTAGCCCCACCCGCCAg  <  8:410232/147‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TTGTTCGATGACTGGCGCGCCTTTCTCCAGGCAATACTGACGCTCTTTCATGCCCAGACGAAAGCGAGAGCGAAATTTTGAACGGGATAAACGCGCGAAGAGTGGATCGAGAGACTGCATTAATAAACCGAACGCCCTAACTGTTTCGTCAATTGTTCAAGAACGGCAATACCCGCAAGCGAGTTACCTGCATCATCCAGTTCCGGACTCCAGACAGCGATGGCCATTTCATGCGGAACAATCGCCACAATACCGCCACCAACGCCAGATTTCGCCGGTAGCCCCACCCGCCAG  >  NC_000913/1612206‑1612499

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: