Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,540,347 C→T 5.9% G748G (GGG→GGA narZ ← nitrate reductase 2 (NRZ), alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,540,3470CT 5.9% 218.3 / 2.5 68G748G (GGG→GGAnarZnitrate reductase 2 (NRZ), alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (34/30);  new base T (2/2);  total (36/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.88e-01

CATATCGTCTTTTTCGTACCAGGTGGCGGTGGGCAGAACGATATCGGAGAACAGGCAGGTACTGGACATGCGGAAGTCGAGCGTCACCAGCAGGTCGAGCTTGCCTTCAATCGCTGCAGTTTGCCACTCGACTTCTTCCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCT  >  NC_000913/1540200‑1540494
                                                                                                                                                   |                                                                                                                                                   
cATATCGTCTTTTTCGTACCAGGTGGCGGTGGGCAGAACGATATCGGAGAACAGGCAGGTACTGGACATGCGGAAGTCGAGCGTCACCAGCAGGTCGAGCTTGCCTTCAATCGCTGCAGTTTGCCACTCGACTTCTTCCGGTTTGATc                                                                                                                                                     <  1:27624/148‑1 (MQ=255)
   atCGTCTTTTTCGTACCAGGTGGCGGTGGGCAGAACGATATCGGAGAACAGGCAGGTACTGGACATGCGGAAGTCGAGCGTCACCAGCAGGTCGAGCTTGCCTTCAATCGCTGCAGTTTGCCACTCGACTTCTTCCGGTTTGATCCCGt                                                                                                                                                 >  6:12977/1‑149 (MQ=255)
       tCTTTTTCGTACCAGGTGGCGGTGGGCAGAACGATATCGGAGAACAGGCAGGTACTGGACATGCGGAAGTCGAGCGTCACCAGCAGGTCGAGCTTGCCTTCAATCGCTGCAGTTTGCCACTCGACTTCTTCCGGTTTGATCCCGTCGCt                                                                                                                                             <  8:351315/149‑1 (MQ=255)
         tttttCGTACCAGGTGGCGGTGGGCAGAACGATATCGGAGAACAGGCAGGTACTGGACATGCGGAAGTCGAGCGTCACCAGCAGGTCGAGCTTGCCTTCAATCGCTGCAGTTTGCCACTCGACTTCTACCGGTTTGATCCCGTCGCTgg                                                                                                                                           >  1:46332/1‑149 (MQ=255)
                                  aGAACGATATCGGAGAACAGGCAGGTACTGGACATGCGGAAGTCGAGCGTCACCAGCAGGTCGAGCTTGCCTTCAATCGCTGCAGTTTGCCACTCGACTTCTTCCGGTTTGATCCCGTCGCTGGCACCGAGt                                                                                                                                   >  7:257186/1‑132 (MQ=255)
                                  aGAACGATATCGGAGAACAGGCAGGTACTGGACATGCGGAAGTCGAGCGTCACCAGCAGGTCGAGCTTGCCTTCAATCGCTGCAGTTTGCCACTCGACTTCTTCCGGTTTGATCCCGTCGCTGGCACCGAGt                                                                                                                                   <  8:257186/132‑1 (MQ=255)
                                       gATATCGGAGAACAGGCAGGTACTGGACATGCGGAAGTCGAGCGTCACCAGCAGGTCGAGCTTGCCTTCAATCGCTGCAGTTTGCCACTCGACTTCTTCCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCttt                                                                                                             <  7:96787/149‑1 (MQ=255)
                                              gagaACAGGCAGGTACTGGACATGCGGAAGTCGAGCGTCACCAGCAGGTCGAGCTTGCCTTCAATCGCTGCAGTTTGCCACTCGACTTCTTCCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTccc                                                                                                      >  5:143303/1‑149 (MQ=255)
                                                          gTACTGGACATGCGGAAGTCGAGCGTCACCAGCAGGTCGAGCTTGCCTTCAATCGCTGCAGTTTGCCACTCGACTTCTTCCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACtt                                                                                          >  8:299495/1‑149 (MQ=255)
                                                              tGGACATGCGGAAGTCGAGCGTCACCAGCAGGTCGAGCTTGCCTTCAATCGCTGCAGTTTGCCACTCGACTTCTTCCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATccc                                                                                                                  >  8:353988/1‑121 (MQ=255)
                                                              tGGACATGCGGAAGTCGAGCGTCACCAGCAGGTCGAGCTTGCCTTCAATCGCTGCAGTTTGCCACTCGACTTCTTCCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATccc                                                                                                                  <  7:353988/121‑1 (MQ=255)
                                                              tGGACATGCGGAAGTCGAGCGTCACCAGCAGGTCGAGCTTGCCTTCAATCGCTGCAGTTTGCCACTCGACTTCTTCCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTc                                                                                         <  5:274224/146‑1 (MQ=255)
                                                              tGGACATGCGGAAGTCGAGCGTCACCAGCAGGTCGAGCTTGCCTTCAATCGCTGCAGTTTGCCACTCGACTTCTTCCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTc                                                                                         >  6:274224/1‑146 (MQ=255)
                                                                   aTGCGGAAGTCGAGCGTCACCAGCAGGTCGAGCTTGCCTTCAATCGCTGCAGTTTGCCACTCGACTTCTTCCGGTTTGATTCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGc                                                                                      <  8:254312/144‑1 (MQ=255)
                                                                   aTGCGGAAGTCGAGCGTCACCAGCAGGTCGAGCTTGCCTTCAATCGCTGCAGTTTGCCACTCGACTTCTTCCGGTTTGATTCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGc                                                                                      <  8:248428/144‑1 (MQ=255)
                                                                   aTGCGGAAGTCGAGCGTCACCAGCAGGTCGAGCTTGCCTTCAATCGCTGCAGTTTGCCACTCGACTTCTTCCGGTTTGATTCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGc                                                                                      >  7:248428/1‑144 (MQ=255)
                                                                   aTGCGGAAGTCGAGCGTCACCAGCAGGTCGAGCTTGCCTTCAATCGCTGCAGTTTGCCACTCGACTTCTTCCGGTTTGATTCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGc                                                                                      >  7:254312/1‑144 (MQ=255)
                                                                   aTGCGGAAGTCGAGCGTCACCAGCAGGTCGAGCTTGCCTTCAATCGCTGCAGTTTGCCACTCGACTTCTTCCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCAtata                                                                                 >  4:94626/1‑149 (MQ=255)
                                                                    tGCGGAAGTCGAGCGTCACCAGCAGGTCGAGCTTGCCTTCAATCGCTGCAGTTTGCCACTCGACTTCTTCCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATAc                                                                                <  5:12977/149‑1 (MQ=255)
                                                                       ggAAGTCGAGCGTCACCAGCAGGTCGAGCTTGCCTTCAATCGCTGCAGTTTGCCACTCGACTTCTTCCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCg                                                                             >  1:144683/1‑149 (MQ=255)
                                                                        gAAGTCGAGCGTCACCAGCAGGTCGAGCTTGCCTTCAATCGCTGCAGTTTGCCACTCGACTTCTTCCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGt                                                                            <  6:282863/149‑1 (MQ=255)
                                                                        gAAGTCGAGCGTCACCAGCAGGTCGAGCTTGCCTTCAATCGCTGCAGTTTGCCACTCGACTTCTTCCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGt                                                                            >  2:325681/1‑149 (MQ=255)
                                                                            tCGAGCGTCACCAGCAGGTCGAGCTTGCCTTCAATCGCTGCAGTTTGCCACTCGACTTCTTCCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCcagca                                                                                                 >  6:127949/1‑124 (MQ=255)
                                                                            tCGAGCGTCACCAGCAGGTCGAGCTTGCCTTCAATCGCTGCAGTTTGCCACTCGACTTCTTCCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCcagca                                                                                                 <  5:127949/124‑1 (MQ=255)
                                                                              gAGCGTCACCAGCAGGTCGAGCTTGCCTTCAATCGCTGCAGTTTGCCACTCGACTTCTTCCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCtt                                                                      >  5:53892/1‑149 (MQ=255)
                                                                              gAGCGTCACCAGCAGGTCGAGCTTGCCTTCAATCGCTGCAGTTTGCCACTCCACTTCTTCCGGTTTGAACCCGTCGCTGGCACCGAGTTCCTCGCCATGAATCCCGCTTTAGGTCCACAGCAGATACTTCTGCTAATACTCGTGGCCtt                                                                      >  6:20568/1‑149 (MQ=255)
                                                                                   tCACCAGCAGGTCGAGCTTGCCTTCAATCGCTGCAGTTTGCCACTCGACTTCTTCCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCg                                                                 >  6:107786/1‑149 (MQ=255)
                                                                                   tCACCAGCAGGTCGAGCTTGCCTTCAATCGCTGCAGTTTGCCACTCGACTTCTTCCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCg                                                                 >  6:102722/1‑149 (MQ=255)
                                                                                     aCCAGGAGGGCGAGCTTGCCTTCAATCGCTGCAGTTTGCCACTCGACTTCTTCCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCAGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCgga                                                               <  5:20568/149‑1 (MQ=255)
                                                                                      ccAGCAGGTCGAGCTTGCCTTCAATCGCTGCAGTTTGCCACTCGACTTCTTCCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCggag                                                              >  7:257382/1‑149 (MQ=255)
                                                                                       cagcagGTCGAGCTTGCCTTCAATCGCTGCAGTTTGCCACTCGACTTCTTCCGGTTTGATCCCGTCGCt                                                                                                                                             <  6:138735/69‑1 (MQ=255)
                                                                                         gcagGTCGAGCTTGCCTTCAATCGCTGCAGTTTGCCACTCGACTTCTTCCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGgaggag                                                           >  5:118978/1‑149 (MQ=255)
                                                                                         gcagGTCGAGCTTGCCTTCAATCGCTGCAGTTTGCCACTCGACTTCTTCCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGgaggag                                                           >  2:294520/1‑149 (MQ=255)
                                                                                         gcagGTCGAGCTTGCCTTCAATCGCTGCAGTTTGCCACTCGACTTCTTCCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCAGCATATACTCGTGGCCTTTGCCGgaggag                                                           <  7:299495/149‑1 (MQ=255)
                                                                                                gAGCTTGCCTTCAATCGCTGCAGTTTGCCACTCGACTTCTTCCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAg                                                      >  1:83025/1‑147 (MQ=255)
                                                                                                      gCCTTCAATCGCTGCAGTTTGCCACTCGACTTCTTCCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGACTCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGt                                                 >  3:303644/1‑146 (MQ=255)
                                                                                                        cTTCAATCGCTGCAGTTTGCCACTCGACTTCTTCCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGc                                                                                                                >  5:277908/1‑81 (MQ=255)
                                                                                                        cTTCAATCGCTGCAGTTTGCCACTCGACTTCTTCCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGc                                                                                                                <  6:277908/81‑1 (MQ=255)
                                                                                                        cTTCAATCGCTGCAGTTTGCCACTCGACTTCTTCCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTagag                                            <  6:118978/149‑1 (MQ=255)
                                                                                                         ttCAATCGCTGCAGTTTGCCACTCGACTTCTTCCGGTTTGATCCTGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAgc                                           >  6:183937/1‑149 (MQ=255)
                                                                                                         ttCAATCGCTGCAGTTTGCCACTCGACTTCTTCCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCg                                                                                                           <  1:374548/85‑1 (MQ=255)
                                                                                                         ttCAATCGCTGCAGTTTGCCACTCGACTTCTTCCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCg                                                                                                           >  2:374548/1‑85 (MQ=255)
                                                                                                         ttCAATCGCTGCAGTTTGCCACTCGACTTCTTCCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAgc                                           >  8:76645/1‑149 (MQ=255)
                                                                                                         ttCAATCGCTGCAGTTTGCCACTCGACTTCTTCCGGTTTGAGCCCGTCGCTGGCACCGAGTTGCTCGCCCTGAATCCCGCTTTCGGTCACCAGCAGATACTTCTGCGTATACTCGTGGCCTTTGGCGGTGGAGCCAAGCAGGTTAGAgc                                           >  5:207613/1‑149 (MQ=255)
                                                                                                            aaTCGCTGCAGTTTGCCACTCGACTTCTTCCGGTTTGATCCCGTCGCTGGCACCGAGGTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAgcgc                                         >  5:45907/1‑148 (MQ=255)
                                                                                                              tCGCTGCAGTTTGCCACTCGACTTCTTCCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCgga                                                               <  1:214031/124‑1 (MQ=255)
                                                                                                              tCGCTGCAGTTTGCCACTCGACTTCTTCCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCgga                                                               >  2:214031/1‑124 (MQ=255)
                                                                                                                gCTGCAGTTTGCCACTCGACTTCTTCCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGac                                    <  1:325681/149‑1 (MQ=255)
                                                                                                                gCTGCAGTTTGCCACTCGACTTCTTCCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGac                                    <  2:83025/149‑1 (MQ=255)
                                                                                                                  tGCAGTTTGCCACTCGACTTCTTCCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGaa                                                                                                                      <  4:86035/65‑1 (MQ=255)
                                                                                                                  tGCAGTTTGCCACTCGACTTCTTCCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGaa                                                                                                                      >  3:86035/1‑65 (MQ=255)
                                                                                                                          gCCACTCGACTTCTTCCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAAtt                           <  6:143303/148‑1 (MQ=255)
                                                                                                                           ccACTCGACTTCTTCCGGTTCGATCCCGTCGCTGGCACCGAGTGCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCGTATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTAc                         <  6:207613/149‑1 (MQ=255)
                                                                                                                            cACTCGACTTCTTCCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCaa                                                       <  3:115912/118‑1 (MQ=255)
                                                                                                                            cACTCGACTTCTTCCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCaa                                                       >  4:115912/1‑118 (MQ=255)
                                                                                                                             aCTCGACTTCTTCCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAg                                                   >  3:276967/1‑121 (MQ=255)
                                                                                                                             aCTCGACTTCTTCCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAg                                                   <  4:276967/121‑1 (MQ=255)
                                                                                                                              cTCGACTTCTTCCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTAcgcg                      <  3:94626/149‑1 (MQ=255)
                                                                                                                              cTCGACTTCTTCCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTAcgcg                      >  6:93106/1‑149 (MQ=255)
                                                                                                                                cGACTTCTTCCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGc                                                     >  2:364059/1‑116 (MQ=255)
                                                                                                                                cGACTTCTTCCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGc                                                     <  1:364059/116‑1 (MQ=255)
                                                                                                                                 gACTTCTTCCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGgat                   <  5:107786/149‑1 (MQ=255)
                                                                                                                                 gACTTCTTCCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGgat                   <  5:102722/149‑1 (MQ=255)
                                                                                                                                       ttcCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAg                                                   <  4:181113/111‑1 (MQ=255)
                                                                                                                                       ttcCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAg                                                   >  3:181113/1‑111 (MQ=255)
                                                                                                                                          cGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATtgctgc          <  6:53892/149‑1 (MQ=255)
                                                                                                                                          cGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATtgctgc          <  8:257382/149‑1 (MQ=255)
                                                                                                                                           ggTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGTAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATtgctgct         <  2:144683/149‑1 (MQ=255)
                                                                                                                                                  tCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATTCTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCt  >  6:118088/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CATATCGTCTTTTTCGTACCAGGTGGCGGTGGGCAGAACGATATCGGAGAACAGGCAGGTACTGGACATGCGGAAGTCGAGCGTCACCAGCAGGTCGAGCTTGCCTTCAATCGCTGCAGTTTGCCACTCGACTTCTTCCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCT  >  NC_000913/1540200‑1540494

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