Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,606,094 (G)8→7 noncoding (72/87 nt) leuV ← tRNA‑Leu

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,606,0870G.96.8% 274.9 / 0.2 63noncoding (79/87 nt)leuVtRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base G (0/1);  major base . (33/27);  minor base C (0/2);  total (33/30)
Fisher's exact test for biased strand distribution p-value = 2.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCTGTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAG‑GAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTT  >  NC_000913/4605940‑4606219
                                                                                                                                                    |                                                                                                                                    
tgCTGTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAG‑GAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAgggg                                                                                                                                   >  1:202878/1‑151 (MQ=255)
    gTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAG‑GAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAgggggg                                                                                                                                 <  2:510844/149‑1 (MQ=255)
      aTCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAG‑GAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAgggggggact                                                                                                                             >  1:768273/1‑148 (MQ=255)
        cacaAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAG‑GAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTg                                                                                                                          <  2:769220/151‑1 (MQ=255)
        cacaAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTCATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAG‑GAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAgggggggactt                                                                                                                            >  1:338860/1‑146 (MQ=255)
             aCAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAG‑GAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAAccc                                                                                                                     >  2:742138/1‑151 (MQ=255)
              cAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAG‑GAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAAcccc                                                                                                                    >  1:824264/1‑151 (MQ=255)
              cAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAG‑GAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAAcccc                                                                                                                    >  1:255120/1‑151 (MQ=255)
                gTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAG‑GAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCa                                                                                                                  <  2:137994/151‑1 (MQ=255)
                  tGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAG‑GAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCAc                                                                                                                 <  2:315857/150‑1 (MQ=255)
                  tGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAG‑GAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACg                                                                                                                >  1:281160/1‑151 (MQ=255)
                    cGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAG‑GAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTc                                                                                                              >  1:374586/1‑151 (MQ=255)
                    cGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAG‑GAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTc                                                                                                              <  2:36875/151‑1 (MQ=255)
                     ggCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAG‑GAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTcc                                                                                                             <  2:253169/151‑1 (MQ=255)
                     ggCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAG‑GAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTcc                                                                                                             >  2:332922/1‑151 (MQ=255)
                           aCATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAG‑GAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAgg                                                                                                       >  1:834952/1‑151 (MQ=255)
                           aCATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAG‑GAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAg                                                                                                        <  1:252512/150‑1 (MQ=255)
                             aTTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAG‑GAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGac                                                                                                     <  1:389345/151‑1 (MQ=255)
                             aTTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAG‑GAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGac                                                                                                     <  1:738071/151‑1 (MQ=255)
                                   cTACAATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAG‑GAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAac                                                                                               >  1:660001/1‑151 (MQ=255)
                                       cATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAG‑GAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCt                                                                                           >  1:649232/1‑151 (MQ=255)
                                                 ggCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAG‑GAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAg                                                                                   <  1:629244/149‑1 (MQ=255)
                                                  gcagaaaTTTCTGGCATCATGCTAATGATAATGAATTTGAAGCAGGGACCAAAAA‑CAAGAGAAGAAACAACAAAAGAAGTAAAGACAATTAGTTCGACGGAGGGGACTTGAACCCCCACGTCCGTAAGGACAATAACGCCTGAAGCTa                                                                                    <  1:779636/142‑1 (MQ=19)
                                                     cTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAG‑GAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTc                                                                             >  1:89782/1‑151 (MQ=255)
                                                        aTTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAG‑GAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTAc                                                                          <  2:438225/151‑1 (MQ=255)
                                                        aTTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAG‑GAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTAc                                                                          >  2:268832/1‑151 (MQ=255)
                                                        aTTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAG‑GAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTAc                                                                          <  1:44910/151‑1 (MQ=255)
                                                              ggCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAG‑GAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAAtt                                                                     <  2:733305/150‑1 (MQ=255)
                                                                    atGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAG‑GAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCAc                                                              >  2:672033/1‑151 (MQ=255)
                                                                         gatgatAATGAATTTGAATCAGGGATGAAAAATCAAG‑GAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTc                                                          <  1:596837/150‑1 (MQ=255)
                                                                         gatgatAATGAATTTGAATCAGGGATGAAAAATCAAG‑GAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCATCTACCAATTCCGCCACCtt                                                           <  1:412839/149‑1 (MQ=255)
                                                                            gatAATGAATTTGAATCAGGGATGAAAAATCAAG‑GAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcac                                                      >  2:803890/1‑151 (MQ=255)
                                                                            gatAATGAATTTGAATCAGGGATGAAAAATCAAG‑GAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGca                                                       <  1:460757/150‑1 (MQ=255)
                                                                            gatAATGAACTTGAATCGGGGAGGAAGAACAGAG‑GACTTACCAATAAAGGAATCAAAGAGAA‑TGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCtt                                                           <  1:107699/146‑1 (MQ=11)
                                                                             atAATGAATTTGAATCAGGGATGAAAAATCAAG‑GAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcac                                                      >  2:637329/1‑150 (MQ=255)
                                                                               aaTGAATTTGAATCAGGGATGAAAAATCAAG‑GAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAg                                                    <  1:742095/150‑1 (MQ=255)
                                                                                    aTTTGAATCAGGGATGAAAAATCAAG‑GAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCt                                              <  2:649274/151‑1 (MQ=255)
                                                                                    aTTTGAATCAGGGATGAAAAATCAAG‑GAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCt                                              <  2:429932/151‑1 (MQ=255)
                                                                                       tGAATCAGGGATGAAAAATCAAG‑GAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTa                                            <  1:122161/150‑1 (MQ=255)
                                                                                        gAATCAGGGATGAAAAATCAAG‑GAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTGATCtt                                             >  2:127891/1‑148 (MQ=255)
                                                                                         aaTCAGGGATGAAAAATCAAG‑GAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACtt                                         >  2:151726/1‑151 (MQ=255)
                                                                                         aaTCAGGGATGAAAAATCAAG‑GAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACtt                                         >  1:20887/1‑151 (MQ=255)
                                                                                          aTCAGGGATGAAAAATCAAG‑GAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACttt                                        >  2:481102/1‑151 (MQ=255)
                                                                                            cAGGGATGAAAAATCAAG‑GAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACttttt                                      <  2:116812/151‑1 (MQ=255)
                                                                                            cAGGGATGAAAAATCAAG‑GAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACttttt                                      >  2:346958/1‑151 (MQ=255)
                                                                                            cAGGGATGAAAAATCAAG‑GAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACttttt                                      <  1:505680/151‑1 (MQ=255)
                                                                                            cAGGGATGAAAAATCAAG‑GAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACttt                                        >  1:637491/1‑149 (MQ=255)
                                                                                                 aTGAAAAATCAAG‑GAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATATTACTTTTTTTGa                                  >  2:114534/1‑150 (MQ=255)
                                                                                                 aTGAAAAATCAAG‑GAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGCCATCTTACTTTTATTGa                                  >  1:465766/1‑150 (MQ=255)
                                                                                                   gAAAAATCAAG‑GAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGatat                               <  1:729064/151‑1 (MQ=255)
                                                                                                    aaaaaTCAAG‑GAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATc                              <  2:703979/151‑1 (MQ=255)
                                                                                                     aaaaTCAAG‑GAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGACGGTGGGTACTTGAATCCCCACGTCCGTAAGCACACTAACACCTGACGCTATCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGata                                <  1:13054/149‑1 (MQ=25)
                                                                                                     aaaaTCAAG‑GAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCg                             >  1:825608/1‑151 (MQ=255)
                                                                                                      aaaTCAAG‑GAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGatat                               >  1:363987/1‑148 (MQ=255)
                                                                                                          cAAG‑GAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCTCAATACtttttaagatttcctatc                         >  1:168820/1‑137 (MQ=18)
                                                                                                          cAAG‑GAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTc                         <  1:757846/150‑1 (MQ=255)
                                                                                                                 agaaACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGt                  <  2:617780/151‑1 (MQ=255)
                                                                                                                     aCAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCg               >  2:266228/1‑150 (MQ=255)
                                                                                                                       aaGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAg             >  2:473115/1‑150 (MQ=255)
                                                                                                                         gAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAgggg          >  2:731581/1‑151 (MQ=255)
                                                                                                                          aaaGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAggggg         <  2:530878/151‑1 (MQ=255)
                                                                                                                           aaGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAgggggg        <  1:305473/151‑1 (MQ=255)
                                                                                                                            aGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTATTGATATCGCCTCGCTTGGTGCGAggcggg        >  1:559579/1‑150 (MQ=255)
                                                                                                                                 gTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACtt  >  2:33578/1‑151 (MQ=255)
                                                                                                                                  tAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACtt  >  2:266080/1‑150 (MQ=255)
                                                                                                                                    aaGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTTCTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACtt  >  2:589428/1‑148 (MQ=255)
                                                                                                                                    aaGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTATCCTTACTTTAAATGAAATCCTATCGCTGGGGGCGAGTCGGGGGACtt  >  1:555518/1‑148 (MQ=14)
                                                                                                                                                    |                                                                                                                                    
TGCTGTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAG‑GAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTT  >  NC_000913/4605940‑4606219

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